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d ' Abstract 

' A diagonalizable non-Hermitian Hamiltonian having a real spectrum may be used to define a unitary quantum 

system, if one modifies the inner product of the Hilbert space properly. We give a comprehensive and essentially 
JO ■ self-contained review of the basic ideas and techniques responsible for the recent developments in this subject. We 

provide a critical assessment of the role of the geometry of the Hilbert space in conventional quantum mechanics 
i | , to reveal the basic physical principle motivating our study. We then offer a survey of the necessary mathematical 

' tools, present their utility in establishing a lucid and precise formulation of a unitary quantum theory based 

on a non-Hermitian Hamiltonian, and elaborate on a number of relevant issues of fundamental importance. In 
. particular, we discuss the role of the antilinear symmetries such as VT, the true meaning and significance of the 

so-called charge operators C and the CPT-inner products, the nature of the physical observables, the equivalent 
description of such models using ordinary Hermitian quantum mechanics, the pertaining duality between local- 
|"~ | non-Hermitian versus nonlocal-Hermitian descriptions of their dynamics, the corresponding classical systems, the 

i i ' pseudo-Hermitian canonical quantization scheme, various methods of calculating the (pseudo-) metric operators, 

subtleties of dealing with time-dependent quasi-Hermitian Hamiltonians and the path-integral formulation of the 
C*") ' theory, and the structure of the state space and its ramifications for the quantum Brachistochrone problem. 
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PACS number: 03.65.-w, 03.65.Ca, ll.30.-j 

Keywords: pseudo-Hermitian, quasi-Hermitian, metric operator, unitary-equivalence, biorthonormal system, VT- 
■ symmetry, projective space, complex potential, quantization, observable 

> : 

X 1 Contents 

b : 

- - i 1 Introduction and Overview 3 

2 Mathematical Tools and Conceptual Foundations 7 

2.1 Hilbert Spaces and Riesz Bases 7 

2.2 Bounded, Invertible, and Hermitian Operators 9 

2.3 Unitary Operators and Unitary-Equivalence 11 

2.4 Biorthonormal Systems 13 

2.5 Metric Operators and Conventional QM 14 

3 Pseudo-Hermitian QM: Ingredients and Formalism 16 

3.1 Quasi-Hermitian versus Pseudo-Hermitian QM 16 

3.2 Pseudo-Hermitian and Pseudo-Metric Operators 17 

3.3 Spectral Properties of Pseudo-Hermitian Operators 22 

3.4 Symmetries of Pseudo-Hermitian Hamiltonians 23 

3.5 A Two-Level Toy Model 26 



1 



4 Calculation of Metric Operator 29 

4.1 Spectral Method 29 

4.1.1 PT-symmetric infinite square well 29 

4.1.2 'PT-symmetric barrier 29 

4.1.3 Delta-function potential with a complex coupling 30 

4.1.4 Other Models 30 

4.2 Perturbation Theory 31 

4.2.1 'PT-symmetric cubic anharmonic oscillator 32 

4.2.2 Imaginary cubic potential 33 

4.2.3 Other Models 33 

4.3 Differential Representations of Pseudo-Hermiticity 33 

4.3.1 Field equation for the metric from Moyal product 34 

4.3.2 Universal field equation for the metric 34 

4.4 Lie Algebraic Method 35 

4.4.1 Lie algebras and their representations 36 

4.4.2 General outline of the method 36 

4.4.3 Swanson Model: Q = su(l, 1) 38 

5 Systems Defined on a Complex Contour 39 

5.1 Spectral problems defined on a contour 39 

5.2 Equivalent spectral problems defined on R 41 

5.3 Wrong-sign quartic potential 43 

6 Complex Classical Mechanics versus Pseudo-Hermitian QM 46 

6.1 Classical-Quantum Correspondence and Observables 46 

6.2 Complex Classical Systems and Compatible Poisson Brackets 47 

6.3 Real Description of a Complex Classical System 49 

7 Time-Dependent Hamiltonians and Path-Integral Formulation 51 

7.1 Time-dependent Quasi- Hermitian Hamiltonians 51 

7.2 Path-Integral Formulation of Pseudo-Hermitian QM 51 

8 Geometry of the State-Space and the Quantum Brachistochrone 54 

8.1 State-Space and Its geometry in Conventional QM 54 

8.2 State-Space and Its geometry in Pseudo-Hermitian QM 57 

8.3 Optimal-Speed Evolutions 59 

9 Physical Applications 62 

9.1 Earlier Applications 62 

9.1.1 Dyson Boson Mapping 62 

9.1.2 Complex Scaling and Resonances 63 

9.1.3 Vortex Pinning in Superconductors 63 

9.2 Relativistic QM, Quantum Cosmology, and QFT 63 

9.3 Electromagnetic Wave Propagation 65 

9.4 Other Applications and Physical Manifestations 67 

9.4.1 Atomic Physics and Quantum Optics 67 

9.4.2 Open Quantum Systems 67 

9.4.3 Magnetohydrodynamics 67 

9.4.4 Quantum Chaos and Statistical Mechanics 68 

9.4.5 Biophysics 68 

10 Summary and Conclusions 68 
Appendix 69 



2 



1 Introduction and Overview 



General Relativity (GR) and Quantum Mechanics (QM) are the most important achievements of the twentieth century 
theoretical physics. Their discovery has had an enormous impact on our understanding of Nature. Ironically, these 
two pillars of modern physics are incompatible both conceptually and practically. This has made their unification 
into a more general physical theory the most fundamental problem of modern theoretical physics. The unification of 
Special Relativity and QM, which is by far an easier task, has been the subject of intensive research since late 1920's. 
It has led to the formulation of various quantum field theories. A most successful example is the Standard Model 
which provides a satisfactory description of all available observational data in high energy particle physics. In spite 
of the immensity of the amount of research activity on the subject and the fact that this is conducted by the most 
capable theoretical physicists of our times, the attempts at quantizing gravity have not been as successful. In fact, 
one can claim with confidence that these attempts have so far failed to produce a physical theory offering concrete 
experimentally verifiable predictions. This state of affairs has, over the years, motivated various generalizations of 
GR and QM. Although none of these generalizations could be developed into a consistent physical theory capable of 
replacing GR or QM, the hope that they might facilitate the discovery of a unified theory of quantum gravity still 
motivates research in this direction. 

The development of the special relativistic quantum theories has also involved attempts at generalizing QM. 
Among these is an idea initially put forward by Dirac in 1942 [77] and developed by Pauli [190] into what came to 
be known as the indefinite-metric quantum theories [223, 183, 185]. This is a rather conservative generalization of 
QM in which one considers in addition to the physical states of the system a set of hypothetical states, called ghosts, 
whose function is mainly to improve the regularity properties of the mathematical description of the physical model. 
The indefinite-metric quantum theories lost their interest by mid 1970's and perhaps unfortunately never found a 
detailed coverage in standard textbooks on relativistic quantum mechanics and quantum field theory. 1 

A more recent attempt at generalizing QM is due to Bender and his collaborators [37, 38] who adopted all its 
axioms except the one that restricted the Hamiltonian to be Hcrmitian. They replaced the latter condition with the 
requirement that the Hamiltonian must have an exact PT-symmetry. Here P and T are the parity and time-reversal 
operators whose action on position wave functions ip{ x ) is given by (P^>)(x) '■= an d (T0)( x ) := ^{x)* 2 The 

exact PT-symmetry of a Hamiltonian operator H means that it has a complete set of PT- invariant eigenvectors ip n , 
i.e., VTipn = Obn^n for some complex numbers a n . This condition assures the reality of the spectrum of H. A class 
of thoroughly studied examples is provided by the PT-symmetric Hamiltonians of the form 

H v = \ P 2 -{ixy, (i) 

where v is a real number not less than 2, and the eigenvalue problem for H v is defined using an appropriate contour 
r in the complex plane C which for v < 4 may be taken as the real line R. A correct choice for V assures that 
the spectrum of H is discrete, real, and positive, [32, 33, 79, 209, 157]. Another example with identical spectral 
properties is the PT-symmetric cubic anharmonic oscillator, 

H=^p 2 + ^ 2 x 2 +zex 3 , (2) 

whose coupling constants fi and e are real and its eigenvalue problem is defined along the real axis (T = R), [40, 158]. 

This PT-symmetric modification of QM leads to an indefinite- metric quantum theory [116, 225, 160], if one 
endows the Hilbert space with the indefinite inner product, 

(•!•>„ :=<•!?>•>, (3) 

known as the PT-inner product [37]. The symbol (-|-) that appears in (3) stands for the L 2 -inner product: (<f>\ip) := 
J r (j>{z)*tj}{z)dz that defines the Hilbert space L 2 (T) of square-integrable functions ip : T — > C, where T is the contour 
in complex plane that specifies the PT-symmetric model [157]. 

The main advantage of the indefinite inner product (3) over the positive-definite inner product (-|-) is that the 
former is invariant under the time-evolution generated by the Schrodinger equation [240, 143, 155], i.e., if (f>(t) and 
ifj(t) are solutions of the Schrodinger equation for the PT-symmetric Hamiltonian H, (^>(t)\t/j(t)) v does not depend 
on time. 

In order to employ the standard formulation of indefinite metric quantum theories [183] for a PT-symmetric 
model we proceed as follows [116]. 

1 In fact, Pauli who had made fundamental contributions to the very foundations of the subject had developed strong critical views 
against it. For example in his Nobel Lecture of 1946, he identifies a "correct theory" with one that does not involve "a hypothetical 
world", [191]. For a comprehensive critical assessment of indefinite-metric quantum field theories, sec [185]. 

2 We use asterisk to denote complex-conjugation. 
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1. We split the space H of state-vectors into the subspaces %± := {ip £ H | sgn((ip\tp) v ) = ±}, where sgn(a) := 
a/\a\ if a is a nonzero real number and sgn(O) := 0. T-L± are orthogonal subspaces in the sense that for all 

i'± €H±, i^-\i>+)-P =«- 3 

2. We impose a superselection rule that forbids interactions mixing the elements of H- and T~L+ and try to devise 
a solution for the difficult problem of providing a physical interpretation of the theory [183, 185]. The simplest 
way of dealing with this problem is to identify the elements of T~L+ with physical state- vectors [230, 110] and 
effectively discard the rest of state-vectors as representing unphysical or ghost states. 

An alternative formulation of the theory that avoids the interpretational difficulties of indchnite-metric theories 
is the one based on the construction of a genuine positive-definite inner product on H.. This construction was initially 
obtained in [144] as a by-product of an attempt to derive a necessary and sufficient condition for the reality of the 
spectrum of a general Hamiltonian operator H that possesses a complete set of eigenvectors [143, 144, 145]. Under 
the assumption that the spectrum of H is discrete, one can show that it is real if and only if there is a positive- 
definite inner product (-|-)+ that makes it Hermitian, i.e., (-\H-)+ = (H ■ |-) + , [144, 145]. The proof of this statement 
involves an explicit construction of (-|-) + . This inner product depends on the choice of the Hamiltonian through its 
eigenvectors. Hence it is not unique [146, 148, 93, 162]. 

In [37] the authors propose a different approach to the problem of identifying an appropriate inner product for 
the 'PT-symmetric Hamiltonians such as (1). They introduce a generic symmetry of these Hamiltonians which they 
term as C-symmetry and construct a class of positive-definite inner products, called the CVT -inner products, that, 
as we show in Section 3.4, turn out to coincide with the inner products (-|-) + , [148, 157]. The approach of [37] may 
be related to a much older construction originally proposed in the context of the indefinite-metric quantum theories 
[183, 185]. It involves the following two steps. 

1. Suppose that T~L = H+ ffi H- where H± are the orthogonal subspaces we defined above, and that both H± 
admit a basis consisting of the eigenvectors of H. 

2. Let n : W — > H be the projection operator onto H+, so that for all ip £ H, ip+ '■= £ H+ and ip_ := 
ip-Tp+£H-. Clearly U 2 =11 and Jhj}_ = 0. 

3. Endow H with the positive-definite inner product: (4>,ip) ■= (4>+\ip+) v — ((f>-\ip-) v . 

4. Let C : H —> H. be defined by Cip := ip+ — ip-- Then, in view of the orthogonality of H±, 

(^V) = («) P =(C^) V , (4) 

It is not difficult to see that C = 211 — /, where I stands for the identity operator acting in "H, i.e., the operator 
that leaves all the elements of T~L unchanged. Obviously for all ip £ H±, Cip = ±ip. Hence, C is a grading operator 
associated with the direct sum decomposition % = "H+ ffi %- of %. Furthermore, in view of the assumption 1 
above, the eigenvectors of H have a definite grading. This is equivalent to the condition that C commutes with the 
Hamiltonian operator, [C, H] = 0, [8]. It turns out that the CPT-inner product introduced in [37] coincides with the 
inner product (4). 

In [37, 38], the authors use the CPT-inner product to formulate a unitary quantum theory based on the VT- 
symmctric Hamiltonians (1). They identify the observables O of the theory with the CPT-symmetric operators , in 
particular 

CTTOCVT=0. (5) 

This definition is motivated by the demand that the structure of the theory must not involve mathematical operations 
such as Hermitian conjugation and be determined only using physical conditions. The definition (5) does fulfil this 
demand 5 , but it suffers from a serious dynamical inconsistency in the sense that in the Heisenberg picture an operator 
0(t) := e ltH O(0)e~ ltH that commutes with CVT at t = may not commutes with this operator at t > 0. Therefore, 
in general, under time-evolution an observable can become unobservable [152, 153]! This inconsistency rules out (5) 

3 The assumption of the existence of such an orthogonal decomposition is referred to as "dccomposability" of the model in indefinite- 
metric theories, [183]. For 'PT-symmetric systems considered in the literature, this assumption is valid if there is a complete basis of 
common eigenvectors of the Hamiltonian and VT . 

4 To ensure that the spectrum of such a CPT-symmetric operator O is real, one demands that O has an exact CPT-symmetry, i.e., its 
eigenstates are left invariant under the action of CVT. This does not however ensure the completeness of the eigenvectors of O. 

5 The authors of [37] emphasize this point by stating that: "In effect, we replace the mathematical condition of Hermiticity, whose 
physical content is somewhat remote and obscure, by the physical condition of space-time and charge-conjugation symmetry." One must 
note however that in this theory the usual coordinate operator x no longer represents a physical observable. As a result V does not affect 
a reflection in the physical space, and there is no reason why one should still refer to PT'-symmetry as the physical "space-time reflection 
symmetry" as done in [37, 38]. 
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as an acceptable definition of a physical observable. As noticed in [39, 42], it can be avoided, if one replaces the 
symmetry condition (5) with: 

CVTOCTT=O t , (6) 

where all operators are identified with their matrix representation in the coordinate-basis and O t stands for the 
transpose of O. In particular, (x\O t \x') := (x'\0\x). 

The presence of the mathematical operation of transposition in (6) shows that apparently the theory could 
not be defined just using "the physical condition of space-time and charge-conjugation symmetry" as was initially 
envisaged, [37, 38]. Note also that (6) puts an implicit and difficult-to-justify restriction on the Hamiltonian H. 
Being an observable commuting with CVT, H must satisfy H T = H, i.e., it is necessarily symmetric! 6 Therefore (6) 
cannot be used to determine the observables of a theory that has a nonsymmetric Hamiltonian. The restriction to 
symmetric Hamiltonians may be easily lifted, if one is willing to adopt the conventional definition of the observables, 
namely identifying them with the operators that arc Hermitian with respect to the CPT-inner product [152, 154]. 

(;0-) = (0-,-), (7) 

Indeed this definition is forced upon us by a well-known mathematical theorem that we present a detailed proof 
of in the appendix. It states that if a linear operator O has real expectation values computed using a given inner 
product, then O is necessarily Hermitian with respect to this inner product. This shows that the requirement of the 
Hcrmiticity of observables and in particular the Hamiltonian has a simple "physical" justification. 7 

An important motivation for considering this so-called PT -symmetric Quantum Mechanics is provided by an 
interesting idea that is rooted in special rclativistic local quantum field theories (QFT). Among the most celebrated 
results of QFT is the CPT-theorem. It states that every field theory satisfying the axioms of QFT is CPT-invariant, 
[98], where C is the charge-conjugation operator. Clearly replacing the axiom that the Hamiltonian is Hermitian 
with the statement of the CPT-theorem might lead to a generalization of QFT. The implementation of this idea 
in a nonrelativistic setting corresponds to the requirement that the Hamiltonian possesses an exact PT-symmetry. 
This in turn allows for the construction of a C-operator that similarly to the charge-conjugation operator C of QFT 
squares to identity and generates a symmetry of the system. The idea that this establishes a nonrelativistic analog 
of the CPT-theorem is quite tempting. Yet there are clear indications that this is not really the case. For example, 
unlike the charge-conjugation operator of QFT, C depends on the choice of the Hamiltonian. It turns out that in 
fact this operator does not play the role of the relativistic charge-conjugation operator, it is merely a useful grading 
operator for the Hilbert space. 8 In this sense the adopted terminology is rather unfortunate. 

One of the aims of the present article is to show that PT-symmetric QM is an example of a more general class of 
theories, called Pseudo-Hermitian Quantum Mechanics, in which PT-symmetry does not play a basic role and one 
does not need to introduce a C-operator to make the theory well-defined. The analogs of PT and C operators can 
nevertheless be defined in general [148], but they do not play a fundamental role. All that is needed is to determine 
the class of non-Hermitian Hamiltonians that are capable of generating a unitary evolution and a procedure that 
associates to each member of this class a positive-definite inner product that renders it Hermitian. It turns out that 
there are always an infinite class of positive-definite inner products satisfying this condition. Each of these defines 
a separate physical Hilbert space with a complete set of observables. In this way one obtains a set of quantum 
systems that have the same Hamiltonian but different Hilbert spaces. Therefore, they are dynamically equivalent 
but kincmatically distinct. 

In order to elucidate the conceptual foundations of Pseudo-Hermitian QM wc will next examine some of the 
basic properties of the mathematical notions of the "transpose" and "Hermitian-conjugatc" of a linear operator. For 
clarity of the presentation we will only consider the operators that act in the space of square-integrable functions 
L 2 (R). The discussion may be generalized to square-integrable functions defined on a complex contour [157]. 

In the literature on PT-symmetric QM, notably [37, 38, 42], the transpose O t and Hcrmitian-conjugate of a 
linear operator O are respectively defined with respect to the coordinate-basis, {\x)}, according to 

(xi|0 T |x 2 ) := (x 2 \0\xi), (8) 
(si|Ot|x a ) := {x 2 \0\x 1 )\ (9) 

where x\,x 2 arc arbitrary real numbers. Therefore the terms "symmetric" and "Hermitian" respectively refer to 
the conditions (xi|0|x2) = (x2|0|xi) and (xi|0|x2)* = (x2|0|xi). These definitions reflect the inclination to treat 

6 In mathematical literature the term "symmetric operator" is usually used for a different purpose, as discussed in footnote 19 below. 
To avoid possible confusion we will not adopt this terminology. 

7 What is however not dictated by this requirement is the choice of the inner product. 

8 One can more appropriately compare C with the chirality operator (7 s ) for the Dirac spinors. 
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operators as matrices. This is certainly permissible provided that one specifies the particular basis one uses for this 
purpose. In this sense the following equivalent definitions are more preferable. 



O 1 := J dan J dx 2 (x 2 \0\ Xl ) \xi){x 2 \, T := / dxi J dx 2 (x 2 \0\ Xl )* \x 1 )(x 2 \. (10) 

A nice feature of (9) that is not shared with (8) is that it is invariant under the basis transformations that map 
{\x}} onto any orthonormal basis. For example one can easily show that if (9) holds, so do 

(pi\Ct\p2) = fa\0\pi)* and Ot = / 'd Pl J ' dp 2 {p 2 \0\ Pl y\vi){p 2 \. (11) 

This invariance under orthonormal basis transformations stems from the fact that O** admits a basis-independent 
definition: It is the linear operator satisfying 

MOty) = (O^), (12) 

i.e., the adjoint operator for O. 9 

The notions of "transpose" and "symmetric operator" introduced above and employed in [37, 38, 39, 42] do 
not share this invariance property of "Hermitian-conjugate" and "Hermitian operator" . For example, it is easy to 
see that (xi\(ip)\x 2 ) = — (x 2 \(ip)\xi) while (pi\(ip)\p 2 ) = (p 2 \(ip)\pi) . Therefore, ip is represented by a symmetric 
matrix in the p-basis while it is represented by an antisymmetric matrix in the x-basis. This shows that there is no 
basis- independent notion of the transpose of an operator or a symmetric operator. 10 

Obviously once we specify a basis, there is no danger of using definition (8). But we must keep in mind that any 
theory in which one uses the notion of transposition in the sense of (8) involves the implicit assumption that the 
coordinate-basis is a preferred basis. The use of the notion of Hcrmitian-conjugation as defined by (9) does not rely 
on such an assumption. As we will explain in the following section the choice of an orthonormal basis is equivalent 
to the choice of an inner product. This is why one can define using its basis-independent defining relation (12) 
which only involves the inner product (-|-). In summary, while the use of the terms "transpose" and "symmetric 
operator" involves making a particular choice for a preferred basis, the use of the term "Hermitian-conjugate" and 
"Hermitian operator" involves making a particular choice for an inner product. 

In conventional QM the inner product is fixed from the outset. Hence the notions of Hermitian-conjugation and 
Hermitian operator are well-defined. The opposite is true about the notions of transposition and symmetric operator. 
This does not cause any difficulty, because they never enter into quantum mechanical calculations, and in principle 
one does not need to introduce them at all. We will see that the same is the case in Pseudo-Hermitian QM. In 
particular, in the discussion of 'PT-symmetric systems, there is no need to identify physical observables using (6). 

The main reason for making a universal and preassigned choice for the inner product in QM is the curious fact 
that up to unitary-equivalence there is a unique inner product. 11 This means that using different inner products 
leads to physically identical theories, or more correctly to different representations of a single physical theory. In 
conventional QM, one eliminates the chance of employing these alternative representations by adopting the usual 
(L 2 -) inner product as the only viable choice. The situation resembles a gauge theory in which one fixes a gauge from 
the outset and then forgets about the gauge symmetry. This will have no effect on the physical quantities computed 
using such a theory, but it is clearly not recommended. It is quite possible that an alternative choice of gauge would 
facilitate a particular calculation. 

We wish to argue that because no one has ever made an independent measurement of the inner product of the 
Hilbert space, it must be kept as a degree of freedom of the formulation of the theory. This is the basic principle 
underlying Pseudo-Hermitian QM. 

We will see that any inner product may be defined in terms of a certain linear operator 77+. It is this so-called 
metric operator that determines the kinematics of pscudo-Hcrmitian quantum systems. The Hamiltonian operator 
H that defines the dynamics is linked to the metric operator via the pscudo-Hermiticity relation, 

fft = 77+U77+ 1 . (13) 

We will explore some of the consequences of this equation whose significance has not been fully noticed or appreciated 
in its earlier investigations, notably in the context of the indefinite-metric quantum theories [190]. 



9 A rigorous definition of the adjoint operator is given in Subsection 2.2. 

10 One can define a notion of the transpose of a linear operator O acting in a complex inner product space V in terms of an arbitrary 
antilinear involution r : V — > V according to O t = tO^t, [6]. A function r : V — ► V is called an antilinear operator if T(atp + b<f>) = 
a*Ttj) + b*T(f> for all complex numbers a, b and all elements t/j, (f> of V. It is called an involution if t 2 = I, where / : V — > V is the identity 
operator. The choice of a basis to define O t is equivalent to the choice of an antilinear involution r. The notion of the transpose used 
in [37, 38, 42] corresponds to choosing the time-reversal operator T as r. 

xl This will be explained in detail in Subsection 2.3. 
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We wish to point out that there is a very large number of publications on the subject of this article. Many of these 
focus on the mathematical issues related to the investigation of the spectrum of various non-Hermitian operators or 
formalisms developed to study such problems. Here we will not deal with these issues. The interested reader may 
consult the review article [80]. Another related line of research is the mathematical theory of linear spaces with 
an indefinite metric and its applications [53, 18]. This is also beyond the aim and the scope of the present article. 
There are a series of review articles by Bender and his collaborators [38, 42, 28, 29, 30] that also address the physical 
aspects of the subject. The approach pursued in these articles is based on the use of the CPT-inner products and the 
definition of observables given by (6). This restricts the domain of validity of their results to symmetric Hamiltonians. 
The discussion of the classical limit of 'PT-symmetric systems offered in these articles is also not satisfactory, because 
it does not involve a quantization scheme that would relate classical and quantum observables. 

It is our opinion that to gain a basic understanding of the subject demands a careful study of the underlying 
mathematical structures without getting trapped in the physically irrelevant mathematical details and technicalities. 
The need for a comprehensive and readable treatment of basic mathematical notions and their physical consequences 
has not been met by any of the previously published reviews. In the first part of the present article (Sections 2 
and 3), we intend to address this need. Here we only discuss the mathematical tools and results that are necessary 
for addressing the conceptual issues of direct relevance to the physical aspects of the subject. In section 3, wc 
use the mathematical machinery developed in Section 2 to present a complete formulation of pscudo-Hcrmitian 
QM and its connection with VT- and C-symmetries. The second part of the article (Sections 4-9) aims to survey 
various recent developments. In Section 4 wc survey different methods of computing metric operators. In Sections 
5-8 we explore systems defined on complex contours, the classical limit of pscudo-Hermitian quantum systems, the 
subtleties involving time-dependent Hamiltonians and the path-integral formulation of the theory, and the quantum 
Brachistochrone problem, respectively. In Section 9 we discuss some of the physical applications of pseudo-Hermitian 
Hamiltonians, and in Section 10 we present our concluding remarks. 

2 Mathematical Tools and Conceptual Foundations 

In this section we survey the necessary mathematical tools and elaborate on a number of conceptual issues that 
are helpful in clarifying various existing misconceptions on the subject. We also offer a through discussion of the 
motivation for considering a more general formulation of QM. 

One of the axioms of QM is that pure physical states of a quantum system are rays in a Hilbert space H. Each 
ray may be determined in a unique manner by a nonzero element ip of "H which we call a state-vector. The physical 
quantities associated with a pure state is computed using a corresponding state- vector and the inner product of the 
Hilbert space. We begin our discussion by a precise description of Hilbert spaces, inner products, bases, Hermitian 
and unitary operators, biorthonormal systems, and their relevance to our investigation. 

We will use the following notations and conventions: R, R + , C, Z, Z + , N denote the sets of real numbers, positive 
real numbers, complex numbers, integers, positive integers, and nonnegative integers (natural numbers), respectively. 
The symbol " := " means that the left-hand side is defined to be the right-hand side, " =: " means that the 
converse is true, and " € " and " C " respectively stand for "is an element of" and "is a subset of" . Throughout 
this paper we will only consider time-independent Hamiltonian operators unless otherwise is explicitly stated. 

2.1 Hilbert Spaces and Riesz Bases 

Consider a complex vector space V and a function (-|-) : V x V — > C that assigns to any pair ip, <p of elements of V a 
complex number (ip\(p). Suppose that (-|-) has the following properties. 

(i) It is positive-definite, i.e., for all nonzero elements ip of V, (ip\ip) is a positive real number, and it vanishes 
if and only if ip = 0, where we use also for the zero vector. 

(it) It is Hermitian, i.e., for any pair ip, <p of elements of V, (ip\<p)* = ((p\ip)- 

(Hi) It is linear in its second slot, i.e., for all ip, <p, x € V and all a, b € C, (ip\a<p + bx) = a(ip\<f>) + k(V'lx)- 

Then (-|-) is called an inner product 12 on V, and the pair (V, (•]■)) is called an inner product space. 

An inner product (-|-) on V assigns to each element ip of V a nonnegative real number, || ip ||:= y/ (ip\ip), that 
is called the norm of ip. We can use the norm to define a notion of distance between elements of V, according to 
|| ip — <p || j and develop analysis and geometry on the inner product space (V, (-|-)). 

A Hilbert space % is an inner product space which fulfills an additional technical condition, namely that its 
norm defines a complete metric space, i.e., for any infinite sequence {ipk} of elements tpk of T-L, the condition that 

12 We use the terms "inner product" and "positive-definite inner product" synonymously. 
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linij.fc-s.oo || ipj — 4>k ||= implies that {V'fc} converges to an element ijj of %: limi._5.0c, || V ; — 11= 0. In other words, 
a Hilbert space is a complete inner product space. 

A subset S of a Hilbert space H is said to be dense, if every element of T-L may be obtained as the limit of a 
sequence of elements of S. A Hilbert space is said to be separable, if it has a countable dense subset. It turns out 
that T-L is separable if and only if it has a countable basis. The latter is a sequence {xn} of linearly independent 
elements of Ti such that the set of its finite linear combinations, 

£({xn}) = jE c «X" \KeZ+, c„ e C j , (14) 

is a dense subset of Ti. For an infinite-dimensional Hilbert space Ti, {Xn} is an infinite sequence and the assertion 
that C{{xn}) is a dense subset means that every element tj) of Ti is the limit of a convergent series of the form 
Y^Li c nXn whose coefficient c n are assumed to be uniquely determined by ip- 13 

It is not difficult to show that any finite-dimensional inner product space is both complete and separable. In this 
sense infinite-dimensional separable Hilbert spaces are natural generalizations of the finite-dimensional inner product 
spaces. In the following we will use the label N to denote the dimension of the Hilbert space in question. N = 00 
will refer to an infinite-dimensional separable Hilbert space. 

An important difference between finite and infinite-dimensional Hilbert spaces is that the definition of a basis in 
a finite-dimensional Hilbert space does not involve the inner product while the opposite is true about the infinite- 
dimension Hilbert spaces. The requirement that (14) be a dense subset makes explicit use of the norm. Therefore 
whether a given sequence of linearly independent vectors forms a basis of an infinite-dimensional Hilbert space Ti 
depends in a crucial manner on the inner product of T-L. 

Given a basis {xn} of a separable Hilbert space Ti, one can apply the Gram-Schmidt process [198, 66] to construct 
an orthonormal basis, i.e., a basis {£ n } satisfying 

(£ml£n) = <5 m „ for all m,ne {1,2,3," ■ ,N}, (15) 

where 5 mn denotes the Kronecker delta symbol: 8 mn :— if m 7^ n and S nn := 1 for all n. For an orthonormal basis 
{Cn}, the coefficients c„ of the basis expansion, 

N 

i> = Y] c„£„, (16) 



of the elements tpofH are given by 
Furthermore, in view of (15) - (17), 



N 

W#) = ^<<^n)<£„|V) for all M e H. (18) 

71=1 



In particular, 

Ikl-ElW)! 3 forall^eH. (19) 



N 



Eq. (18) implies that the operator X defined by lip := X)ri=i(£™l^) equals the identity operator / acting in 
%; I = J. 14 This is called the completeness relation which in Dirac's bra-ket notation takes the familiar form: 

ElllU(Cn|=/- 

Next, we wish to examine if given a basis {Cn} of a separable Hilbert space "H, there is an inner product (-|-) on 
% with respect to which {Cn} is orthonormal. Because {Cn} is a basis, for all ip, <f> € T-L there are unique complex 
numbers c n ,d n such that ip = Y] ^ c„Cn an d 4> — 5_ n =i d n Cn- We will attempt to determine {4>\4>) in terms of c n 
and d n . 



13 This notion of basis is sometimes called a Schauder basis [236]. It is not to be confused with the algebraic or Hamcl basis which for 
infinite-dimensional Hilbert spaces is always uncountable [99]. 

14 Let ipi,ip2 S H be such that for all <j> S H, = (</>\ip2)- Then {<f>\^)i — »/_} = for all <j> S H. Setting <j> = tpi — ip2, one then 

finds || ipi — ip2 || 2 = which implies ipi = ip2- I n view of this argument, {(fr\Tip) = {(fr\ip) for all ip, (f> S H implies Xifi = ip for all tp S H. 
Hence X = I. 
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First, consider the finite-dimensional case, N < oo. Then, in view of Eq. (18), the condition that {Cn} is 
orthonormal with respect to (-|-) defines the latter according to 



We can easily check that (-|-) possesses the defining properties (i) - (Hi) of an inner product and satisfies (Cm|Cn) = 
5 mn . It is also clear from Eq. (18) that any other inner product with this property must satisfy (20). This shows 
that (-I-) is the unique inner product that renders {Cn} orthonormal. 

The case N = oo may be similarly treated, but in general there is no guarantee that the right-hand side of 
(20) is a convergent series. In fact, it is not difficult to construct examples for which it is divergent. Therefore, an 
inner product that makes an arbitrary basis {Cn} orthonormal may not exist. The necessary and sufficient condition 
for the existence of such an inner product (-|-) is that ip = Y^Li c n(n implies Y^Li \ c n\ 2 < oo for all ip € H. 15 
Furthermore, we shall demand that the inner product space H' obtained by endowing the underlying vector space of 
H with the inner product (-|-) is a Hilbert space. As we will discuss in Subsection 2.3, any two infinite-dimensional 
separable Hilbert spaces, in particular H and H' , have the same topological properties, i.e., the set of open subsets 
of H coincide with those of H'. This restricts (-|-) to be topologically equivalent to (-|-), i.e., there are positive real 
numbers C\ and Ci satisfying ci(ip\ip) < (ip\ip) ^ c 2(ip\ip) for a U ip € H. It turns out that the inner product (20) 
that renders the basis {Cn} orthonormal and is topologically equivalent to (-|-) exists and is unique if and only if it 
is obtained from an orthonormal basis {£ n } through the action of an everywhere-defined bounded invcrtiblc linear 
operator A : H — > H, i.e., Cn = A£ n . A basis {Cn} having this property is called a Riesz basis, [94, 236]. In summary, 
we can construct a new separable Hilbert space H' in which {Cn} is orthonormal if and only if it is a Riesz basis. We 
will give a derivation of a more explicit form of the inner product of H' in Subsection 2.4. 

2.2 Bounded, Invertible, and Hermitian Operators 

Consider two Hilbert spaces Hi and Hi with inner products and(-|-) 2 , respectively, and a linear operator A that 
maps Hi to Hi- The domain T>(A) of A is the subset of Hi such that the action of A on any element of T>(A) yields 
a unique element of Hi- The range IZ(A) of A is the subset of Hi consisting of elements of the form Aipi where ipi 
belongs to V(A). If V(A) = Hi, one says that A has full domain, or that it is everywhere- defined. If 1Z(A) = Hi, 
one says that A has full range, i.e., it is an onto function. As an example consider the momentum operator p acting 
in Hi = Hi = L 2 (R), (pip)(x) := —ih-fj^ip(x). Then V(p) consists of the square-intcgrable functions that have a 
square- integrable derivative, and 7Z(p) is the set of square- integrable functions ipi such that ipi(x) = S-oo ipi(u)du is 
also square- integrable. In particular p is not everywhere-defined, but its domain is a dense subset of L 2 (M). Such an 
operator is said to be densely- defined. All the operators we encounter in this article and more generally in QM are 
densely-defined. 

A linear operator A : Hi — > Hi is said to be bounded if there is a positive real number M such that for all 
ip € V (A), || Atp \\i< M || ijj ||i, where || ■ ||i and || • \\i are respectively the norms defined by the inner products (-1^ 
and (-|-) 2 . The smallest M satisfying this inequality is called the norm of A and denoted by || A \\. A characteristic 
feature of a bounded operator is that all its eigenvalues a are bounded by its norm, \a\ <\\ A \\. Furthermore, a linear 
operator is bounded if and only if it is continuous [195]. 16 Linear operators relating finite-dimensional Hilbert spaces 
are necessarily bounded. Therefore the concept of boundedness is only important for infinite-dimensional Hilbert 
spaces. 

A : Hi — > Hi is called an invertible operator if it satisfies both of the following conditions [105]. 17 

1. A is one-to-one and onto, so A -1 : Hi Hi exists and has a full domain; 

2. A" 1 is a bounded operator. 

If A is bounded, one-to-one and onto, then according to a theorem due to Banach its inverse is also bounded; it is 
invertible, [127]. 18 An important class of bounded invertible operators that play a fundamental role in QM is the 
unitary operators. We will examine them in Subsection 2.3. 

15 If (-|-) exists, (20) must hold because {Cn} is orthonormal with respect to (-|-)- Hence the left-hand side of (20) is well-defined and its 
right-hand must be convergent. In particular, for if> = ip we find that (ip\ip) = JZil'Li l c «| 2 must be finite. This establishes the necessity 
of the above condition. Its sufficiency follows from the inequality: For all K £ Z+, | c n^»i| 2 ^ E^Li \ c n\ 2 Mm| 2 < oo. 

16 A function / : Hi — > Hz is said to be continuous if for all ip £ f(/) and every sequence {ip n } in i>(f) that converges to ip, the 
sequence {/(V^i)} converges to f(ip). 

17 Some authors do not require the second condition. The above definition is the most convenient for our purposes. 

18 A continuous one-to-one onto function with a continuous inverse is called a homeomorphism. The domain and range of a homeomor- 
phism have the same topological properties, li f : Hi — > H2 is & homeomorphism relating two Hilbert spaces Hi and Hi, a, sequence 
{in} of elements x n of H converges to an element x £ H if and only if the {/(x n )} converges to f(x) in H2- 



N 




(20) 



n=l 
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Next, consider a linear operator A : Hi — > H2 that has a dense domain T>(A). Let V be the subset of H2 
whose elements ip2 satisfy the following condition: For all ipi £ D (A), there is an element <f>\ of Hi such that 
(?p2\Aipi) 2 = {(f>i\^pi) 1 . Then there is a unique linear operator A^ : H 2 — > Hi fulfilling [195] 

(^ 1 \A^ 2 ) 1 = (Mi\^2) 2 for all 6 V{A) and i> 2 £ V . (21) 

This operator is called the adjoint or Hermitian- conjugate of A. By construction, V is the domain of A'; D(j4T) = T>' . 

For the case H 2 = Hi =: H, where H is a Hilbert space with inner product (-|-), a linear operator A : H — > % 
having a dense domain 2?(^4) is called a self-adjoint or Hermitian operator 19 if = A. In particular, 2?(A^) = 2? (A) 
and 

(ih\Aih) = {A^ 2 ) for all Vi,V>2 tV(A). (22) 

An occasionally useful property of Hermitian operators is that every Hermitian operator having a full domain is 
necessarily bounded. This is known as the Helling er-Toeplitz theorem [195]. 

Hermitian operators play an essential role in QM mainly because of their spectral properties [229]. In particular, 
their spectrum 20 is real, their eigenvectors with distinct eigenvalues are orthogonal, and they yield a spectral resolution 
of the identity operator I. For a Hermitian operator A with a discrete spectrum the latter takes the following familiar 
form, if we use the Dirac bra-ket notation. 

N 

I=J2\»n)(a n \, (23) 

71=1 

where {a n } is an orthonormal basis consisting of the eigenvectors a n of A whose eigenvalues a n are not necessarily 
distinct, 

Aa n = a n a n , for all n £ {1, 2, 3, • • • , N}. (24) 
Eqs. (23) and (24) imply that A is diagonalizable and admits the following spectral representation 

N 

A = ^2a„\a„)(a„\. (25) 

n=l 

Well-known analogs of (23) and (25) exist for the cases that the spectrum is not discrete, [123]. 

An important property that makes Hermitian operators indispensable in QM is the fact that for a given densely- 
defined 21 linear operator A with T> (A) = D(A^), the expectation value (i/j\Ai[>} is real- valued for all unit state- vectors 
ip £ T>(A) if and only if the Hermiticity condition (22) holds, [123]. 22 This shows that in a quantum theory 
that respects von Neumann's measurement (projection) axiom, the observables cannot be chosen from among non- 
Hcrmitian operators even if they have a real spectrum. The same conclusion may be reached by realizing that 
the measurement axiom also requires the eigenvectors of an observable with distinct eigenvalues to be orthogonal, 
for otherwise the reading of a measuring device that is to be identified with an eigenvalue of the observable will 
not be sufficient to determine the state of the system immediately after the measurement [175]. This is because 
an eigenvector that is the output of the measurement may have a nonzero component along an eigenvector with a 
different eigenvalue. This yields nonzero probabilities for the system to be in two different physical states, though 
one measures the eigenvalue of one of them only! 

Let {£ ra } be an orthonormal basis of H and A be a Hermitian operator acting in H, then according to property 
(ii) of the inner product and the Hermiticity condition (22) we have 

Aran '■= (S,m\A^ n ) = (£n|^4£ m )* = A* nm . (26) 

This shows that A is represented in the basis {£, n }, according to 

N 

A = J2 A mn\U)^n\, (27) 
n=l 

using the N x N Hermitian matrix 23 A := (A mn ). 



19 In some mathematics texts, e.g., [195], the term "Hermitian" is used for a more general class of operators which satisfy (22) but have 
T>(A) C D(A^). A more commonly used term for such an operator is "symmetric operator". We will avoid using this terminology which 
conflicts with the terminology used in the literature on PT-symmetric QM, [37, 38, 42]. 

20 For a precise definition of the spectrum of a linear operator see Subsection 3.3. 

21 Observables must have dense domains, for otherwise one can construct a state in which an observable cannot be measured! 
22 A simple proof of this statement is given in the appendix. 
23 A square matrix M is called Hermitian if its entries M mn satisfy M m 
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It is essential to realize that a Hermitian operator can be represented by a non-Hermitian matrix in a non- 
orthonormal basis. This implies that having the expression for the matrix representation of an operator and knowing 
the basis used for this representation are not sufficient to decide if the operator is Hermitian. One must in addition 
know the inner product and be able to determine if the basis is orthonormal. Referring to an operator as being 
Hermitian or non-Hermitian (using its matrix representation) without paying attention to the inner product of the 
space it acts in is a dangerous practice. 

For example, it is not difficult to check that the Hermitian matrix (7% = f J * J represents the operator 

L : C 2 —> C 2 defined by J ■— ( Zi z 2 z? J m the basis B := {( J j < ( \ j\- The same operator is represented 

in the standard basis £> :={ ( J ) '( 1 )} usm S tne non-Hermitian matrix -l ) ' ^ e w * s ^ to stress tnat this 

information is, in fact, not sufficient to ascertain if L is a Hermitian operator, unless we fix the inner product on the 
Hilbert space C 2 . For instance, if we choose the standard Euclidean inner product which is equivalent to requiring 
Bq to be orthonormal, then L is a non-Hermitian operator. If we choose the inner product that makes the basis B 
orthonormal, then L is a Hermitian operator. We can use (20) to construct the latter inner product. It has the form 
(z\w) := z*(wi - W2) + z^i—wi + 2w 2 ), where z = f V and w = ( ™* J. 

The above example raises the following natural question. Given a linear operator H that is not represented by a 
Hermitian matrix in an orthonormal basis, is there another (non-orthonormal) basis in which it is represented by a 
Hermitian matrix? This is equivalent to asking if one can modify the inner product so that H becomes Hermitian. 
The answer to this question is: "No, in general." As we will see in the sequel, there is a simple necessary and sufficient 
condition on H that ensures the existence of such an inner product. 

2.3 Unitary Operators and Unitary-Equivalence 

Equations (17) and (18) may be employed to derive one of the essential structural properties of separable Hilbert 
spaces, namely that up to unitary-equivalence they are uniquely determined by their dimension. To achieve this we 
first explain how one compares inner product spaces. Two inner product spaces Hi and H2 with inner products 
and (-|-) 2 are said to be unitary- equivalent, if there is an everywhere-defined onto linear operator U : Hi — > H2 such 
that for every </>i, m Hi we have 

(C7^i|C/Vi) 2 = (^IV'i) 1 - (28) 
Such an operator is called a unitary operator. In view of (12) and (28), we have 

tfU = h, (29) 

where Ii denotes the identity operator acting on Hi. One can use (28) and the ontoness property of U to show that 
U is an invertible operator and its inverse that equals U* , is also unitary. 

Unitary-equivalence is an equivalence relation. 2 Therefore to establish the unitary-equivalence of all iV-dimcnsional 
separable Hilbert spaces it suffices to show that all of them are unitary-equivalent to a chosen one. The most conve- 
nient choice for the latter is the Hilbert space 



niN _ 
H — 



: N for N ± 00 
t 2 for N = 00, 



where C N is the set of TV-dimensional complex column vectors endowed with the standard Euclidean inner prod- 
uct (w\z) := w* ■ z, a dot denotes the usual dot product, I 2 is the set of square-summable sequences, £ 2 := 
{ {c n } I c„ e C , Yln=i M 2 < 00 }> equipped with the inner product: ( {c„} | {c„} ) := J2^Li 5 ,*Cn, and {c„}, {c„} G 
£ 2 . Now, let H be any ^-dimensional separable Hilbert space with inner product (-|-), {£ T i} be an orthonormal basis of 
H, and U : H — >• Hq be defined by U(ip) := { (£ n \'4>) } for all ip in H. It is not difficult to sec that U is an everywhere- 
defined and onto linear map. Furthermore in view of (18) it satisfies, (U(j)\U'ip) = X^=i(£n|</')*(£n|V') = (0|VOj f° r 
all <f), tp £ H. Hence, (28) holds, U is a unitary operator, and H is unitary-equivalence to Hq . 

For a quantum system having R as its configuration space, one usually uses the coordinate wave functions ip( x ) to 
represent the state- vectors. The latter are elements of the Hilbert space L 2 (E) of square- integrable complex- valued 
functions ip : R — > C. The inner product of i 2 (R) has the form 

poo 

(<^) = / dx 4>{x)*^{x). (30) 



24 This means that every inner product space is unitary-equivalent to itself; liHi is unitary-equivalent to H2, so is H2 to Hi; if "Hi 
unitary-equivalent to H2 and H2 is unitary-equivalent to "Hz, then Hi is unitary-equivalent to Hs- 



11 



A concrete example for an orthonormal basis {£ n } for L 2 (M.) is the basis consisting of the standard normalized 
eigenfunctions £„ = ip n -i of the unit simple harmonic oscillator Hamiltonian [139], or equivalently that of the 
operator —£j + x 2 , i.e., ip m (x) = ^^(mW^^e-* 2 / 2 H m (x), where H m (x) = e* 2 ' 2 (x- -^) m e' x2 ' 2 are Hermite 
polynomials and m £ N. For this example, switching from the coordinate-representation of the state-vectors to their 
representation in terms of the energy eigenbasis of the above harmonic oscillator corresponds to affecting the unitary 
operator U. 

Unitary operators have a number of important properties that follow from (28). For example if U : Hi — > H2 is 
a unitary operator relating two separable Hilbcrt spaces Hi and H2, then 



1. U is an everywhere-defined, bounded, and invertible operator. 

2. If {£n} is an orthonormal basis of Hi, then {[/£„} is an orthonormal basis of H.2- 



25 



3. Let A\ : Hi — > Hi be a Hermitian operator with domain T>(Ai). Then UAiW : H2 — > H2 is a Hermitian 
operator with domain U(V(Ai)) = {Z7-01 G H 2 \ipi £ T>(A{)}. 

A direct consequence of statement 1 above and the fact that the inverse of a unitary operator is unitary is that 
unitary operators are homeomorphisms 26 . As a result, all iV-dimcnsional separable Hilbcrt spaces have identical 
topological properties. In particular, if two separable Hilbert spaces Hi and H2 share an underlying vector space, a 
sequence {x n } converges to a; in Hi if and only if it converges to 1 in W2. 

Next, we recall that every quantum system s is uniquely determined by a separable Hilbcrt space H that determines 
the kinematic structure of s and a Hamiltonian operator H : H —> H that gives its dynamical structure via the 
Schrodinger equation. Let si and S2 be quantum systems corresponding to the Hilbert spaces Hi, H2 and the 
Hamiltonians Hi : Hi — > Hi, H2 : H2 — > Hi- By definition the observables of the system Sj, with i £ {1,2}, are 
Hermitian operators Oi : Hi — > Hi- si and S2 are physically equivalent, if there is a one-to-one correspondence 
between their states and observables in such a way that the physical quantities associated with the corresponding 
states and observables are identical. Such a one-to-one correspondence is mediated by a unitary operator U : Hi —> H2 
according to 

^^.^--UiPi Oi -> 2 := UOiUK (31) 

In particular, if ipi(t) is an evolving state- vector of the system Si, i.e., it is a solution of the Schrodinger equation 
ih4rtpi(t) = Hiipiit), we have ip2(t) = Utpi(t). The necessary and sufficient condition for the latter is H2 = U H-JJ^ . 
This observation motivates the following theorem. 

Theorem 1: As physical systems si = S2, if there is a unitary operator U : Hi — > H2 satisfying 
H 2 = U HiW. 

To prove this assertion we recall that all physical quantities in QM may be expressed as expectation values of 
observables. Suppose that such a unitary operator exists, and let us prepare a state of si that is represented by 
V>i G Hi and measure the observable Oi. The expectation value of this measurement is given by 

(VilV-i)i (uH 2 \Wip 2 )i (^2^2)2 ' 1 ' 

where is the inner product of Hi, and we have used (31) and the fact that U' is also a unitary operator. Eqs. (32) 
show that the above measurement is identical with measuring O2 in a state of S2 represented by the state- vector i\>2- 
This argument is also valid for the case that ipi is an evolving state-vector. It shows that the existence of U implies 
the physical equivalence of si and S2. 27 

If the hypothesis of Theorem 1 holds, i.e., there is a unitary operator U : Hi —> H2 satisfying H 2 = U HiW , we 
say that the pairs (Hi, Hi) and (^2,^2) and the quantum systems they define arc unitary- equivalent. 

In conventional QM one mostly considers unitary operators that act in a single Hilbcrt space H. These generate 
linear transformations that leave the inner product of the state- vectors invariant. They form the unitary group 14(H) 
of the Hilbert space which includes the time-evolution operator e r ttH / h as a one-parameter subgroup. 28 Given a 
quantum system s with Hamiltonian H, one may use the unitary operators U £ U(H) to generate unconventional 
unitary-equivalent systems s v having the same Hilbert space. If x and p are the standard position and momentum 
operators of s, the position, momentum, and Hamiltonian operators for s v are respectively given by x u := U x U', 
p v := U pU\ and H v := U H . The transformation x — > x u , p — > p u , and H — > H u is the quantum analog of a 
classical time- independent canonical transformation. Therefore, unitary transformations generated by the elements 
oiU(H) play the role of canonical transformations. 

25 The converse is also true in the sense that given an orthonormal basis {£' H } of H2 there is a unique unitary operator U : Hi — > H2 
such that i' n = U£ n for all n £ {1, 2, • • • , N}. 
26 See footnote 18 for the definition. 

27 The converse of Theorem 1 can be formulated similarly to the Wigner's symmetry theorem [231, p91]. 
28 All possible symmetry, dynamical, and kinematic groups arc also subgroups of U(H), [142]. 
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2.4 Biorthonormal Systems 

Let {ipn} be a basis of an TV-dimensional separable Hilbert space "H, with N < oo, and be the orthonormal 
basis obtained by performing the Gram-Schmidt process on {ip n }- Because {ip n } is a basis, there are unique complex 
numbers B mn £ C such that for all m £ {1, 2, 3, • • • , A} 

N 



and because {£ n } is an orthonormal basis, 



n=l 



N 



k=l 

for all n £ {1,2,3, ••■ ,A}. If we substitute (34) in (33) and use the orthonormality of £„, we find Yl n =i Bnm(£j\ipn) = 
S m j for all m,j £ {1, 2, 3, • • • , A}. This shows that the A x A matrix B = (B mn ) is invcrtible and the entries of B _1 
are given by = (£ m \' , Pn}- We can use this relation to express (34) in the form 

N 

^n = 5>^&- (35) 
fe=l 

This equation suggests that the vectors cp m defined by 

N 

4> m :=Y,B* mj ^ for all me {1,2,3,- ,N}, (36) 

3=1 

fulfil 

(<p m \ipn) = S mn for all m,n £ {1, 2, 3,-- • ,N}. (37) 

Furthermore, employing (16), (35), and (36), we can show that for all ip £%, '52, n - 1 (<j>\ip)ipn = ip- We can use Dirac's 
bra-ket notation to express this identity as 

N 

J2^n){ct>n\=I- (38) 

n=l 

This is a generalization of the more familiar completeness relation (23). 

A sequence {(ipn, fin)} of ordered pairs of elements of "H that satisfy (37) is called a biorthonormal system 
[201, 205, 236]. A biorthonormal system satisfying (38) is said to be complete. 

Let {ipn} be a basis of a separable Hilbert space H, and {4> n } be a sequence in H such that {(ip n , <Pn)} is a 
complete biorthonormal system. Then, one can show that {<p n } is the unique sequence with this property and that it 
is necessarily a basis of %, [236]. {</> n } is called the biorthonormal basis associated with {ip n }, and the biorthonormal 
system {(^„,0„)} is called a biorthonormal extension of {ip n }- 

If A < oo, the right-hand side of (36) is well-defined and the above construction yields the biorthonormal basis 
{</>„} associated with every basis {ip n } of H. UN — oo, {<p m } can be constructed provided that the right-hand side 
of (36) converges. This is the case if and only if 

oo 

^T\B mj \ 2 <oo forallmeZ+. (39) 

3=1 

A theorem due to Bari [94] states that given a basis {ip n }, the biorthonormal system {(tp n , <p n )} exists and J2n°=i I (^n 1^)1 
and J2n°=i K^"!^)! 2 both converge for all ip £ H, ii and only if {ipn} is a Riesz basis, i.e., there are an orthonormal 
basis {Xn} and an everywhere-defined bounded invertible operator A : H — > % such that ip n = Ax n - 29 In this case 
{4> n } is also a Riesz basis, and {(ip n , 4>n)} is the unique biorthonormal extension of {ip n }, [94, 201, 236]. 

The coefficients of the expansion of a vector ip £ H in a Riesz basis {ip n } can be expressed in terms of its 
biorthonormal basis {<p n }- Given ip = X^=i c nipn we have, in light of (37), c n = (<p n \ip). Hence, for all ip £ H, 

N 



2 



n=l 



29 As any pair of orthonormal bases are related by a unitary operator U : H — > H which is an everywhere-defined bounded invertible 
operator, one can take \ n = £ n without loss of generality. This allows for the identification of the infinite matrix B with the matrix 
representation of A -1 in the basis {£n}, for we have B mn = (^m|A _1 g n ). 
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Clearly the roles of {4> n } and {ip n } are interchangeable. In particular, X) n =i l&^iipn] = I and tp = X^^i^nlV') 0" 
for all V € "H. 

Let {V'n} be a Riesz basis and {(i/vi, </>«)} be its biorthonormal extension. As we explained in Subsection 2.1, we 
can construct a unique inner product (-|-) on H that makes a Riesz basis orthonormal. We can use (20) and (40) to 
obtain the following simplified expression for this inner product. 

N 

(V#) := X^I^X^) = M»7+$> for all ^,0 e^, (41) 

ra=l 

where we have introduced the operator rj + : H — > H according to 

N 

7 ?+ V:=^(0n|V')0n- (42) 
n=l 

Using Dirac's bra-ket notation we can write it in the form 

N 

V+ = J2\<i>n)(<f>n\. (43) 

71 = 1 

Because {i[>n} is a Riesz basis, the inner product (41) is defined for all ip, 4> € H. This shows that rj + is everywhere- 
defined. Furthermore, it is not difficult to see, by virtue of (37), that it has an everywhere-defined inverse given 

by 

N 

^ X := X>nX^I- (44) 

n=l 

This shows that 77+ is a one-to-one onto linear operator. As suggested by (43) it is also Hermitian, which in particular 
implies that both 77+ and 77+ 1 are bounded. Finally, in view of (43), (tp\r]+'ip) = J2 n =i K^nl^)] 2 , for all ip e H. 
Therefore r/ + is a positive operator. Moreover because it is an invcrtiblc operator, its spectrum is strictly positive. 
A Hermitian operator with this property is called a positive-definite operator. The operator i] + constructed above 
is an everywhere-defined, bounded, positive-definite, invcrtiblc operator. A linear operator with these properties is 
called a metric operator. 



2.5 Metric Operators and Conventional QM 

Consider a pair of separable Hilbcrt spaces Hi and H2 that are identical as vector spaces but have different inner 
products. We will denote the inner products of Hi and Hi by and (-|-) 2 , respectively, and view (-|-) 2 as an 

alternative inner product on Hi. Our aim is to find a way to express (-|-) 2 in terms of (-1^. We will first consider 
the case that the underlying vector space V of both Hi and H2 is finite-dimensional, i.e., N < 00. 

Let {£„} be an orthonormal basis of Hi . As we argued above, it satisfies the completeness relation: Yl n =i \% n ) 1 I = 
/, where / is the identity operator of V . In general, is not an orthonormal basis of H2 and the operator 

N 

V+--=52\tn) 22 (U (45) 

n=l 

does not coincide with I. Eq. (45), which we can express in the more conventional form: 

N 

77+ := X^IVX £«, for all ^ e V, (46) 



71 = 1 



defines 77+ as a linear operator 77+ : V — > V having a full domain. 

Now, let <p be an arbitrary element of V. In view of (16), we can express it as 



N 



Eftml^&n. (47) 



771=1 



Using Eqs. (46), (47) and the properties (it) and (Hi) shared by both the inner products, we can easily show that 
77+ fulfils 

(</#) 2 = (#7+ i>), for all tA, G V. (48) 
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Employing properties (i) and (ii) of (-|-)2, we can also verify that 

(#7+^ = <#7+0>i for all </;,</> eV, (49) 
(^1??+ ip) 1 > for all nonzero ip £ V. (50) 

Eq. (49) shows that 77+ is a Hermitian operator acting in Hi, in particular it has a real spectrum. Eq. (50) implies 
that indeed the spectrum of 77+ is strictly positive; 77+ is a positive-definite operator. 

If Hi and H2 are infinite-dimensional, the sum appearing in (46) stands for an infinite series and we must address 
the issue of its convergence. The convergence of this series is equivalent to the requirement that K&ilV^I < 00 
for all ip £ V. According to the above-mentioned theorem of Bari this requirement is fulfilled provided that {£, n } is a 
Ricsz basis of H2- Under this assumption (46) defines a linear operator 77+ acting in "Hi and satisfying (48). It is a 
positive-definite (in particular Hermitian) operator with a full domain. Being Hermitian and everywhere-defined it 
is also necessarily bounded. 

Next, we exchange the roles of Hi and H2- This gives rise to an everywhere- defined bounded positive-definite 
operator r)' + acting in H2 such that (<p\ip) 1 = (</ , \v'+ f° r a ^ V'i 4> in V. Combining this equation and (48) yields 
(<p\ip) 1 = (<j)\ri-^rf, ip) j for all <p in V. This in turn implies that 77+77+ = /. 30 Because 77+ is a bounded operator 
with a full domain, 77+ is an invertible operator with inverse 77^ 1 = ?;+. 31 Therefore 77+ is a metric operator. 

The above construction of the metric operator 77+ is based on the choice of an orthonormal basis {£ n } of Hi. Wc 
will next show that indeed 77+ is independent of this choice. Let 77+ : V — > V be an everywhere-defined linear operator 
satisfying 

(M>> 3 = (#7W>i for all 7/;, 0e V. (51) 

Eqs. (48) and (51) show that (4>\(r]' + — V+)' l P)i = f° r an ip,<p £V. In view of the argument presented in footnote 14, 
this implies rj' + ip = 77+7/1 for all ip £ V, i.e., 77^ = 77+. This establishes the uniqueness of the metric operator 77+ 
which in turn means that the inner products that make a complex vector space V into a separable Hilbert space are 
in one-to-one correspondence with the metric operators 77+ acting in one of these Hilbert spaces. 

To employ the characterization of the inner products in terms of metric operators we need to select a Hilbert 
space H. We will call this Hilbert space a reference Hilbert space. We will always fix a reference Hilbert space and 
use its inner product (-|-) to determine if a given linear operator acting in V is Hermitian or not. The following are 
some typical examples. 

• For systems having a finite number (N) of linearly independent state- vectors, i.e., when V — C N , H is defined 
by the Euclidean inner product, (<p\ip) := <p* ■ ip, for all <p,ip £ C . 

• For systems whose configuration space is a diffcrcntiable manifold M with an integral measure d/j,(x), V is 
the space L 2 (M) of all square- integrable functions ip : M — > C and H is defined by the L 2 -inner product: 
(<p\ip) := J M (p(x)*ip(x) dfi(x), for all <p,ip £ L 2 (M). The systems whose configuration space is a Euclidean 
space (M = R d ) or a complex contour (M = V) belong to this class. We will discuss the latter systems in 
Section 5. 

In summary, given a separable Hilbert space H with inner product (-|-), we can characterize every other inner 
product on H by a metric operator 77+ : H — > H according to 

(•|->„ + ■=(-\v+-)- (52) 

Each choice of 77+ defines a unique separable Hilbert space H v+ . Because 77+ is a positive-definite operator, we can 
use its spectral representation to construct its positive square root p := \pr\,- As a linear operator acting in H, p is 
a Hermitian operator satisfying p 2 = 7/+ . We can use this observation to establish 

{cP\ip) n+ = {,p\r 1+ iP) = {p^\pip) = {p<t>\pi>), faraU^GU (53) 

This relation shows that as a linear operator mapping H onto H, p is a unitary operator. 32 It provides a realization 
of the unitary-equivalence of the Hilbert spaces H v+ and H. 

In ordinary QM one fixes the physical Hilbert space of the system to be one of the reference Hilbert spaces listed 
above and develops a theory based on this preassigned Hilbert space. The argument that the unitary-equivalence of 

30 The proof uses the argument given in footnote 14. 

31 If we use the prescription of Subsection 2.1 to obtain the inner product (-|-) on "H% that renders the Riesz basis {£1} orthonormal, 
we find by its uniqueness that (-|-) = (•]•)!• According to (41), there is an everywhere-defined bounded invertible operator 77+ such that 

=(-|t?+->2- :t is § iven b y V+ = V+ 1 - 

32 Strictly speaking (53) shows that p is an isometry. However, in view of the fact that 77+ is invertible, so does p. This implies that 
p : — > H is a genuine unitary operator. 
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all separable Hilbert spaces justifies this convention is not quite acceptable. For example, although for all d G Z + , 
L 2 (R d ) is unitary-equivalent to L 2 (R), we never use L 2 (M.) to describe a system having more than one real degree 
of freedom. The choice of the particular Hilbert space should in principle be determined by physical considerations 
or left as a freedom of the formulation of the theory. In view of lack of a direct measurement of the inner product 
or the associated metric operator, we propose to adopt the second option. We will see that this does not lead to 
a genuine generalization of QM, but it reveals a set of alternative and equally useful representations of QM which 
could not be utilized within its conventional formulation. Furthermore, the introduction of the freedom in choosing 
the metric operator may be used as an interesting method of extending QM by relaxing some of the restrictions put 
on the metric operator. Indeed the indefinite-metric quantum theories are examples of such a generalization. 

3 Pseudo-Hermitian QM: Ingredients and Formalism 

3.1 Quasi-Hermitian versus Pseudo-Hermitian QM 

To the best of our knowledge, the first publication investigating the consequences of the freedom in the choice of 
the metric operator is the article by Scholtz, Geycr, and Hahnc [207] in which the choice of the metric operator is 
linked to that of an irreducible set of linear operators. The latter is any (minimal) set S of operators O a acting in 
a vector space V that do not leave any proper subspace of V invariant, i.e., the only subspacc V' of V that satisfies: 
" O a e S and ip G V imply O a ip G V, " is V. 

The approach pursued in [207] involves using the physical characteristics of a given system to determine an 
irreducible set of operators (that are to be identified with the observables of the theory) and employing the latter to fix 
a metric operator and the associated inner product of the Hilbert space. We will call this formalism Quasi-Hermitian 
Quantum Mechanics. The main problem with this formalism is that it is generally very difficult to implement. This 
stems from the fact that the operators belonging to an irreducible set must in addition be compatible, i.e., there 
must exist an inner product with respect to which all the members of the set are Hermitian. In order to use this 
formalism to determine the inner product, one must in general employ a complicated iterative scheme. 

• Select a linear operator 0\ with a complete set of eigenvectors and a real spectrum; 

• Find the set ill of all possible metric operators that make 0\ Hermitian, and select a linear operator O2, 
linearly independent of 0\, from among the linear operators that are Hermitian with respect to the inner 
product defined by some r] + G Hi ; 

• Find the set H2 of all possible metric operators that make O2 Hermitian, and select a linear operator O3, 
linearly independent of 0\ and O2, from among the linear operators that are Hermitian with respect to the 
inner product defined by some n+ G ill nil2, where D stands for the intersection of sets; 

• Repeat this procedure until the inner product (respectively metric operator r)+) is fixed up to a constant 



As we see, in trying to employ Quasi-Hermitian QM, one needs a procedure to compute the most general metric 
operator associated with a given diagonalizable linear operator with a real spectrum. This is the main technical tool 
developed within the framework of Pseudo-Hermitian QM. 

Pseudo-Hermitian QM differs from Quasi-Hermitian QM in that in the former one chooses 0\ to be the Hamilto- 
nian, finds ill and leaves the choice of O2 arbitrary, i.e., identifies all the operators O that arc Hermitian with respect 
the inner product associated with some unspecified metric operator 77+ belonging to iti . The metric operator 77+ fixes 
a particular inner product and defines the "physical Hilbert space" "Hphys of the system. The "physical observables" 
are the Hermitian operators O acting in H. p hys- We can use the unitary-equivalence of "Hphys and H realized by 
p := \/v + to construct the physical observables O using those of the conventional QM, i.e., Hermitian operators o 
acting in the reference Hilbert space H. This is done according to [152, 175] 



Note that as an operators mapping H p h ys onto H, p is a unitary operator. Hence O is a Hermitian operator acting 
in "Hphys if and only if o is a Hermitian operator acting in TL. 

For instance, let H = L 2 (R d ) for some d G Z + . Then we can select the usual position xt and momentum pi 
operators to substitute for o in (54). This defines a set of basic physical observables, 



coefficient. 



O = p^op. 



(54) 



X t := p 1 x l p, 



Pi ■= P 1 PiP, 



(55) 
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which we respectively call the n+ -pseudo-Hermitian position and momentum operators. They furnish an irreducible 
unitary representation of the Weyl-Heisenberg algebra, 



[X,, X,} = [P u Pj) = 0, [Xi, Pj] = HSijI, for all i,j e {1, 2, ■ • • , d}. (56) 

In principle, we can express the Hamiltonian H as a function of Xi and Pi and attempt to associate a physical 
meaning to it by devising a quantum-to-classical correspondence principle. One way of doing this is to define the 
underlying classical Hamiltonian H c for the system as 

H c (x Cl p c ) := lim H(X,P) , (57) 

n— ±o x ^ x c 

? -> PC 

where w := (iui, u>2, ■ • ■ ,Wd) T for to = X,P,x c ,p c ; and x c and p c stand for classical position and momentum 
variables. Supposing that the right-hand side of (57) exists, one may reproduce the quantum system described by 
the Hilbcrt space "H p hy S arid Hamiltonian H by quantizing the classical system corresponding to H c according to 

X C ^X, Pc^P, {;-} c ^-ih- 1 [;-], (58) 

where {■, -} c and [-, ■] stand for the Poisson bracket and the commutator, respectively. This is called the rj + -pseudo- 
Hermitian canonical quantization scheme [175, 158]. 

The quantum system described by 'Hphys and H admits a representation in conventional QM in which the Hilbcrt 
space is H = L 2 (R d ), the observablcs arc Hermitian operators acting in %. and the Hamiltonian is given by 

h := pHp- 1 . (59) 

Because p : "H p hy S - ► H is unitary, so is its inverse p^ 1 : TL — > "Hphys- This in turn implies that h is a Hermitian 
operator acting in T~i. We will call the representation of the quantum system that is based on the Hermitian 
Hamiltonian h the Hermitian representation. In this representation, we can proceed employing the usual prescription 
for identifying the underlying classical Hamiltonian, namely as 

H c {x c ,Pc) ■— \imh(x,p) . (60) 

h— >o x -> x c 

V — ► Pc 

Note that this relation is consistent with (57), because in view of (55) and (59), h — f(x,p) if and only if H = f(X, P), 
where we suppose that / : R 2d — s- C is a piecewise real-analytic function. 

Each choice of a metric operator 7j + € iti corresponds to a particular quantum system with an associated 
Hermitian representation. One can in principle confine his (her) attention to this representation which can be fully 
understood using the conventional QM. The main disadvantage of employing the Hermitian representation is that in 
general the Hamiltonian h is a terribly complicated nonlocal operator. Therefore, the computation of the energy levels 
and the description of the dynamics are more conveniently carried out in the pseudo-Hermitian representation. In 
contrast, it is the Hermitian representation that facilitates the computation of the expectation values of the physical 
position and momentum operators as well as that of the localized states in physical position or momentum spaces. 
See [175, 158, 159, 162] for explicit examples. 



3.2 Pseudo-Hermitian and Pseudo-Metric Operators 

Definition 1: A linear operator A-.T-L^H acting in a separable Hilbert space 1-L is said to be pseudo- 
Hermitian if T>(A) is a dense subset of %, and there is an everywhere- defined invertible Hermitian 
linear operator rj : H — s- H such that 

A^=riArf x . (61) 



We will refer to an operator n satisfying (61) as a pseudo-metric operator associated with the operator A and denote 
the set of all such operators by DJIa- Clearly, A is pseudo-Hermitian if and only if DJIa is nonempty. Furthermore, 
for every linear operator A, 9JIa C 9Jt/ where / is the identity operator. We will call elements of 971/ pseudo-metric 
operators. Because they are Hermitian and have full domain, they are necessarily bounded. 33 

33 The definition of a pseudo-Hermitian operator given in [143] requires i; to be a Hermitian automorphism. This is equivalent to 
the definition given above because of the following. Firstly, an automorphism is, by definition, everywhere-defined. Hence, if it is 
Hermitian, it must be bounded. Secondly because the inverse of every Hermitian automorphism is a Hermitian automorphism, r; -1 is 
everywhere-defined and bounded, i.e., r\ is invertible. The fact that every everywhere-defined invertible operator is an automorphism is 
obvious. 
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Clearly if rji belongs to 971,4, then so does rj r := rrji, for every nonzero real number r. The scaling 771 — > r\ r of 
the pseudo-metric operators has no physical significance. It signifies a spurious symmetry that we will eliminate by 
restricting our attention to pseudo-metric operators 77 whose spectrum o~(A) is bounded above by 1, i.e., max[cr(j4)] = 
l. 34 The latter form a subset of 971 a which we will denote by 51 a- 

In general, either 11,4 is the empty set and A is not pseudo-Hcrmitian or Ua is an infinite set 35 ; the pseudo-metric 
operator associated with a pseudo-Hermitian operator is not unique. To see this, let r\ € Ha, B : % — > H be an 
everywhere-defined invertible bounded operator commuting with A, and 77 := B^r\B. Then is an everywhere- 
defined bounded operator [235] commuting with A' . These in turn imply that 77 is an everywhere-defined invertible 
Hcrmitian operator which in view of (61) satisfies rjArj^ 1 = Bh]BAB~ 1 i]~ 1 B~ 1 ^ = B^A^B^ 1 ^ = A\ i.e., 77 € 971,4. 
Clearly there is an infinity of choices for B defining an infinite set of pseudo- metric operators of the form B^r/B. It 
is not difficult to observe that upon making appropriate scaling of these pseudo-metric operators one can construct 
an infinite set of elements of Ha, he., Ha is an infinite set. 

The non-uniqueness of pseudo-metric operators associated with a pseudo-Hermitian operator motivates the fol- 
lowing definition. 

Definition 2: Let A be a pseudo-Hermitian operator acting in % and 77 G 971/ be a pseudo-metric 
operator. Then A is said to be rj-pseudo-Hermitian if rj € 9JIa ■ 

Clearly, in order to determine whether a pseudo-Hermitian operator A is 77-pseudo-Hermitian one must know both 
A and 77. It is quite possible that a pseudo-Hermitian operator fails to be 77-pseudo-Hermitian for a given 77 G 971/. 36 
Given a pseudo-Hermitian operator A, the set 971,4 may or may not include a positive-definite element 77+ . If such 
a positive-definite clement exists, we can use it to construct the inner product 

(■IV :=<-kr) (62) 

that renders A Hermitian, 

(.\A.) rl+ =(A.\-) v+ . (63) 

This means that if we endow the underlying vector space V of H with the inner product (62), we find a separable 
Hilbert space H ri such that A : H — > H is Hermitian. In particular the spectrum a(A) of A is real. If cr(A) 
happens to be discrete, we can construct an orthonormal basis {ipn} of H r consisting of the eigenvectors of A. As 
a sequence of elements of H, {ip n } is a Riesz basis. Hence A : H — > T~L is diagonalizable. This shows that for a 
densely-defined operator having a discrete spectrum, the condition that it is a diagonalizable operator having a real 
spectrum is necessary for the existence of a metric operator 77+ among the elements of 971,4 . The existence of 77+, in 
particular, implies that 971,4 is non-empty. Hence, A is necessarily pseudo-Hcrmitian. 

It is not difficult to show that the same conditions arc also sufficient for the inclusion of a metric operator in 971a , 
[144, 145]. Suppose that A : H — s- W is a densely-defined diagonalizable operator having a real spectrum. Let {ipn} 
be a Riesz basis consisting of the eigenvectors of A and {(V>n, 4>n)} be its biortho normal extension. Then, in view of 
the spectral representation of A, i.e., A = Y] ^ o n \ip n ) {4>n\ where a n are eigenvalues of A, and the basic properties 
of the biorthonormal system {(ip n , n )}, we can easily show that 77+, as defined by 

N 

77+ :=£>„) (0„|, (64) 

71=1 

is a positive-definite operator belonging to 971,4. As we explained in Subsection 2.4, this is the unique metric operator 
whose inner product makes {ip n } orthonormal. 

Again if DJXa includes a metric operator 77+, then 77+ = B^ij+B for any everywhere-defined, bounded, invertible 
operator B commuting with A is also a metric operator belonging to 971,4. This shows that the subset 97t^ of VJIa 
that consists of metric operators is either empty or has an infinity of elements, [146, 148]. The same holds for 
il^ := 971^ n In summary, for a Hilbert space with dimension N > 1, either there is no metric operator 77+ 
satisfying (63) or there is an infinite set of such metric operators that in addition fulfil max[cr(?7 + )] = 1. 

One may generalize the notion of the inner product by replacing Condition (j) of Subsection 2.1 by the following 
weaker condition. 

34 For a given 771 £ OTyt there always exists ?** £ R such that the spectrum of 7] rt is bounded above by 1 This follows from the fact that 
because r/i is an invertible bounded self-adjoint operator, its spectrum cr(r]i) is a compact subset of R excluding zero, [127]. Let ai and 
0:2 be respectively the minimum and maximum values of o-(r/i). If 02 > 0, we take r* := ; if 02 < 0, we take := —a^ 1 . 

35 This is true unless the Hilbert space is onc-dimcnsional. 

36 The term "77-pseudo-Hermitian operator" is used to emphasize that one works with a particular pseudo-metric operator. It coincides 
with the notion of a "J-Hermitian operator" used by mathematicians [192, 18] and the old notion of a "pseudo-Hermitian operator" used 
occasionally in the context of indefinite-metric quantum theories [65, 223]. 
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(i) (-|-) is nondegenerate, i.e., given ip £ % the condition "(</>|V>) = for all <j) € W implies :: ip = 0". 

A function -< -\- >~: % x H — > C which satisfies conditions (i), (ii) and (Hi) is called a pseudo-inner product. Clearly 
every inner product on % is a pseudo-inner product. The converse is not true, because in general there are pseudo- 
inner products -<-|->- that fail to satisfy (i). This means that there may exist nonzero ip € H such that -<^\i^>< 0. 
Such a pseudo-inner product is called an indefinite inner product. It is not difficult to see that given a pseudo-metric 
operator 77 S 971/, the following relation defines a pseudo-inner product on %. 

M-y=(-W) =:(■[) n . (65) 

Let A : H — > H be a densely-defined operator, rj : H — > "H be a pseudo-metric operator, and ^ := 1 ^4^^ Then 
for all i^i G ^(^4) and ^2 S ^^r)); we have (il>\\A\.ij)2) = (^1 IV^),, ■ In particular if A is 77-pseudo-Hermitian, 
A\. = A and 

= (A^ 2 ) v . (66) 

This means that every pseudo-Hcrmitian operator A is Hermitian with respect to the pseudo-inner product 
defined by an arbitrary element rj of 9JIa- It is not difficult to see that the converse is also true: If r/A and A^rj have 
the same domains and A satisfies (66) for some 77 G 9Jtj, then A is pseudo-Hcrmitian and r\ e 9JIa- 

An indefinite-metric space is a complex vector space V endowed with a function -< -| - >~: % x % — > C satisfying 
(i), (ii) and (Hi), [53, 18]. One can turn a Hilbert space into an indefinite-metric space by endowing the underlying 
vector space with an indefinite inner product whose pseudo-metric operator rj is not positive-definite. The latter 
is called an indefinite metric operator. It is important to realize that the study of general indefinite-metric spaces is 
not the same as the study of consequences of endowing a given Hilbert space with an indefinite inner product. The 
latter, which is known as the ^-formalism, avoids a host of subtle questions such as the nature of the topology of 
indefinite-metric spaces [183, 185]. 

The indefinite-metric quantum theories [190, 223, 183, 185] involve the study of particular indefinite-metric spaces 
having a fixed indefinite inner product. In this sense they share the philosophy adopted in conventional QM; the 
indefinite-inner product is fixed from the outset and the theory is built upon this choice. The situation is just the 
opposite in pseudo-Hcrmitian QM where the space of state- vectors is supposed to have the structure of a (separable) 
Hilbert space with an inner product which is neither indefinite nor fixed a priori. 

In pseudo-Hermitian QM, the physical Hilbert space is constructed using the following prescription. First one 
endows the vector space of state- vectors with a fixed auxiliary (positive-definite) inner product. This defines the 
reference Hilbert space H in which all the relevant operators act. Next, one chooses a Hamiltonian operator that 
acts in Ji, is diagonalizable, has a real spectrum, but needs not be Hermitian. Finally, one determines the (positive- 
definite) inner products on % that render the Hamiltonian Hermitian. Because there is an infinity of such inner 
products, one obtains an infinite class of kinematically different but dynamically equivalent quantum systems. The 
connection to indefinite-metric theories is that for the specific 'PT-symmctric models whose study motivated the 
formulation of pseudo-Hcrmitian QM, there is a simple and universal choice for an indefinite inner product, namely 
the 'PT-inner product (3), which makes the Hamiltonian Hermitian. But this indefinite inner product does not define 
the physical Hilbert space of the theory. 

Clearly the basic ingredient in both the indefinite-metric and pseudo-Hcrmitian QM is the pseudo-metric operator. 
In general the spectrum of a pseudo-metric operator need not be discrete. However, for simplicity, we first consider a 
pseudo-metric operator 77 £ 9Jt/ that has a discrete spectrum. We can express it using its spectral representation as 

N 

r] = ^2e n \e n )(e n \, (67) 

ra=l 

Because rj is a bounded invertible Hermitian operator, its eigenvalues e„ are nonzero and its eigenvectors |e„) form 
a complete orthonormal basis of the Hilbert space T~L. Therefore, we can define 

N 

B:=Y J \^r 1,2 \e n )(e n \^\r 1 r 1 l 2 , (68) 

71=1 

and use it to obtain a new pseudo-metric operator, namely 

N 

77 := B^rjB = £sgn(e n ) \e n )(e n \. (69) 

n=l 

37 Failure to pay attention to this point is responsible for confusing pseudo-Hermitian QM with indefinite-metric quantum theories. See 
for example [126]. 
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The presence of a continuous part of the spectrum of 77 does not lead to any difficulty as far as the above construction 
is concerned. Because 77 is Hermitian, one can always define \r]\ := y/rf and set B := I??!" 1 / 2 . Both of these operators 
are bounded, positive-definite, and invertible. Hence rj := B^r/B is an element of 371a whose spectrum is a subset of 
{-!,!}• 

If we perform the transformation 77 — > 77 on a (positive-definite) metric operator ry + , wc find fj = I and (-|-)- = (-|-). 
This observation is used by Pauli to argue that we would only gain "something essentially new if we take into 
consideration indefinite bilinear forms • • • ," [190]. 38 To provide a precise justification for this assertion, let 77+ be 
a (positive-definite) metric operator and H v+ be the Hilbcrt space having the inner product Then, for all 

(B^Bfh)^ = (Byji\v+B^ 2 ) = (^B^B^) = l^i\f, + ^) = (^2), (70) 

where we have used the identities B := \r/ + \~ 1 / 2 = 77 + -1 / 2 and 774. := Bh]+B = I. Eq. (70) shows that B is a 
unitary operator mapping H onto H v+ . As a result, the quantum system s v whose state- vectors belong to H v+ is 
unitary-equivalent to the quantum system s r whose Hilbert space is the reference Hilbert space H. They describe 
the same physical system. This is the conclusion reached by Pauli in 1943, [190]. There is a simple objection to this 
argument. It ignores the dynamical aspects of the theory. As we will see below, the description of the Hamiltonian 
and the time-evolution operator can be very complicated in the " Hermitian representation" of the physical system. 
Therefore, although considering H v+ with a (positive-definite) metric operator i] + generally yields an equivalent 
"pseudo-Hermitian representation" of the conventional QM, a clever choice of 77+ may be of practical significance 
in deriving the physical properties of the system under investigation. As we discuss in Subsection 9.2, it turns out 
that indeed these new representations play a key role in the resolution of one the oldest problems of modern physics, 
namely the problem of negative probabilities in relativistic QM of Klein-Gordon fields [150, 151, 161, 179] and certain 
quantum field theories [35] . 

We end this subsection with the following remarks. 

• Strictly speaking Pauli's above-mentioned argument does not hold, if one keeps 77, to be positive-definite but 
does not require it to be invertible or bounded [123, 211]. For example one might consider the case that 
exists but is unbounded. In this case, r/ + is not onto and B fails to be a unitary operator. This type of 
generalized metric operators and the corresponding non-unitary transformations B have found applications in 
the description of resonances [133, 12, 196]. They also appear in the application of pscudo- Hermitian quan- 
tum mechanics for typical 'PT-symmetric and non-'PT-symmetric models. For these models the Hamiltonian 
operator is a second order differential operator H acting in an appropriate function space T that renders the 
eigenvalue problem for H well-posed. As discussed in great detail in [160], if H is to serve as the Hamiltonian 
operator for a unitary quantum system, one must construct an appropriate reference Hilbert space % in which 
H acts as a quasi-Hermitian operator. This in particular implies the existence of an associated metric operator 
that satisfies the boundedness and other defining conditions of the metric operators. 

Suppose H' is a differential operator acting in a function space T and having a discrete spectrum, i.e., there is 
a countable set of linearly-independent eigenfunctions of H' with isolated non-degenerate or finitely degenerate 
eigenvalues. We can use H' and T to define a unitary quantum system as follows [149, 175]. 

First, we introduce & to be the subset of T that contains the eigenfunctions of H' with real eigenvalues, and 
let C be the span of j^, i.e., C :— j J2m=i Cmipm | -A/ € %> + , c m <= C, ip m £ & |. Next, we endow C with the 
inner product [129] 

J K min(J,K) 

' c ^ I dk ^' k ) := c* m d m , (71) 



k=l 



and Cauchy-complete 39 the resulting inner product space into a Hilbert space JC. We can then identify the 
restriction of H' onto C, that we denote by H, with the Hamiltonian operator of a quantum system. It is a 
densely-defined operator acting in /C, because its domain has a subset C that is dense in IC, [195]. In fact, in 
view of (71) and the fact that C is dense in IC, the eigenvectors of H form an orthonormal basis of IC. Moreover, 
H : IC — > K, has, by construction, a real spectrum. Therefore, it is a Hermitian operator. 

This construction is quite difficult to implement in practice. Instead, one takes the reference Hilbert space % 
to be an L 2 -space, ensures that the given differential operator that is now denoted by H has a real spectrum, 
and that the set of its eigenfunctions & is dense in 7i. Then, one constructs an invertible positive operator 77+ 



38 Pauli uses the "term bilinear form" for what we call a "pseudo-inner product". 

39 Every separable inner product space N can be extended to a separable Hilbert space K., called its Cauchy completion, in such a way 
that N is dense in /C and there is no proper Hilbcrt subspace of IC with the same properties, [195]. 
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satisfying the pseudo-Hermiticity condition, 



H^= V+ H V -\ (72) 

and uses it to construct the physical Hilbert space and the Hermitian representation of the system. 

For most of the concrete models that have so far been studied the obtained r/ + turns out not to be everywhere- 
defined or bounded. But, these qualities are highly sensitive to the choice of the reference Hilbert space that 
may also be considered as a degree of freedom of the formalism. The mathematical data that have physical 
content are the eigenfunctions of H and their linear combinations, i.e., elements of C Therefore, the only 
physical restriction on the reference Hilbert spaces TL is that £ be a dense subset of TL. This means that the 
question of the existence of a genuine metric operator associated with H requires addressing the problem of 
the existence of a (reference) Hilbert space TL such that 

1. £ is a dense subset of TL, and 

2. there is metric operator ij + : TL — >• TL satisfying (72). 

It is also possible that given a Hilbert space TL fulfilling the first of these conditions and an unbounded invertible 
positive operator rj + : TL — > TL satisfying (72), one can construct a genuine bounded metric operator fulfilling the 
latter condition. These mathematical problems require a systematic study of their own. Following physicists' 
tradition, we shall ignore mathematical subtleties related to these problems when we deal with specific models 
that allow for an explicit investigation. 

Let A be a densely-defined linear operator with a non-empty il^ and r) + G Because B^ 1 := rj+ 2 : TL V — > % 
is a unitary operator and A : TL — > T~L, 1+ is Hermitian, the operator a:~ B~~ l AB is a Hermitian operator 
acting in TL. This shows that A : TL — >• TL is related to a Hermitian operator a: TL — >• TL via a similarity 
transformation, 

A = BaB-\ (73) 
Such an operator is called quasi- Hermitian, [207]. 40 

Let B := i] + X ^ 2 , B' : TL — > TL n be an arbitrary unitary operator, and a 1 := B' 1 AB' . Then in view of (73), 
a' = U~ 1 aU, where U : TL — > TL is defined by U := B~ l B' . Because both B and B' are unitary operators 
mapping TL onto TL , 0! and U are respectively Hermitian and unitary operators acting in %. Conversely, for 
every unitary operator U : % — > % the operator B' := BU is a unitary operator mapping TL onto and 

a' := B' AB' is a Hermitian operator acting in Ji. These observations show that the most general Hermitian 
operator a' : TL — > % that is related to A via a similarity transformation, 

A = B'a'B'~ 1 , (74) 

has the form 

B' = BU = i]+ 1/2 U, (75) 

where U is an arbitrary unitary transformation acting in TL, i.e., U G U(TL). If we identify A with the 
Hamiltonian operator for a quantum system and employ the formalism of pseudo-Hcrmitian QM, the metric 
operator n + defines the physical Hilbert space as "Hphys '■— %, J+ • Being Hermitian operators acting in TL p hy S , the 
observables O can be constructed using the unitary operator B : TL — > TL n and Hermitian operators o : TL — > TL 
according to 

O = B0B- 1 = ? 7 + 1/2 or 7 + /2 . (76) 

One can use any other unitary operator B' : TL — > TL n for this purpose. Different choices for B' correspond to 
different one-to-one mappings of the observables O to Hermitian operators o. According to (75) if o = B~ x O B, 
then d := B' l O B' = U~ 1 oU. Therefore making different choices for B' corresponds to performing quantum 
canonical transformations in TL. This in turn means that, without loss of generality, we can identify the physical 
observables of the system in its pseudo-Hcrmitian representation using (76). 



40 A linear densely-defined operator A : "H — > H acting in a Hilbert space T-L is said to be quasi- Hermitian if there exists an everywhere- 
defined, bounded, invertible linear operator B •.TL—^'H and a Hermitian operator Q-.'H—t'H such that A = BCLB" 1 . The above analysis 
shows that A is quasi-Hermitian if and only if it is pseudo-Hermitian and is nonempty. In mathematical literature, the term quasi- 
Hermitian is used for bounded operators A satisfying A^T = TA for a positive but possibly non-invcrtible linear operator T, [76]. These 
and their various generalizations and special cases have been studied in the context of symmetrizable operators [197, 131, 213, 239, 125]. 
For a more recent review see [109]. 
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3.3 Spectral Properties of Pseudo-Hermitian Operators 

Consider a pseudo-Hermitian operator A acting in an TV-dimensional separable Hilbert space H, with N < oo, and 
let n £ 9JLa. The spectrum of A is the set <j(A) of complex numbers A such that the operator A — XI is not invertiblc. 
Let A £ cr(A), then A — XI is not invertible and because rj is invertible, r)(A — A/)n _1 = A^ — XI must not be 
invertible. This shows that A £ a(A'). But the spectrum of A^ is the complex-conjugate of the spectrum of A, 
i.e., A £ tr(A^) if and only if A* £ tr(A), [123]. This argument shows that as a subset of the complex plane C, the 
spectrum of a pseudo-Hermitian operator is symmetric under the reflection about the real axis, [18, 53]. In particular, 
the eigenvalues a n of A (for which A — a n I is not one-to-one) are either real or come in complex-conjugate pairs, 
[190, 143]. 

Let A be a diagonalizablc pseudo-Hermitian operator with a discrete spectrum [195]. 41 Then, one can use a Ricsz 
basis {ipn} consisting of a set of eigenvectors of A and the associated biorthonomal basis {<p n } to yield the following 
spectral representation of A and a pseudo-metric operator ?] G 9JLa, [143, 146]. 

Wo Af 

A = E anolV'noX^nol +^{ a v\ll>v)(<l>v\ + {<t>-v\) , (77) 

no— 1 f—1 
N Af 
T) := E CT «oKo)(0nol+ + \4>-v)(<t>»\), (78) 

no — 1 v=\ 

where no labels the real eigenvalues a no (if any), v labels the complex eigenvalues a v with positive imaginary part 
(if any), — v labels the complex eigenvalues ol- v = a* with negative imaginary part, the eigenvalues with different 
spectral labels n £ {no, v, — v} need not be distinct, < Nq,M < 00, er„ £ { — 1,1} are arbitrary, and we have 

A\ip no ) = a no \^no), A\i/}± v ) = a± v \ip± v ), (79) 
(k |VO = <W«o for au m o, no £ {1, 2, 3, • • • , N }, (80) 
(<?WIVv) = S 5,t) <W for all g, []£{-,+} and u, ^ £ {1,2,3, ••• ,AT}. (81) 

Consider the set £({tp n }) of finite linear combinations of V'n's, as defined by (14). According to (77), elements of 
£({ip n }) belong to the domain of A, i.e., £({ip n }) C V{A). But because {ip n } is a basis, £({?/>«}) is a dense subset 
of H. This implies that V{A) is a dense subset of H. 

Next, we show that the operator rj defined by (78) does actually define a pseudo-metric operator, i.e., it is an 
everywhere-defined, bounded, invertible, Hermitian operator. 

Let {xn} be an orthonormal basis of W and B : T-L — > H be the everywhere-defined, bounded, invertible operator 
that maps {xn} onto the Riesz basis {ip n }, i.e., ip n = Bx-n for all n. It is not difficult to see that the biorthonormal 
basis {4>n} may be mapped onto {x-n} by B', Xn — B'(j) n for all n. We can use this relation and (78) to compute 

No Af 

fj := B^B = o no |xno)(Xnol + E (M(X-,\ + \X-v)(xA) ■ (82) 

no— 1 f— 1 

It is not difficult to show that fj is a pseudo-metric operator, fj £ 9Jtj. To see this, let tp € % be arbitrary. Then 
because {xn} is orthonormal, 

N No Af 

$ = ^{Xn\i>) Xn = YliXno^Xno +^{(xAi>) Xv + (X-v\i>) X-v) ■ (83) 
n— 1 no — 1 u=\ 

In view of (82) and (83), we have 

N Af 

f)lp = E O'no(XnolV') Xn + E i(X-u\lp) X v + (Xu\i>) X-u) ■ (84) 
no— 1 u—1 

In particular, 

No Af 

ii w n 2 = E Kxnd^i 2 + E (\(x-»m 2 + \(x,w\ 2 ) = u ii 2 • (85) 

no — 1 u=\ 

This shows that not only fj is everywhere-defined but it is bounded. Indeed, we have || fj ||= 1. Because fj is a 
bounded everywhere-defined operator, ff is also everywhere-defined and as is obvious from (82), it coincides with fj, 



'This in particular implies that the eigenvalues of A have finite multiplicities. 
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i.e., f] is Hermitian. Finally, in view of (82), we can easily show that fj 1 = I. In particular, rj~ l = fj is bounded and 
fj is invertible. This completes the proof of 77 £ DJli. 

Next, we observe that 77 £ SOT/ implies 77 £ DJlj. This is because according to (82), 77 = i)B~ l and B^ 1 
and £t" are bounded everywhere-defined invertible operators. Therefore, 77 as defined by (78) is a pseudo-metric 
operator. Let us also note that the inverse of r\ is given by 

N Af 

^ X =B^ = ^ 0noWv»oXVv.ol+X! (I^X^M + |V-,)(^D- (86) 

rig — 1 ^— 1 

We can easily show that 77 belongs to 931a by substituting (77), (78), and (86) in r\Ar\~ x and checking that the result 
coincides with A> . 

In Ref. [146], it is shown that any element 77 £ 371,4 can be expressed in the form (78) where {tpn} is the 
biorthonormal basis associated with some (Riesz) basis {ip n } consisting of the eigenvectors of A. As any two Ricsz 
bases are related by a bounded everywhere-defined invertible operator L : % — > % one may conclude that the elements 
of DJIa have the following general form. 

r/ = L^L, (87) 

where 77 is the pseudo-metric operator (78) that is defined in terms of a fixed (but arbitrary) (Riesz) basis {tpn} 
consisting of the eigenvectors of A and set of arbitrary signs. The operator L appearing in (87) maps 

eigenvectors tjj n of A to eigenvectors Lip n of A in such a way that ijj n and Lip n have the same eigenvalue, [146]. This 
in particular implies that L commutes with A. 

Suppose A has a complex-conjugate pair of nonreal eigenvalues a±„, let ip± u be a corresponding pair of eigen- 
vectors, 77' £ 9JIa be an arbitrary pseudo-metric operator associated with A, L be an everywhere-defined, bounded, 
invertible operator commuting with A and satisfying (87), and £ := L~ lr ip v . Then (£|r/£) = (ipv\Tj tpv) = 0. Because 
L is invertible £ 7^ 0, this is an indication that 77' is not a positive-definite operator. Similarly suppose that one of 
the signs cr„ appearing in (78) is negative and let ( := L" 1 i/'n . Then (CWO = (V'nol^^no) = a n a = — 1, and again 
t]' fails to be positive-definite. These observations show that in order for 77' to be a positive-definite operator the 
spectrum of A must be real and all the signs a no appearing in (78) must be positive. In this case, we have 77' = rj' + , 
where 

N 

V+-=Lh 1+ L, r ]+ = Y,\<t>n)(4>n\- (88) 

n=l 

The choice of the signs a no is not dictated by the operator A itself. Therefore it is the reality of the spectrum of 
A that ensures the existence of positive-definite elements of DJIa ■ By definition, such elements are metric operators 
belonging to 0Jl\. Their general form is given by (88). 



3.4 Symmetries of Pseudo-Hermitian Hamiltonians 

Consider a pseudo-Hermitian operator A : % — > Ti. and let 771 and 772 be a pair of associated pseudo- metric operators; 
rjiAr/i 1 — A^ = ^Arj^ 1 ■ Then it is a trivial exercise to show that the invertible linear operator S := 77^ 1 771 commutes 
with A, [143]. If we identify A with the Hamiltonian of a quantum system, which we shall do in what follows, S 
represents a linear symmetry of A. 

Next, consider a diagonalizable pseudo-Hermitian operator A with a discrete spectrum, ip n be eigenvectors of A, 
and {(tp n , (f> n )} be the complete biorthonormal extension of {ipn} so that A admits a spectral representation of the 
form (77): 

Wo Af 
A =Y. an„\i^ a )(KA + ^{av\^v)(^\ + at\iP-„}((l)-„\) . (89) 

no — 1 u=l 

Moreover, for every sequence a = (<x no ) of signs (<r no £ {— 1, +1}), let 

N Af 

Tic ■= ^ol0n o )(^nol + X]O^)(^l + l^)(^l)> ( 90 ) 

no — 1 u=l 

N Af 

C a := J ^ol^n )(^ol + E (l^><<M + \i>-»){<t>-A), (91) 
rin — 1 u—1 
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and 771 , 6 : H — > H be denned by 

Wo Af 

no — 1 f=l 
Wo Af 

© := £ l^o>*^nol+E(lVv)*<<M + lV>-,}*^-,|), (93) 

no — 1 1/— 1 

where for every cf> 6 H the symbol * {<j>\ denotes the following antilinear operator acting in H. 

W) * (MC ■= (t\<t>) 1> = <0IC>* V>, forallCGH. (94) 

As we discussed in Section 3.3, r\ a and 771 are pseudo-metric operators associated with A. The operators C a and 6 
have the following remarkable properties [148]. 

• In view of the fact that rf[ 1 = Y.nl=i IV'noXVvJ + Y^Li (IVv) (ip- v \ + \i>-„){il)„\), we have 

Ca = C V- (95) 
Therefore, is a linear invertible operator that generates a symmetry of A: 

[C (r ,A]=0. (96) 

• Using (94) and the biorthonormality and completeness properties of {(i/i,,,^)}, we can check that 6 is an 
invertible antilinear operator that also commutes with A, 

[6,A]=0. (97) 

• C a and <3 are commuting involutions, i.e., 

[C a ,©]=0, Cl = G 2 =I. (98) 

In summary, we have constructed an involutive antilinear symmetry generator & and a class of involutivc linear 
symmetry generators C CT that commute with @. 

It turns out that if a given diagonalizable operator A with a discrete spectrum commutes with an invertible 
antilinear operator, then A is necessarily pseudo-Hermitian. Therefore, for such operators pseudo-Hermiticity and 
the presence of (involutive) antilinear symmetries are equivalent conditions [145, 219]. Furthermore, each of these 
conditions is also equivalent to the pseudo-reality of the spectrum of A. The latter means that the complex-conjugate 
of every eigenvalue of A is an eigenvalue with the same multiplicity [145, 219]. These observations are the key for 
understanding the role of VT symmetry in the context of our study. They admit extensions for a certain class of 
non-diagonalizablc operators with discrete spectrum [147, 208, 218, 62] and some operators with continuous spectrum 
[150, 159, 165]. 

The reality of the spectrum of A is the necessary and sufficient condition for the existence of an associated metric 
operator and the corresponding positive-definite inner product that renders A Hermitian [144]. For the case that the 
spectrum of A is real, the expressions for •q a ,rj\,C ai and © simplify: 

N N 



■= ^2^n\<t>n){^n\, Vl = ^2 \<t>n) (<Pn\ = V+ > (99) 

n— 1 n— 1 

N N 

= £>nW>„XM © = X>»> (100) 



n=l 



and wc find 



Co-lpn = Vnlpn, &4>n = tpn- (101) 



Hence, C a and 6 generate exact symmetries of A. 

In order to make the meaning of C a more transparent, we use the basis expansion of an arbitrary ip G H, namely 
V> = J2n=i c ni'n, to compute 

Af N 
n— 1 n— 1 
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where ip± := J2 n em± c >iV'n and 91± := jn G {1, 2, 3, • • • , N} | er„ = According to (102), for every ip G % there 

are unique state- vectors tp± belonging to the eigenspaccs H± := {tp € % \ C a ip = ±ip } of C a such that V = tp+ — ip- ■ 
This identifies C a with a (Z2-) grading operator for the Hilbert space. 

If a n = ±1 for all n, we find C a = ±1. In the following we consider the nontrivial cases: C ^ ±1. Then T~L± are 
proper subspaces of H satisfying 

n = n+®U-, (103) 

C a is a genuine grading operator, both r\ a and —r\ a fail to be positive-definite, and (•| , )ij <r is indefinite. Furthermore, 
in view of (102), the operators II± '■= \ (I ± Co-), satisfy H±i/j = ip±, i.e., they are projection operators onto H±. 

Next, consider computing (tpg\(f>b)ri a f° r arbitrary g, f) G { — ,+} and ip±,<f>± G "H±. Using the basis expansion of 
ip± and 0±, i.e., 0± = J2nevi ± d ^n and ^± = J2nem ± c n^Pn, and (99), we have r] a (f>± = ±J2nzvt ± and 

(V'ol^f))^ = (V^k<r0f)) =0<5 3 ,f) XI c « d «- ( 104 ) 

Therefore, with respect to the indefinite inner product (-\-) ntT , the subspaces H+ and H- are orthogonal, and (103) 
is an orthogonal direct sum decomposition. 

Another straightforward implication of (104) is that, for all ijj, 4> G H, 

{^\4>)n. = {i , + \4>+)r, a + {ip-\(j)-)r, a = E c* n d n - E c* n d n = (tp + \<t> + ) n+ - (i/>_\</>_) v+ , (105) 

n£<n + n£Vl- 

where ip± := H±ip, 4>± '■= n±0, c n := (</>„|?/>) and d n := (<f> n \<fi). Similarly, we have 

N 

mg^>„«, = (tm&v - w-tt-K = E < d « + E < d " = E < d » = w>^- ( 106 ) 

nem + ne<rt_ n=i 

This calculation shows that the positive-definite inner product (•| , )^ + can be expressed in terms of the indefinite 
inner product (•| , ) r ; ( , and the grading operator C a according to 



■\C«)v„- (107) 



Conversely, one can use (107) to define a positive-definite inner product that makes H Hcrmitian. The latter scheme 
may be traced back to a similar construction developed in the 1950's in the context of indefinite-metric quantum 
theories [187]. See [183, 185] for reviews. In the context of "PT-symmetric quantum mechanics, it was proposed (with 
a specific choice for the sequence a) in [37] and coined the name CPT-inncr product. 

To see the connection with the treatment of [37], consider the case that A, which is now viewed as the Hamiltonian 
operator for a quantum system, is symmetric, i.e., A T := T A^T — A, where T is the time-reversal operator. 42 Then 
given the spectral representation of A and A\ we can easily choose a biorthonormal system {(ip n , 4>n)} such that 

<t>n = VnTlpn- (108) 

Using this relation together with the biorthonormality and completeness properties of {(ipn, </>n)}, we can obtain the 
following spectral representation of T . 

N 

T= E CT «I^"MM (109) 
n=l 

In view Eqs. (99), (100), (109), and T 2 = I, we have 

n) — <3*m &n \ fin \<M, (no) 
G=TVa- (Ill) 

Clearly, we could use (109) to define an invertible antilinear operator satisfying (108) for an arbitrary possibly 
non-symmetric A, namely 

N 

T a := E a n\4>n) * (<Pn\- (112) 
n=l 

But in this more general case, (110) may not hold and T a may not be an involution. 



42 Recall that T is an antilinear Hermitian (and unitary) involution [233, 139]. 
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In fact, condition (110) is not only a necessary condition for T£ = / but it is also sufficient [148]. An analogous 
necessary and sufficient condition for rfe = I is [148] 



If both (110) and (113) hold, 
and as a result 

By virtue of this relation and 6 2 = /, 
Therefore, 



\<t>n), (H3) 

r ff 2 = ill = i, (114) 

© = (115) 

[n*,%] = o. (ii6) 

© = VaTa- (117) 

For the symmetric and "PT-symmetric Hamiltonians (1) that arc considered in [37], we can find a biorthonormal 
system {(ip n ,4>n)} satisfying (110) - (113), [148]. Moreover, setting er„ := (— 1)" +1 for all n G Z + , we have 

r]a=V- (118) 

Therefore, in light of (111) and (116), the antilinear symmetry generator © coincides with VT, 

© = VT, (119) 

and the linear symmetry generator C a is the "charge-conjugation" operator C of [37]. In view of (95), (97), (99), 
(98), and (118), it satisfies 

C 2 = I, [C,A] = [C,PT] =0, C = ? ? + 1 7>. (120) 

Furthermore, because in the position representation of the state vectors (Tip)(x) = tp(x)*, the positive-definite inner 
product (107) coincides with the CPT-hmer product. 

This completes the demonstration that the CVT- inner product is an example of the positive-definite inner products 
that we explored earlier. 

We conclude this subsection by noting that although in general we can introduce a pair generalized time-reversal 
and parity operators, r\ a and To, they may fail to be involutions. 



3.5 A Two-Level Toy Model 

In this subsection we demonstrate the application of our general results in the study of a simple two- level model which, 
as we will see in Subsection 9.2, admits physically important infinite-dimensional generalizations [150, 161, 151]. 

Let T-L be the (reference) Hilbert space obtained by endowing C 2 with the Euclidean inner product and {ei,e 2 } 
be the standard basis of C 2 , i.e., e\ :=( * V e 2 :=( " V Then we can represent every linear operator K acting in 

T-L in the basis {ei,e2} by a 2 x 2 matrix which we denote by Kj, the entries of K_ have the form K_ irj := (ei\Kej) 
where i, j G {1, 2}. 

Now, consider a linear operator A : C 2 — > C 2 represented by 

A: =K-dYi -d-\)' 

where D is a real constant. A : H — > T-L is a Hermitian operator if and only if D = 1. We can easily solve the 
eigenvalue problem for A. Its eigenvalues a n and eigenvectors tp n have the form 

ai = -a 2 = D>, Vi = ci I 1 __ D i/2 J' ^ = 62 1 + D i /2 I, (122) 

where C\,C2 are arbitrary nonzero complex numbers. Clearly, for D = 0, a\ = a 2 = 0, the eigenvectors become 
proportional, and A is not diagonalizable. For D ^{), A has two distinct eigenvalues and {^1,^2} forms a basis of 
C 2 . This shows that D = marks an exceptional spectral point [123, 103], for D > the eigenvalues are real and 
for D < they are imaginary. 

It is an easy exercise to show that A is E3-pseudo-Hermitian where S3 : C 2 — > C 2 is the linear operator represented 
in the standard basis by the Pauli matrix (T3 :=( J -1 )" H ence includes £3, and A is £3-pseudo-Hcrmitian 
for all DeM. This is consistent with the fact that £3 is not a positive-definite operator, because otherwise A could 
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not have imaginary eigenvalues. According to our general results, for D > 0, 971,4 must include positive-definite 
operators. To construct these we first construct the biorthonormal basis {4>i,(j) 2 } associated with {ipi,ip 2 }. For 
D > 0, the basis vectors 4> n are given by 

/ 1 -i- n- 1 / 2 \ / 1 _ n- 1 / 2 \ 

0i = (4CD" 1 [ J _ B-1/2 J , 02 = (4c;)" 1 ^ J + 2- V 2 j • (123) 

Inserting these relations in (78) we find the following expression for the matrix representation of the most general 
pscudo- metric operator 77 G 9)1 a- 



B = rioi^ 1 _ D _ 1 {1 _ D . 1/2)2 )+r 2 a 2 \^ 1 _ D _, {1 + D - 1/2)2 ), (124) 

where r\ := 1 4ci | 2 and r 2 := |4c2 1~ 2 are arbitrary positive real numbers and o\,a 2 are arbitrary signs. 

The choice o\ = —o 2 = 1 and r± = r 2 = D 1 / 2 /4 yields 77 = S3. The choice o\ = a 2 = 1 yields the form of 
the most general positive-definite element of StJt^. A particularly simple example of the latter is obtained by taking 
ci = c 2 = D~ 1 / 4 /2 which implies n = r 2 = Z? 1 / 2 /4. It has the form 

1 / D l ' 2 + D- 1 ' 2 D^-D- 1 ' 2 \ 
R + ~ 2 \ D 1 ' 2 - D- 1 ' 2 D^ + D- 1 / 2 J' (125j 

We can simplify this expression by introducing 43 9 := ^ In D, and using the fact that the Pauli matrix <J\ := ( J 1 
squares to the identity matrix J. This yields 

/ cosh6* sinh^ \ . . T . , . rr, 

n+ = ^ s . nh , cogh , j =cosh^ + smh^ l=e ^ 1 . (126) 

The metric operator represented by (125) and (126) defines the following (positive-definite) inner product on C 2 
with respect to which A is a Hermitian operator. 

(z\w} v+ := (z\r)+w) = (z* W\ + z 2 w 2 ) cosh# + {z{w 2 + z 2 w{) sinh 6, (127) 

where z = (z\, z 2 ) T ,w — (wi,w 2 ) T £ C 2 are arbitrary. The inner product (127) has a more complicated form than 
both the reference (Euclidean) inner product, (z\w) = z*w\ + z^wi^ and the indefinite inner product defined by S3, 

(_?|mJ) S3 := (z\T, 3 w) = z*W\ - z* 2 w 2 . (128) 

Furthermore, unlike (-|-) and , the inner product (-|-)^ depends on 6 and consequently D. In particular, 

as D — > it degenerates. In fact, a quick inspection of Eq. (124) shows that every D-indepcndcnt pseudo-metric 
operator is proportional to £3 and hence necessarily indefinite; ±£3 are the only /^-independent elements of Ua- 

Now, suppose that A is the Hamiltonian operator of a two-level quantum system. If we employ the prescription 
provided by the indefinite-metric quantum theories, we should endow C 2 with an indefinite inner-product from the 
outset. The simplest choice that is historically adopted and viewed, according to the above-mentioned argument due 
to Pauli, as being the unique choice is 77 = S3. In this case the system has a physical state corresponding to the state- 
vector ei and a hypothetical state or ghost corresponding to e 2 . Unfortunately, the subspace of physical state- vector, 
i.e., the span of {ei}, is not invariant under the action of A unless D = 1. Hence, for B ^ 1, such an indefinite- metric 
quantum theory suffers from interpretational problems and is inconsistent. In contrast, pseudo-Hcrmitian quantum 
mechanics provides a consistent description of a unitary quantum theory based on the Hamiltonian A. This is done 
by endowing C 2 with the (positive-definite) inner product (-l - )^' , where ri' + is given by the right-hand side of (124) 
with (Ji = a 2 = 1. It involves the free parameters r\ and r 2 that can be fixed from the outset or left as degrees of 
freedom of the formulation of the theory. We can represent this most general metric operator by 

/ cosh0 + s sinhfl \ 
1J -+- 7 \ sinhe* cosh6>- Sy !' { J} 

where r := 2{r\ + r-^D^ 1 ! 2 € M + and s := e ( — -*-> •"•) are arbitrary. We can use (129) to determine the most 

general inner product on C 2 that makes A Hermitian. This is given by 

{z\w) v , := (z\r]' + w) = r (z\w) n+ + s {z\w) Sa , (130) 



s This was pointed out to me by Professor Haluk Beker. 
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where (•|-)„ + and (-|-) S3 are respectively defined by (127) and (128). 

We can relate rj' + to rj + using a linear operator L : C 2 — > C 2 commuting with A via r/' + = L^r] + L. This operator 
has the following general form L = 2Z?~ 1 / 4 [y/rie^ 1 \tpi)(^>i \ + ^/rje 1 ^ 2 , where <fi,ip2 G [0, 2tt) are arbitrary. 

We can represent it by 

A- cosh 9 + A+ A_ sinh 9 
—A_ sinh 9 —A- cosh 9 + A4 



(131) 



where A± := D" 1 / 4 (V? T 1 e^ 1 ± y/¥ 2 e lip2 ). With the help of (126), (129), and (131), we have checked that indeed 

nL+ =L} r+l. 

Next, we wish to establish the quasi- Hermiticity of A for D > 0, i.e., show that it can be expressed as A = BaB 1 
for an invertible operator B : C 2 — > C 2 and a Hermitian operator a : H — > H. As we explained in Section 3.2, we can 
identify B with the inverse of the positive square root of a metric operator belonging to 9Jtj^ . A convenient choice is 

— 1/2 

B = ?7 + , for in light of (126) we have 

B ±x = gT f a 1 = ( cosh f =F sinh f 
\ =F sinh I cosh | 

This leads to the following remarkably simple expression for the matrix representation of (X. 

a = B-'AB={^ _° 1/2 ) = D^a 3 . (132) 

Hence, a = L> 1/2 £ 3 . 

Because every Hermitian operator o : % — > % is a linear combination, with real coefficients, of the identity 
operator / and the operators £1, £2, and £3 that are respectively represented by Pauli matrices (J 1; cr 2 , and C3, we 
can express every physical observable O : —> H v as O = a^I + Ylj=i a j where Oq, cti, ct2, 03 are some real 

numbers and Sj := BY^jB^ 1 = 77, 1 £j r) 1 / 2 for all j G {1,2,3}. The observables Sj are represented by Si = u u 
S?. = cosh 9 G\ — i sinh 9 (T3 , and S^_ = i sinh 9 0~2 + cosh (T3 . 

Next, we repeat the calculation of the Hermitian Hamiltonian and the physical observables for the case that we 
choose the general metric operator rj' + to construct the physical Hilbcrt space, i.e., set "Hphys = H v , ■ The matrix 

representation of the Hermitian Hamiltonian a! \H — > % is then given by 

' - W ^ & n> -V2 _ r.1/2 ( u(9 , s) v(6,s) \ _ 1/2 



gL = r£S"Al£>r'-" = DV*( l^l> _^ s) )=D^[v(0,s)a 1 + u(e,s)a 3 ], (133) 



O' = OqI + y^djS'j, (134) 



where u(9,s) := e e + a p 3S^l and v{9,s) := sJs2 ^^^ }kl - We can similarly express the physical 

observables in the form 

3 

X]u? up 

3=1 

where S' 3 := V+~ 1/2 ^ 3 ti' + 1 ' 2 . 

As seen from (133) the Hermitian Hamiltonian a describes the interaction of a spin i particle with a magnetic 
field that is aligned along the unit vector (u, 0, v) T in K 3 . As one varies s the magnetic field rotates in the x-z plane. 
It lies on the z-axis for s = which corresponds to using ij + to define the physical Hilbert space. Clearly, there is 
no practical advantage of choosing s ^ 0. Furthermore, for all s G ( — 1, 1) and in particular for s = 0, the Hermitian 
representation of the physical system is actually less complicated than its pscudo-Hcrmitian representations. This 
seems to be a common feature of a large class of two-level systems [149]. 

Next, we compute the symmetry generator C a for o\ = — a% = 1. Denoting this operator by C for simplicity, 
realizaing that C = |?Ai)(0i| — IV^X^I, and using (122), (123), and D = e 2e , we find 

( cosh 6» sinh6» \ 
- ^ -sinb.0 -coshtf ) ■ 

Observing that in view of (121), A 2 = DI, and making use of this relation and (135) we are led to the curious relation 
[161] 

C = ^L. (136) 

Va 2 

It is important to note that in performing the above calculation we have not fixed the normalization constants 
Ci and C2 appearing in (122) and (123). Therefore, up to an unimportant sign, C is unique. 
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As we mentioned above setting c\ — c 2 = D~ 1 / 4 /2 and o\ = —02 = 1, we find the pseudo-metric operator 
Vcr = \4>i)(4>i \ — I <?^2 ) (02 1 = S3. Because S3 is a linear involution we can identify it with P. It is a simple exercise to 
show that V \C and 77+ actually satisfy 

C = rtfV. (137) 

We can similarly construct the antilinear symmetry generator 6. It turns out that unlike C, & depends on (the 
phase of) normalization constants Ci and c 2 that appear in (122) and (123). Setting C\ = c 2 = D~ 1 ' 4 /2, we find 

6 = T, (138) 

where T denotes complex conjugation, Tz = z* . In view of (138), the symmetry condition [<5,A] = corresponds 
to the statement that A is a real operator, i.e., A is a real matrix, which is a trivial observation. 

Similarly, we can introduce an antilinear operator Ta according to (112). This yields T a '■= |0i)*(</>i| — |02)*(</>2|, 
which in general depends on the choice of ci and C2. For c\ = c 2 = D~ 1 / 4 /2, we have Ta = VT. Combining this 
relation with (138) yields & = VTa- It is not difficult to see that indeed Ta is an involution. But it differs from the 
usual time-reversal operator T ■ Let us also point out that we could construct the pseudo-Hcrmitian quantum system 
defined by A without going through the computation oi C, V and T a operators. What is needed is a metric operator 
that defines the inner product of the physical Hilbert space. 

The construction of the C operator for two-level systems with a symmetric Hamiltonian has been initially un- 
dertaken in [50]. The C operator for general two-level systems and its relation to the metric operators of pseudo- 
Hcrmitian quantum mechanics are examined in [149]. See also [243]. A comprehensive treatment of the most general 
pseudo-Hermitian two-level system that avoids the use of the C operator is offered in [178]. 



4 Calculation of Metric Operator 

A pseudo-Hermitian quantum system is defined by a (quasi- Hermitian) Hamiltonian operator and an associated metric 
operator r] + . This makes the construction of rj + the central problem in pseudo-Hermitian quantum mechanics. There 
are various methods of calculating a metric operator. In this section, we examine some of the more general and useful 
of these methods. 

4.1 Spectral Method 

The spectral method, which we employed in Section 3.5, is based on the spectral representation of the metric operator: 

N 

?7 + =^|0n)(^|- (139) 
1 

It involves the construction of a complete set of eigenvectors <p n of A^ and summing the series appearing in (139) (or 
performing the integrals in case that the spectrum is continuous). 



4.1.1 'PT-symmetric infinite square well 

The first pseudo-Hermitian and "PT-symmetric model with an infinite-dimensional Hilbert space that has been treated 
within the framework of pseudo-Hcrmitian quantum mechanics is the one corresponding to the PT-symmetric square 
well potential [241, 19]: 

v(x) = ( ^ Sgn(x) f 1 ' \ X \ < L T % (140) 
I 00 for \x\ > L/2, 

where C and L are real parameters, L is positive, and x takes real values. This was achieved in [175] using the 
spectral method combined with a certain approximation scheme that allowed for a reliable approximate evaluation 
of a metric operator as well as the corresponding equivalent Hermitian Hamiltonian and pseudo-Hermitian position 
and momentum operators. A more recent treatment of this model that makes use of the spectral method and obtains 
a perturbative expansion for a C operator and the corresponding metric operator rj + in powers of £ is given in [51]. 

4.1.2 PT-symmetric barrier 

In [159], the spectral method has been used for treating a pseudo-Hcrmitian quantum system defined by the scattering 
potential: 

V[X) ~\ for \x\ > L/2, ll41j 
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where again C>L,x £ R and L is positive. This potential was originally used in [199] as a phenomenological tool 
for describing the propagation of electromagnetic waves in certain dielectric wave guides. 44 It is the first example 
of a "PT-symmetric potential with a continuous spectrum that could be studied thoroughly within the context of 
pseudo-Hcrmitian quantum mechanics. 

Application of the spectral method for this potential involves replacing the sum in (139) with an integral over 
the spectral parameter and taking into account the double degeneracy of the energy levels. The extremely lengthy 
calculation of a metric operator for this potential yields the following remarkably simple expression [159]. 

{x\T)+\y) = S(x -y) + *$g£ (2L + 2\x + y\ - \x + y + L\ - \x + y - L|) sgn(z -y) + 0(C 2 ), (142) 

where x,y GM., and 0(C 2 ) stands for terms of order C, 2 and higher. 

An unexpected feature of the scattering potential (141) is that the corresponding equivalent Hermitian Hamilto- 
nian has an effective interaction region that is three times larger than that of the potential (141). In other words, in 
the physical space, which is represented by the spectrum of the pseudo-Hcrmitian position operator, the interaction 
takes place in the interval [— ^] rather than [— |r]- 



4.1.3 Delta-function potential with a complex coupling 

Another complex scattering potential for which the spectral method could be successfully applied is the delta-function 
potential [165]: 

v(x)=iS(x), (143) 

where 3 is a complex coupling constant with a non- vanishing real part. For this system it has been possible to 
compute a metric operator and show that it is actually a bounded operator up to and including third order terms in 
the imaginary part of 3. It is given by 

(x\n+\y) = S(x - y) + ^ [e(xy) e~ K ^ + 6{-xy) e^^} sgn(y 2 - x 2 ) + 0(( 2 ), (144) 

where £ := 3(3), K := mh~ 2 ?R.(%), -ft and 3 denote the real and imaginary parts of their arguments, 9 is the step 
function defined by 0{x) := [1 + sgn(x)]/2 for all x £ K, and we have omitted the quadratic and cubic terms for 
brevity. 

In order to determine the physical meaning of the quantum system defined by the potential (143) and the metric 
operator (144), we should examine the Hermitian representation of the system. The equivalent Hermitian Hamiltonian 
is given by [165] 

h= 2^ + mS{x) + ! W h2 + °^ 3 ^ ( 145 ) 

where 

(M)(a?) := H e ~" N + Mfc). ( 146 ) 

V> £ L 2 (R) and x £ R are arbitrary, := ip(0), and := e~ K ^ip(y) dy. As seen from (145) and (146) the 
nonlocal character of the Hermitian Hamiltonian h is manifested in the "0-dependence of the coefficients and b^. 

A generalization of the delta-function potential (143) that allows for a similar analysis is the double-delta function 
potential: v(x) = $-5(x + a) + %+5(x — a), where 3± and a are complex and real parameters, respectively [177, 138]. 
Depending on the values of the coupling constants 3±, this potential may develop spectral singularities [184, 132]. 
These are the points where the eigenfunction expansion for the corresponding Hamiltonian breaks down [177]. This is 
a well-known mathematical phenomenon [184, 132] with an interesting and potentially useful physical interpretation: 
A spectral singularity is a real energy where both the reflection and transmission coefficients diverge. Therefore it 
corresponds to a peculiar type of scattering states that behave exactly like resonances: They are resonances with 
zero- width [172, 173] 45 



4.1.4 Other Models 

The application of the spectral method for systems with an infinite-dimensional Hilbcrt space is quite involved. If the 
system has a discrete spectrum it requires summing complicated series, and if the spectrum is continuous it involves 
evaluating difficult integrals. This often makes the use of certain approximation scheme necessary and leads to 
approximate expressions for the metric operator. A counterexample to this general situation is the quantum system 

44 The use of complex potential in constructing various phenomenological models and effective theories has a long history. For a 
discussion that is relevant to complex scattering potentials see the review article [181] and [2, 72, 63]. 

45 The single delta-function potential (143) develops a spectral singularity for imaginary values of 3. For other examples of complex 
potentials with a spectral singularity, see [202, 172, 173]. 
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describing a free particle confined within a closed interval on the real line and subject to a set of PT-symmetric 
Robin boundary conditions [128]. For this system the spectral method may be employed to yield a closed formula 
for a metric operator. Other systems for which the spectral method could be employed to give an explicit and 
exact expression for the metric operator are the infinite-dimensional extensions of the two-level system considered 
in Section 3.5 where D is identified with a positive-definite operator acting in an infinite-dimensional Hilbcrt space 
[150, 151]. These quantum systems appear in a certain two-component representation of the Klein-Gordon [140] and 
(minisupcrspace) Wheeler-DeWitt fields [141]. 



4.2 Perturbation Theory 

The standard perturbation theory has been employed in the determination of the spectrum of various complex 
potentials since long ago [60] . 46 In the present discussion we use the term "perturbation theory" to mean a particular 
perturbative method of constructing a metric operator for a given quasi-Hermitian Hamiltonian operator. This 
method involves the following steps. 

1. Decompose the Hamiltonian H in the form 

H = H + eH 1 , (147) 

where e is a real (perturbation) parameter, and Hq and H\ are respectively Hcrmitian and anti-Hermitian 
e-independent operators. 

2. Use the fact that 77+ (being a positive-definite operator) has a unique Hcrmitian logarithm to introduce the 
Hcrmitian operator Q := — ln^ + , so that 

V+ = e- Q , (148) 
and express the pseudo-Hermiticity relation = r/ + Hri^ 1 in the form 

H ] =e- Q He Q . (149) 

In view of the Backer-Campbell- Hausdorff identity [200], 

e-®H e« = H + j } [H, Q] t = H + [H, Q] + ±[[H, Q],Q} + ±[[[H, Q],Q],Q] H , (150) 

1=1 

where [H, Q]i := [[■ ■ • [[H, Qj, Q], ■ ■ ■ ], Q] and I is the number of copies of Q appearing on the right-hand side 
of this relation, (149) yields 



H^ = H + f^j l [H,Q] t . (151) 
3. Expand Q in a power series in e of the form 



ev 



where Qj are e-independent Hcrmitian operators. 

4. Insert (147) and (152) in (151) and equate terms of the same order in powers of e that appear on both sides of 
this equation. This leads to a set of operator equations for Qj which have the form [162] 

llu-(), /.',. (153) 
Here j G Z + and Rj is determined in terms of Hi and Qfc with k < j according to 

J -2Hi for j = l, _ A ^ {-l) n n k m\ 

3 ^ I Ei-2 QkZ ki for j > 2, Qk 2-< 2-, k]2 m- i n!(m _ n) \ > ^ 



Y, [[[■■■[H ,Qsi},Q S2 ],--- ,},Qs h }. (155) 



Si-i hs fc =J 



3 For more recent developments, see [45, 61, 59] and references therein. 
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More explicitly we have 

[H ,Qi] = -2Hi, (156) 
[H ,Q 2 ] = 0, (157) 

[Ho,Q 3 ] = -~[#i,Qi] 2 , (158) 
[H ,Q4 = -~([[H 1 ,Qi],Q 2 ] + [[H 1 ,Q 2 ],Q 1 ]), (159) 

D 

[H ,Qs] = ^[Hi,Qil-^([H 1 ,Q 2 ] 2 + [[H 1 ,Q 1 ],Q 3 ) + [[H 1 ,Q 3 ],Q 1 ]). (160) 

5. Solve the above equations for Qj iteratively by making an appropriate ansatz for their general form. 

A variation of this method was originally developed in [40] to compute the C operator for the following VT- 
symmctric Hamiltonians and some of their multidimensional and ficld-thcoretical generalizations. 

H = -^-p 2 + \^ 2 x 2 + iex 3 , (161) 
2m 2 

H = ^-p 2 +iex 3 , (162) 
2m 

where fi and e are nonzero real coupling constants. 
4.2.1 T'T-symmetric cubic anharmonic oscillator 

A perturbative calculation of a metric operator and the corresponding equivalent Hermitian Hamiltonian and pseudo- 
Hcrmitian position and momentum operators for the Hamiltonian (161) has been carried out in [158, 118]. 

Following [40] one can satisfy the operator equations for Qj by taking Q^i = for all i G Z + and adopting the 
ansatz 

t+i 

Q 2i+ i - {^ J ^ 2fe+1 }, (163) 

j,k=0 

where {•,•} stands for the anticommutator and Cijk are real constants. Inserting (163) in (153), one can determine 
cijk for small values of i [158]. See also [40, 118]. 

Again, to determine the physical content of the system defined by the Hamiltonian (161) and the metric operator 
rj + = e~Q, we need to inspect the associated Hermitian Hamiltonian operator [158]: 

2m 2 mfi \ 2 J \m° 16m z 

81/i 2 2 2 2 33/i 4 r 4 2l 69^i 4 2 2 2 ^ /J,® g\ 4 „, 5, (ica\ 

~^~2P x P TT~ {x ,P } - -5 — x p x —x ° e +C(e D ), (164) 

8m/ 16m 8m 4 / 

and the underlying classical Hamiltonian (60): 



Pi 1 . 2„2 , 3/2 



tt -1-0,-22,^/^ 22,2412 

2m 2 2/^ 4 \ m 



2 ( pi 18 M 2 2 4 51 M 4 42 



- <pi --r^ U 4 + (165) 



// ; -' V '»" c " c 4m " c " c 4 
If we only consider the terms of order e 2 and lower, we can express (165) in the form 

where M(x c ) := m(l + 3/i~ 4 e 2 a; 2 ) -1 = m(l — 3^~ 4 e 2 a; 2 ) + 0(e 4 ). This shows that for small values of e, the 
"PT-symmctric Hamiltonian (161) describes a position-dependent-mass quartic anharmonic oscillator [158]. This 
observation has motivated the use of non-Hermitian constant-mass standard Hamiltonians, H = p 2 /(2m) + v(x), in 
the perturbative description of a class of position-dependent-mass standard Hamiltonians [21]. 

As seen from (164), the 4-th order (in e) contribution to the equivalent Hermitian Hamiltonian h involves p 6 . It is 
not difficult to show that h = z~27Lo £ 2t , where hi is a polynomial in p whose degree is an increasing function of I. 
Therefore, the perturbative expansion of h includes arbitrarily large powers of p. This confirms the expectation that 
h is a nonlocal operator. The same holds for the pscudo- Hermitian position and momentum operators [158, 118]. 
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4.2.2 Imaginary cubic potential 



Rcf. [162] gives a perturbativc treatment of the Hamiltonian 

H=^—+iex 3 , (167) 
2m 

in which the operator equations (153) are turned into certain differential equations and solved iteratively. This method 
relies on the observation that for this Hamiltonian, Hq = p 2 /2m. Therefore, (x\[Hq , Q j]\y) = 

~lm (Px ~ ®y) ( x \Qj\v)- ^ n v i ew °f this identity and (153), we find 

{-d 2 x + d 2 y ) (xlQM = {x\R 3 \y). (168) 

Because Rj is given in terms of Hi and Qi with i < j, one can solve (168) iteratively for (x\Qj\y). Note also that 
this equation is a non-homogeneous (1 + l)-dimcnsional wave equation which is exactly solvable. 

This approach has two important advantages over the earlier perturbativc calculation of the metric operator 
for the imaginary cubic potential [40]. Firstly, it involves solving a well-known differential equation rather than 
dealing with difficult operator equations. Secondly, it is not restricted by the choice of an ansatz, i.e., it yields the 
most general expression for the metric operator. In particular, it reveals large classes of CVT and non-CPT-inner 
products that were missed in an earlier calculation given in [40] . Here we give the form of the equivalent Hermitian 
Hamiltonian associated with the most general admissible metric operator: 



tz + ^ ( {A \} + 22ft 2 {^ 4 , ~ 4 } + a 2 h\x 2 , 1}- 



^,-} + 22h 2 {x\-} + a 2 tf{x 2 , , 
2m 16 \ p A p 4 p° 

(14a 2 + 1680)ft 6 3 3 1 



+ fah i {x i ,-}V)e 2 + 
p° p° 



h 6 (a 3 (H{x 2 , ^} + ) + */3 3 {x\ ^} e 3 + G(e 4 ), (169) 

where a 2 , «3, fa, are free real parameters characterizing the nonuniqueness of the metric operator, and V is the 
parity operator [162]. 

A remarkable feature of the Hamiltonian (167) is that the underlying classical Hamiltonian is independent of the 
choice of the metric operator (to all orders of perturbation). Up to terms of order e 3 it is given by the following 
simple expression. 

The presence of p\ in the denominator of the second term is a clear indication that the equivalent Hermitian 
Hamiltonian is a nonlocal operator (and that this is the case regardless of the choice of the metric operator.) Again 
the classical Hamiltonian (170) clarifies the meaning of the imaginary cubic potential iex 3 . 



4.2.3 Other Models 

The perturbation theory usually leads to an infinite series expansion for the metric operator whose convergence 
behavior is difficult to examine. There are however very special models for which this method gives exact expressions 
for Q and consequently the metric operator i] + = . Examples of such models are given in [35, 48, 120, 49]. The 
simplest example is the free particle Hamiltonian studied in [162]. 

For other examples of the perturbative calculation of a metric operator and the corresponding equivalent Hermitian 
Hamiltonian, see [22] and particularly [88]. 



4.3 Differential Representations of Pseudo-Hermiticity 

In the preceding subsection we show how one can turn the operator equations appearing in the perturbativc calculation 
of the metric operator into certain differential equations. In this subsection we outline a direct application of 
differential equations in the computation of pseudo- metric operators for a large class of pseudo-Hermitian Hamiltonian 
operators H acting in the reference Hilbert space L 2 (K). 

In the following we outline two different methods of identifying a differential representation of the pscudo- 
Herniiticity condition, 

H^tjHt]- 1 . (171) 



33 



4.3.1 Field equation for the metric from Moyal product 

Consider expressing (171) in the form 

r]H = H^r), (172) 
and viewing rj and H as complex- valued functions of x and p that are composed by the Moyal *-product: 

f(x,p) * g{x lP ) := f(x,p) e iM * d » g(x,p) = £ £J- [d k J{x,p)\ d k p g{x,p). (173) 

fc=0 

This yields: r](x,p) * H(x,p) = H(x 7 p)* * r](x,p), [206]. With the help of (173) we can express this equation more 
explicitly as 

E ntr {^ H ^ d * - ifiH&py] d v} i&p) = °- ( 174 ) 

fe=0 

This is a linear homogeneous partial differential equation of finite order only if H(x,p) is a polynomial in x and p. 

2 

For example, for the imaginary cubic potential, i.e., the Hamiltonian H = 2— + iex , it reads 

[eft 3 ^ - 3ie^ 2 x92 _ (2m)' 1 h 2 d 2 x - 3 e hx 2 d p + im^hpdx + 2i e x 3 ] r)(x,p) = 0. (175) 

The presence of variable coefficients in this equation is an indication that it is not exactly solvable. Particular 
perturbative solutions can however be constructed. This applies more generally for other polynomial Hamiltonians. 
Explicit examples are given in [206, 89, 15]. 

We should like to note however that not every solution of (174) defines a pseudo- metric (respectively metric) 
operator. We need to find solutions that correspond to Hermitian (respectively positive-definite) and invertible 
operators rj. Ref. [206] suggests ways to address this problem. 

The above-described method that is based on the Moyal product has two important shortcomings. 

1. If the Hamiltonian is not a polynomial of x and p, then the resulting equation (174) is not a differential equation 
with a finite order. This makes its solution extremely difficult. This is true unless H(x,p) has a particularly 
simple form. A typical example is the exponential potential e %x treated in Ref. [69]. This potential is actually 
one of the oldest 'PT-symmetric potentials whose spectral problem has been examined thoroughly [91]. For 
x £ R, its spectrum includes an infinity of spectral singularities that prevent this potential from defining a 
genuine unitary evolution. 47 

2. For the polynomial Hamiltonians, for which (174) is a differential equation, the general form and even the order 
of this equation depends on the Hamiltonian. In particular, for the standard Hamiltonians of the form 

H=^+v(x), (176) 

they depend on the choice of the potential v(x). 

We shall next discuss a differential representation of the pscudo-Hermiticity that does not suffer from any of these 
shortcomings. 



4.3.2 Universal field equation for the metric 

Consider a pseudo- Hermitian Hamiltonian of standard form (176). Substituting (176) in the pseudo- Hermiticity 
relation r\ H = Wr] and evaluating the matrix elements of both sides of the resulting equation in the coordinate basis 
{|a;)}, we find [163] 

[-d 2 x + d 2 + f i 2 (x,y)] V (x,y)=0, (177) 

where n 2 {x,y) := %r [ v {x)* —v(y)}, and r)(x,y) := {x\rj\y). Eq. (177) is actually a Klein-Gordon equation for (x\r]\y) 
(with a variable mass term.) As such, it is much easier to handle than the equation obtained for the pseudo-metric 
in the preceding subsection, i.e., (174). Moreover, it applies to arbitrary polynomial and non-polynomial potentials. 

It turns out that one can actually obtain a formal series expansion for the most general solution of (177) that 
satisfies the Hermiticity condition: rj(x,y) = rjfy^x)*. This solution has the form [163] 

oo 

V (x,y)=J2^(^y)> (178) 

1=0 

47 For a discussion of biorthonormal systems for this potential with x taking values on a circle (a closed interval with periodic boundary 
condition on the eigenfunctions), see [68]. 
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where K. is the integral operator defined by 



m 



dr 



x+y- 



x—y-\-r 



ds v(r) f(s, r) 



ds 



x+y- 



dr v{s)* /(s, r) 



x+y+s 



(179) 



/ : R 2 — > C is an arbitrary test function, u : M 2 — > C is defined by u{x,y) := u+{x — y) + u-{x + y), and u± are 
arbitrary complex- valued (piecewise) smooth (generalized) functions satisfying u±(x)* = u±(^fx). 

For imaginary potentials, the series solution (178) provides an extremely effective perturbative method for the 
construction of the most general metric operator. For example, the application of this method for the "PT-symmetric 
square well potential that we discussed in Subsection 4.1.1 yields, after a page-long straightforward calculation, [163] 



r)(x,y) = 5(x-y) + ( 



w+(x 



y) + w-(x + y) + -^\x + y\ sgn(x 



V) 



(180) 



where w± : [— |f, -jf] — > C are arbitrary functions satisfying w±(x)* = w±(^fx) and w±(±L) = 0. The metric operator 
associated with the CPT- inner product that is computed using the spectral method in [51] turns out to correspond 
to a particular choice for w± in (180). 

A probably better evidence of the effectiveness of this method is its application in the construction of a metric 
operator for the PT-symmetric barrier potential that we examined in Subsection 4.1.2. Again, a two-pages-long 
calculation yields [163] 



r](x,y) = 5{x-y) + ([w + (x-y)+w-{x + y) + 



^ (2\x + y\-\ x + y + L\-\x + y-L\) sgn(x - y)] + 0(( 2 ). 



(181) 



This expression reproduces the result obtained in [159] using the spectral method (after over a hundred pages of 
calculations), namely (142), as a special case. 

In Ref. [163] this differential representation of pseudo-Hermiticity has been used to obtain a perturbative expres- 
sion for the metric operator associated with the imaginary delta-function potentials of the form 48 



v(x) 



N 

! E 

n=l 



(n 5{X 



where Cn,a n € K. The result is 



N 



v(x,y) = S(x -y) + ^2 



2m£ l; 



h 2 



w n+ (x -y) + w n -(x + y) + -0(x + y- 2a n ) sgn(y - x) 



o(C). 



(182) 



(183) 



For the special case: N = 2, a\ = —a2 > and Ci = — C2 > 0, where (182) is a PT-symmetric potential, a careful 
application of the spectral method yields a positive-definite perturbatively bounded metric operator [24] that turns 
out to be a special case of (183). The general N — 2 case, that depending on the choice of Cfc may or may not posses 
PT-symmetry, has been examined in [177, 138]. 

The main difficulty with the approaches presented in this section (and its subsections) is that they may lead 
to a "metric" operator that is unbounded or non-invertible. 49 For example setting N = 1 in (182), one finds a 
delta function potential with an imaginary coupling that gives rise to a spectral singularity [165]. Therefore, the 
corresponding Hamiltonian is not quasi-Hermitian, and there is actually no genuine (bounded, invertible, positive- 
definite) metric operator for this potential. Yet, one can use (183) to obtain a formula for a "metric operator"! This 
observation suggests that one must employ this method with extra care. 



4.4 Lie Algebraic Method 

In Subsection 4.2, we described a perturbative scheme for solving the pseudo-Hermiticity relation, 

H^=e- Q He Q , (184) 

for the operator Q that yields a metric operator upon exponentiation, ?] + = . In this section we explore a class 
of quasi-Hermitian Hamiltonians and corresponding metric operators for which (184) reduces to a finite system of 
numerical equations, although the Hilbert space is infinite-dimensional. The key idea is the use of an underlying Lie 
algebra. In order to describe this method we first recall some basic facts about Lie algebras and their representations. 

48 The spectral properties of PT-symmetric potentials of this form and their consequences have been studied in [117, 1, 7, 72, 226]. In 
particular, see [177]. 

49 One must also restrict the free functions appearing in the formula for rj(x, y) so that the operator 77 they define is at least densely- 
defined. 
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4.4.1 Lie algebras and their representations 



Consider a matrix Lie group G, i.e., a subgroup of the general linear group GL(N, C) for some N £ Z + , and let Q 
denote its Lie algebra, [83, 108]. A unitary representation of G is a mapping U of G into the group of all unitary 
operators acting in a separable Hilbert space T-L such that the identity element of G is mapped to the identity 
operator acting in H and for all gi, <?2 G G, U{g\92) = U{g\)U{g-i). Such a unitary representation induces a unitary 
representation for Q, i.e., a linear mapping 11 of Q into the set of anti-Hcrmitian linear operators acting in % such 
that for all X l7 X 2 G 5, it([X 1 ,X 2 ]) = [il(Jfi),il(X 2 )], [83, 85]. The mappings U and it are related according to: 
U(e x ) = e^ x \ for all X £ Q. Furthermore, because for all X £ Q, il(X) is an anti-Hcrmitian operator acting in H, 
there is a Hermitian operator K : H — >• 'H such that = «/v . 

Let {Ki, if 2, ■ • • , -fQ} be a set of Hermitian operators acting in H such that {iK\,iK2, • • • , is a basis of H(Q)- 
Then K a with a £ {1, 2, • ■ • , e?} arc called generators of G in the representation If U is a faithful representation, 
i.e., it is a one-to-one mapping, the same holds for il and d coincides with the dimension of G. In this case, we refer 
to the matrices 

Ka^li-'iKa) (185) 

as generators of G in its standard representation. 

Next, consider the set of complex linear combinations of K a , i.e., the complexification of Q: Q € :~ |X)a=i c a K a | c a G 

We can extend the domain of definition of il to Q c by linearity: For all c a G C, il fX)a=i c a ^k) := Ef=i c o ^(j^a) = 

So=i Ca K a - Similarly, we extend the definition of U to the set of elements of GL(N, C) that are obtained by 
exponentiation of those of Q c . This is done according to 

U{e x ) = e u(x \ iora\\Xeg c . (186) 

Next, we recall that according to Backcr-Campbcll-Hausdorff identity (150), for all X, Y £ Qc, e~ x Ye x £ Q c , 
Furthermore, (150) and (186) imply 

U{e- x Ye x ) = e-"Wll(F) =U(e- x )tt(Y)U(e x ) for all X, Y £ Q c . (187) 

This completes our mathematical digression. 



4.4.2 General outline of the method 

Suppose that H : % — > T-L can be expressed as a polynomial in the Hermitian generates K a of G in a faithful unitary 
representation U, 50 i.e., 

n d 

H = A aija2j ... Qfc K ai K a . 2 ■ ■ ■ K ak , (188) 

where n £ Z + , d is the dimension of G, and A aij(l2j ... a|f £ C. Demand that H admits a metric operator of the form 
i] + = e~Q with Q given by 

d 

Q = ^r a K a , for some r a £ R. (189) 

a=l 

Then as we will show below the right-hand side of (184) can be evaluated using the standard representation of Q 
and readily expressed as a polynomial in K a with the same order as H. Upon imposing (184), we therefore obtain a 
(finite) set of numerical equations involving the coupling constants A 0l)02i ... 0fe and the parameters r a that determine 
the metric operator via 



exp ^- ^2 r a K aj 



(190) 



In order to demonstrate how this method works, we introduce the matrix T] + := cxp (^— X)a=i r a K a J, that 
belongs to exp(£ c ) and satisfies ry + = U(i] + ). This together with (185) and (187) imply 

n d 

V + Hi 1+ - X = V,o 3 ,-o* {v + K ai V + ~ 1 )(v + K a2 i 1+ - 1 )---(ri + K ak r/ + - 1 ) 

k—1 ai,02}'" — 1 
n d 

= E E V, 02 ,.. g , ii(y + K ai ■ ■■ii(y + K ak (191) 

k—1 ai ,a2 ,ufc — 1 



50 In mathematical terms, one says that H is an element of the enveloping algebra of Q in the representation U. 
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Because for all a £ {1,2, • • • ,d}, rj + K a rj + 1 belongs to Q c , there are complex coefficients n a b depending on the 
structure constants C a bc of the Lie algebra Q and the coefficients r a such that 



V 



Kaf^- 1 = (192) 



6=1 



As a result, it(f? + KaV+ 1 ) = S&=i K ab il(Kb) = StLi K abKb- Inserting this relation in (191), we find 

n d 

V+H V^ 1 = E E \ bl b 2 -b k K bl K b2 ---K bkl (193) 

k=l bi,b 2 ,— ,b k =l 

where Xbib 2 -b k '■= Z)oi,o 3 ,- ,o fc =i ^a I ,a 3 ,-a fc K<n6i'W>a " ' Note tnat tllc coefficients A6 1 6 3 ...t h depend on the 

parameters r a of the metric operator (190). 

In view of (188) and (193), the pseudo-Hcrmiticity relation H' = i] + Hrj + ~ 1 takes the form 

n d 

E E (K k c k . 1 ...c 2 c 1 -Kc 2 -..c k . 1 c k ) K Cl K C2 ---K Ck =0. (194) 

k=l ci,oa,— ,Ck_i,cjfc=l 

We can use the commutation relations for the generators if a , namely [if , Kb] = i So=i ^abcK c , to reorder the factors 
K Cl K C2 ■ ■ ■ K Ck and express the left-hand side of (194) as a sum of linearly independent operators. Consequently, 
the coefficients of this sum must identically vanish. This yields a system of equations for r a . In general this system 
is over-determined and a solution might not exist. However, there is a class of Hamiltonians of the form (188) for 
which this system has solutions. In this case, each solution determines a metric operator. 
For the particular case that n = 1, so that 

d 

H = Y, KK h , (195) 

o=l 

we may employ a more direct method of deriving the system of equations for r a . This is based on the observation 
that in this case we can obtain a representation of the pseudo-Hermiticity relation = ri + Hri + ~ 1 in Qc, namely 

Hl=V ± Hr }± - 1 , (196) 

where 51 

d d 

H:=Y,^aKa, M^:=J2KEa- (197) 

a—l a—1 

The matrix equation (196) is equivalent to a system of d complex equations for d real variables r a . Therefore, it is 
generally over-determined. 

We can use the above Lie algebraic method to compute the equivalent Hermitian Hamiltonian h for the quasi- 

Hcrmitian Hamiltonians of the form (188). In view of the definitions: h := pHp^ 1 and p := y/rf^ = exp ^X^a=i if > 

h is given by the right-hand side of (192) provided that we use % in place of r a . 

An alternative Lie algebraic approach of determining metric operator and the equivalent Hermitian Hamiltonian 
is the following. First, we use the argument leading to (191) to obtain 

n d 

h = P Hp- 1 =J2 E *»!.«».•••«» Wj>K»x P' 1 ) HP *** P' 1 ) ■ ■ -iXQ? K ak p- 1 ), (198) 

k—l Exi,a2,-" ,afc=l 

where p := exp f^f=i ^f-^-a^j- Then, we evaluate p K a p _1 and express it as a linear combination of K a with 
r a -dependent coefficients. 52 Substituting the result in (198) and using the linearity of il and K a = il(K a ) give 
h = YJk=i Z)oi,oa,- ,o fc =i £ a ua2 ,-a k K ai K a2 ■ ■ ■ K ak , where e ai ,a 2 , - a fc are r Q -dependent complex coefficients. In this 
approach, we obtain the desired system of equations for r a by demanding that h = ftJ . This is the root taken in 
[194] where the Lie algebraic method was originally used for the construction of the metric operators and equivalent 
Hermitian Hamiltonians for a class of quasi-Hcrmitian Hamiltonians of the linear form (195) with underlying su(l, 1) 
algebra. For an application of this approach to Hamiltonians that are quadratic polynomials in generators of SU (1, 1), 
see [17]. 

51 Unless G is a unitary group, U is not a *-representation [85], and ^ HJ . 
52 This is possible, because p 6 exp(5 c ). 
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4.4.3 Swanson Model: Q = su{l, 1) 

In this section we explore the application of the Lie algebraic method to construct metric operators for Swanson's 
Hamiltonian [224]: 

1 2 

H = hbj{a^a +-) + aa 2 + 13 a* , (199) 

where a := , X := x, p := ^= , a,/3,Qj,m are real parameters, m > 0, ui > 0, and fi, 2 oj 2 > 4a/?. The 

latter condition ensures the reality and discreteness of the spectrum of (199). 

The problem of finding metric operators for the Hamiltonian (199) is addressed in [118, 206, 182]. The use of the 
properties of Lie algebras for solving this problem was originally proposed in [194]. 

Swanson's Hamiltonian (199) is an example of a rather trivial class of quasi-Hcrmitian Hamiltonians of the 
standard form 

[p + A{x)} 2 

H = +v{x), (200) 

where A and v are respectively a complex- valued vector potential and a real- valued scalar potential, and M £ M + is the 
mass. It is easy to see that these Hamiltonians admit the cc-dependent metric operator: t] + = cxp (— | J dx ( <s[A(x)]) . 
This in turn yields the equivalent Hermitian Hamiltonian: h = ^- (p + $t[A(x)]) + v(x). The subclass of the 
Hamiltonians (200) corresponding to imaginary vector potentials has been considered in [3]. The Swanson Hamil- 
tonian (199) is a special case of the latter. It corresponds to the choice: M — - m , A(x) — i ^ ™"( a ~|) ^ x ^ anc [ 

v(x) = | ^ * -s^p ) mi ° 2x2 i where 

& --=h *-- = i- (201) 

As shown in [118, 182], the Hamiltonian (199) admits other exactly constructible metric operators. The Lie 
algebraic method considered in this section offers a systematic approach for constructing metric operators for this 
Hamiltonian. In order to describe the details of this construction we begin by recalling that the operators a and a) 
are the usual harmonic oscillator annihilation and creation operators that satisfy 

[a,a f ] = l. (202) 

A well-known consequence of this relation is the possibility of constructing a unitary representation of the Lie algebra 
su(l, 1) using quadratic polynomials in a and . To see this, consider the Hermitian operators: Ki := j(a 2 + a> ), 

K 2 := \{a 2 - a^), and K 3 := \(aa) + a^a) = \{a ] a + ±) that act in H := L 2 (R), [142]. In view of (202), they 
satisfy the su(l,l) algebra: [Ki,K 2 ] = —iK^, [K 2 ,K 3 ] = iK±, \K 3 ,Ki\ = %K 2 . Clearly the Hamiltonian (199) can 
be expressed as a linear combination of K\ , K 2 and K% : 

H = 2 [(a + P)K 1 + i(a - f3)K 2 + hw K 3 ] . (203) 

This relation identifies the Hamiltonian (203) as a special case of the Hamiltonians of the form (195) with G = 
SU(l, 1), d = 3, Ai = 2(a + /?), A 2 = 2i(a — (3), and A3 = 2hw. We can further simplify (203) by introducing 
non-Hcrmitian generators K± := K\ ± iK 2 . In terms of these, we have H = 2 (aK + + (3K^ + Hoj K3). 

In order to apply the above method of constructing a metric operator for (203), we need to find a set of generators 
K a of SU(1, 1) in its standard representation and a faithful unitary representation it of the Lie algebra su(l, 1) such 
that K a = il(K a ) for all a £ {1, 2, 3}. A simple choice is 

Kl--=2^ = - 2 [i oj' ^ : =2 a2 =2{-l oj' ^ := 2 ff3 = H0 (2 ° 4) 
where (7\, (J 2, and (T3 are the Pauli matrices. We also have 

K ± :=K 1 +iK 2 = (^° J J, K_:=K 1 -iK 1 =(^ . °J, (205) 

that fulfill U(K±) = K±. 

Comparing (204) and (205), we see that it is more convenient to work with the generators K± and K3 rather than 
K a with a £ {1,2,3}. This in particular suggests the following alternative parametrization of the metric operator 
(190). 

r) + = exp(zK + ) exp{2rK 3 ) exp( z * K-), with zeC.rel. (206) 
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Clearly, 

77+ = exp(z^+) exp(2rK3) exp(z* K_) = ^ ) , (207) 

= exp(-z*K-) exp(-2rK 3 ) exp(-zK + ) = ( l c -r z * e r -l~r\ z ^ ) ' ( 208 ) 
where we have made use of (204) and (205). Also in view of (197), (201), (203), and K± :~ K\ ±1X2 = K^, we have 

H = 2hu[aK++ f3K_+ K3} =hw ( J- (209) 

H* = 2hu[aK-+ 0K++ Ks\ = Hu> ( 2 J„ . (210) 

Next, we insert (207) - (210) in (196). This yields the following three independent complex equations that are 
more conveniently expressed in terms of s := e r and w :— e~ r z. 

s(as -w) + j3{w 2 - 1) + s\w\ 2 [w + as(\w\ 2 - 2)] = 0, (211) 

fi-as 2 + sw*(l + asw*) = 0, (212) 

i3w + sw*[w + as(\w\ 2 - 1)] = 0. (213) 

To solve these equations, we rewrite (213) as: w + as(\w\ 2 — 2) = -as — and use this relation in (211) to obtain 

as(l- \w\ 2 ) = -+w. (214) 



s 



Substituting this equation back into (213), we find j3(w — w*) = 0. Therefore, either j3 = or w £ R. It is easy to 
show using (212) that the condition $ = implies melas well. Hence, w is real and (214) reduces to a quadratic 
equation whose solution is 



-1 ± \ Aa 2 s 2 + 1 - AaB 

It turns out that (211) - (213) do not impose any further restriction on s. Therefore, (215) is the solution of the 



system (211) - (213). In terms of the original parameters r and z, it reads z = ( — 1 ± y 4a 2 e 2r + 1 — 4af3)/(2a) = 

(—hio± ^/4a 2 e 2r + huj — 4a/3)/(2a). Substituting this formula in (206), we find two one-parameter families of metric 
operators for Swanson's Hamiltonian. 

5 Systems Defined on a Complex Contour 
5.1 Spectral problems defined on a contour 

Consider the Schrddingcr operator + V(x), where V : R — > C is a complex- valued piecewise real-analytic 

,2 

potential. The study of the spectral problem for this operator and its complex generalization, — -^2 + ^( z )> with 
z taking values along a contour 53 T in C, predates the discovery of quantum mechanics. 54 The case of polynomial 
potentials have been studied thoroughly in [212]. For a more recent discussion, see [210]. 

The spectrum of + V{z) depends on the choice of the contour T. In the case that T visits the point at 

infinity, the spectrum is essentially determined by the boundary condition imposed on the solutions ^ : T — > C of 
the following eigenvalue equation at infinity. 

-4o+V(z)]*(z) = E*(z). (216) 

For an extended contour F, that is obtained by a continuous invertible deformation of the real axis in C, the 
eigenvalue problem (216) is well-posed provided that we demand ^(z) to decay exponentially as \z\ — > 00 along T. 
To make this condition more explicit, we identify T with the graph of a parameterized curve £ : R — > C in C, i.e., 

T = {C(s) I a G R }. (217) 



53 Here by the term 'contour' we mean (the graph of) a piecewise smooth simple curve that needs not be closed. 

54 Hermann Weyl's dissertation of 1909 provides a systematic approach to this problem. For a detailed discussion of Weyl's results, see 
[104, §10]. 
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The assumption that T is simple implies that £ is a one-to-one function, and we can express the above-mentioned 
boundary condition as 

l*(C(s))l exponentially as s -> ±00. (218) 



A simple consequence of this condition is 



\y(((s))\ 2 ds < 00. (219) 



If wc view C as a Riemannian manifold, namely M 2 endowed with the Euclidean metric tensor, and consider T 
as a submanifold of this manifold, we can use the embedding map £ : M — ¥ C to induce a metric tensor (g) on T. 
The corresponding line clement is given by d£ := ^J~Qds = ^dx(s) 2 + dy(s) 2 — |C'(s)| ds, where x(s) := 5ft(£(s)) and 
y(s) := 3(£(s)). Therefore, the integral measure defined by g on T is the arc-length element |C'( S )I ds. This in turn 
suggests the following paramctrization-invariant definition of the L 2 -inner product on T. 

/>oo 

^|<1>^:= / *( C ( s) )*$( C ( s ))| C '( s )| ds . (220) 



If we identify s with the arc-length parameter, for which \C{s)\ = 1, and let L 2 (T) := : T ->• C -< ^< 00 j, 
we can express (219) as 

# e i 2 (r). (221) 

This shows that the boundary condition (218) implies the square-intcgrability condition (221) along T. The converse 
is not generally true; the spectrum defined by the boundary condition (218) is a subset of the point spectrum 55 of the 
operator — 4-z + V(z) viewed as acting in the Hilbert space L 2 (T). In the following we shall use the term spectrum 
to mean the subset of the point spectrum that is defined by the boundary condition (218). 

In order to demonstrate the importance of the choice of the contour in dealing with the spectral problem (216), 
consider the imaginary cubic potential V(z) = iz 3 . As wc mentioned in Section 1, the spectrum defined by the 
boundary condition (218) along the real axis (£(s) = s) is discrete, real, and positive [79, 209]. But, the spectrum 
defined by the same boundary condition along the imaginary axis is empty. To see this, we parameterize the imaginary 
axis according to z = £(s) = is with s e M. Then the operator + iz 3 takes the form + s 3 , and we can 

respectively express the eigenvalue equation (216) and the boundary condition (218) as 



d 2 3 ' 

h 8 

ds 2 



^)=B#), (222) 



and 

\ip(s)\ — > exponentially as s — > ±00, (223) 

where ip(s) := ^S(is) for all s <E R. 56 But, it is well-known that (222) does not have any solution fulfilling (223) for 
either real or complex values of E. 57 

The imaginary cubic potential belongs to the class of potentials of the form 

V„(x) = \x 2 (ix) v , i/el, A e IR + . (224) 

As shown in [79, 209], for v > these potentials share the spectral properties of the imaginary cubic potential, if we 
impose the boundary condition (218) along a contour T„ that lies asymptotically in the union of the Stokes wedges 
[32]: 

St ■= {re-^' j r e [0,oo), <p G (-*„,*„) }, (225) 

where 



tf-2^=:',, 0-:=,-6 u , 6 v : =VTl . (226) 
Here by asymptotic inclusion of T„ in S~ we mean that if T v = {C^(s)\s £ K} for a piecewise smooth one-to-one 



55 By definition, the spectrum <r(A) of an operator A acting in a Banach space is the set of complex numbers E for which the operator 
A — EI is not invertible, i.e., one or more of the following conditions hold: (1) A — EI is not one-to-one; (2) A — EI is not onto; (3) 
A — EI is one-to-one so that it has an inverse, but the inverse is not a bounded operator [105] . The point spectrum of A is the subset of 
cr(A) consisting of the eigenvalues E of A, i.e., the numbers E for which A — EI is not one-to-one [195]. 

56 This was pointed out to me by Prof. Yavuz Nutku. 

57 The point spectrum of -j-j +s 3 is C, i.e., (222) admits square-integrable solutions for all E G C (This was pointed out to me by Prof. 
Patrick Dorey.) These solutions do not however satisfy (223). They do not represent physically acceptable bound states, because they 
do not belong to the domain of the observables such as position, momentum or some of their powers. 
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Figure 1: A contour Y v (the solid curve) lying asymptotically in the union of the Stokes wedges (the grey 
region) . The dashed lines are the bisectors of the Stokes wedges S„ . The angle d v between the positive 
real axis and the bisector of is also depicted. 



function : R — > C, then there must exist a positive integer M such that for all s <E R the condition ±s > M implies 
C(s) € Sf. Figure 1 shows Stokes wedges and a typical contour lying in 5" U 5+ asymptotically. 
For v = 2, (224) gives the wrong-sign quartic potential, 



V 2 (x) 



-Xx 4 



A e 



(227) 



which is known to have an empty spectrum (defined by (218) along the real axis.) Setting v = 2 in (225) and (226), 

we have = [re 10 r G [0,oo), 9 6 (-tt,-^) | and S} = |re lS r € [0,oo), 6» £ (-§,0) j. Therefore the 

condition that T2 must lie asymptotically inside U S 1 ^" excludes the real axis as a possible choice for T2. It is not 
difficult to sec that the same holds for all v > 2. 



5.2 Equivalent spectral problems defined on M. 

The fact that the spectrum of the potentials (224) defined by the above mentioned boundary condition along r„ is 
discrete, real and positive is by no means obvious, and its proof is quite complicated [79, 209]. In this subsection we 
will outline a transformation scheme that maps the spectral problem for these and similar potentials to an equivalent 
spectral problem that is defined on the real line [157]. This scheme provides an intuitive understanding of the spectral 
properties of the potentials (224) and in particular allows for a straightforward treatment of the wrong-sign quartic 
potential (227) that we shall consider in the following subsection. 

Given an extended contour T, we can use x := 5R(z) to parameterize it. We do this by setting V = {((x)\x £ IR} 
where ((x) := x + if(x) for all x <G R, and / : R — > R is a piecewise smooth function. This implies that along T, 
dz = d((x) = [1 + if'(x)]dx and the eigenvalue equation (216) takes the form 



-g{ X f^ +l g( X Yf"( X )±+i( X ) 

dx^ dx 



(228) 



where for all x € 



1 + if'(x) 



v(x) := V(x + if{x)), i>{x) := f (x + if{x)), 



(229) 



and a prime stands for the derivative of the corresponding function. 

Next, we examine the consequences of using the arc- length parametrization of T. If we define F : R — > R by 
F(x) := Jq y/l + f'(u) 2 du, we can express the arc-length parameter along T, which we denote by X, as X := F(x). 
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Under the transformation x — > X, the eigenvalue equation (228) takes the form 



-2t£(x) 



d 2 . , d ' 



^(x) = £;v(x), 



where 58 



C(x) 
«(x) 



tan" 1 (/'(£)) 



x=F- 1 (x) 



2ie(x; 



^F-^x)), 



a(x) :=f(x) 
V(x) ^^(F-^x)). 



/"(*) 



x=F- 1 (x') 



(230) 

(231) 
(232) 



The arc-length parametrization of V is achieved by the function G : M — >• T defined by G(x) := F _1 (x) + 
if(F~ 1 (x)) for all xel. This is the invertible function that maps the real line K onto the contour T in such a way 
that the (Euclidean) distance is preserved, i.e., if Xi,X 2 <G K are respectively mapped to z\ := G(Xi) and z 2 := G(x 2 ), 
then the length of the segment of T that lies between z\ and z 2 is given by |Xi — X 2 |. In other words, G : M — > T is 
an isometry. In light of (229) and (232), G relates the solutions VP and if> of the eigenvalue equations (216) and (230) 
according to ip(x) = ^>(G(x)). We can use G to express v in terms of the potential V directly: 



t>(x) = e 2i «( x V(G(x)). 



(233) 



Furthermore, recalling that the arc-length parametrization of T corresponds to identifying s and £(s) of (217) - (219) 
respectively with X and G(x), we have 



|*(G(x))| 2 dx 



|V(x)| 2 dxHV#)- 



(234) 



This observation has two important consequences. Firstly, it implies that ^ € -^ 2 (r) if and only if ip <G L 2 (M.). 
Secondly, it allows for the introduction of an induced unitary operator G* : L 2 (T) — > L 2 (R), namely G*^) := ip if 
V>(x) = *(G(x)) for all X e M. 59 

The above constructions show that a pseudo-Hermitian quantum system that is defined by a Hamiltonian operator 



of the form 



Hamiltonian operator 



dz 1 



V(z) acting in the reference Hilbert space L (T) is unitary-equivalent to the one defined by the 



H 



3 -2if(x) 



d 2 . , , d . ' 
— + ia(x)-+v(x) 



In terms of the unitary operator G* , we have 



dr 

dz' 2 



(235) 
V(z) = 



that is defined in the reference Hilbert space L 2 (B 
G~^ l H G„ . 

A particularly simple choice for a contour is a wedge-shaped contour: := {x + if(x)\x <G R} where /(x) := 
— tan# and € [0, ^). A typical example is the contour obtained by adjoining the bisectors of the Stokes wedges 
(the dashed lines in Figure 1.) For such a contour, 



x = F(x) = sec 9 x, x = F 1 (x) 
G(x) = cos 9 x - i sin 9 |x| = 



= cos 9 x, 

*C for X < 0, 
X for x > 0, 

and in view of (231), (232), and (233) the Hamiltonian operator (235) takes the form 

d 



H 



= 2i0 S gn(x) 



dx 2 



2i9S{x) 



dx 



V(cos9x- i sin 9 \x\ 



(236) 
(237) 

(238) 



The presence of delta-function in (238) has its root in the non-diffcrcntiability of at the origin. One can 
smooth out T^ 6 ' in an arbitrarily small open neighborhood of the origin and show that this delta-function singularity 
amounts to the imposition of a particular matching condition at X = for the solutions of the corresponding eigenvalue 
problem. As shown in [157], these are given by 



where for every function 



^(o+) = v(o) = v«r), < 

.->K, 0(O±) :=lim x ^ 0± <j>(x). 



V(o-) 



V(o+), 



(239) 



58 Because F is a monotonically increasing function, it is one-to-one. In particular it have an inverse that we denote by F _1 . 
59 (234) shows that G* preserves the norm. This is sufficient to conclude that it is a unitary operator, for in an inner product space the 
norm uniquely determines the inner product [123, §6.1]. 
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In view of (238), we can express the eigenvalue equation for H in the form 



where 



and ?/>_ : (—00, 0] 
44(0) :=^(0±). 



H±ip±(x) =Ei/)±(x), for xel 1 , (240) 

H± := - e ±2ie -^ + V(e Tl9 x), (241) 
-0+ : [0,oo) -> C are defined by ^±(x) := ip(x) for all x g R±, V±(0) := -0(0*), and 



In summary, the eigenvalue problem for the Schrodinger operator —-j^i + V( z ) that is defined by the boundary 
condition (218) along is equivalent to finding a pair of functions ip± satisfying 



,±2i0„/,// 



V>±(x) + V(e^x) V±(x) =Eif>±(x) for xe 



— e 

0_(O) = 0+(O), e- 2l V-(0) = e 2i V;(0), 

■0±( x ) ~~ ^ exponentially as X — > ±00. 



(242) 
(243) 
(244) 



To elucidate the practical advantage of this formulation, we explore its application for the potentials V v {z) = A z 2 (iz) 1 ' 
with A g R+. 

As we explained in the preceding subsection, we need to choose a contour that belongs to the union of the Stokes 
wedges asymptotically. We shall choose the wedge-shaped contour T^ 6 '^ that consists of the bisectors of S^. 
Setting V = V U and 9 = 6 V in (241) and using (226) we find the following most surprising result. 



Hi 



+ 2i6 



dx 2 



A|x 



<y+2 



Similarly, (242) becomes 



-^(x) + A|x| l/+ V±(x) = E e* 2ie "ip±(x) for xe 



(245) 



(246) 



The appearance of the real confining potential A|x| ly+2 in (245) and (246) allows for an alternative proof of the 
discreteness of the spectrum of the potentials z 2 (iz) v ■ See [157, Appendix] for details. 60 



5.3 Wrong-sign quartic potential 

In the preceding section we showed how one can transform the spectral problem for a potential defined along a 
complex contour to one defined along R. The form of the transformed Hamiltonian operator depends on the choice 
of the contour and its parametrization. This raises the natural question whether one can choose an appropriate 
parameterized contour so that the transformed Hamiltonian admits an easily constructible metric operator. The 
wrong-sign quartic potential V%{z) = — Az 4 , with A g R + , is a remarkable example for which the answer to this 
question is in the affirmative. 

Let T 2 = {((s)\s g R} be the contour defined by [120] 



C(s) := -2iy/T+ 



is, 



for all s g 



(247) 



If we parameterize T 2 using x = 9?(C(s)), we find T 2 = {x + if(x)\x g R} where / is given by f(x) := — \Jx 2 + 1 for 
all x g R. This shows that T 2 is a hyperbola with asymptotes l± := {±re~ 1 ^ \r g R+}. Because £± lies in the Stokes 
wedge S^, (247) defines an admissible contour for the potential V 2 . 



If we perform the change of variable z — > C( s )i we can express the eigenvalue equation: 
E$>(z), in the form 



ds 2 



1 d 

2 ds 



16A(1 + isf 



(s)=Ecj)( S ), 



(248) 



where 4>(s) := ^'(('(s)). We can identify s with a usual (Hcrmitian) position operator, introduce the corresponding 
wave number operator ^ as (s\& :- 



-i4-(s\, and express (248) as the eigenvalue equation for the Hamiltonian 



H := (1 + is)& 2 + - - 16A (1 + is) 2 



(249) 



that acts in L 2 (l 



60 It would be interesting to see if this approach can lead to an alternative proof of the reality of the spectrum. 
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As shown in [120], the application of the perturbative scheme of Subsection 4.2 yields the following exact expres- 
sions for a metric operator and the corresponding equivalent Hermitian Hamiltonian, respectively. 



h 



cxp 



sr 



v 48A 
.ft 4 ft 
64A 2 



6 

16As 2 . 



(250) 
(251) 



We shall offer an alternative derivation of these formulas momentarily. 

Because h is a Hermitian operator that is isospectral to H, the spectrum of H and consequently the operator 
— -4-z — Xz 4 defined along r 2 , is real. It is also easy to show that the common spectrum of all these operators is 
positive and discrete. To see this, we express h in its ^-representation where eigenvalue equation h& = E$> reads 



16A 



dfi 2 + 64A 



ft 



(252) 



and := = (2ir) x l 2 ds e tRs ^(s). The operator in the square bracket in (252) is a Schrbdinger 

operator with a confining quartic polynomial potential. Therefore its spectrum is positive and discrete [139]. 



We can also use the same approach to treat the quartic anharmonic oscillator, V(z) 
parametrization (247), we find the following generalizations of (249) - (252). 

ft 



,2^2 



H 



(1 + is)& + 



V+ = exp 



48A- (2+ 4A^ 



16\(l + isy -4u 2 (l + is) 

2 



, (£ 2 -4w 2 ) 2 R „, 2 

ft 
2 



-16A 



(fi - 4w 2 ) 



2 ^2 



d& 2 



64A 



Xz . Using the 

(253) 
(254) 

(255) 
(256) 



Therefore, by scaling the eigenvalues according to E — > -fE where 7 := 1/(16A), we can identify the spectrum of the 
operator 



dz 2 



+ lu z — Xz (defined along T 2 ) with that of the operator 



(257) 



that acts in L 2 (R), [120]. This observation has previously been made in [58]. As discussed in [11], the approach of 
[10] also leads to the same conclusion. 

Next, we outline an alternative and more straightforward method of constructing a metric operator for the 
Hamiltonian (253). 

If we separate the Hermitian and anti-Hcrmitian parts of H, we find 



H = 16As 2 + £ 2 - - - (16A + 4w 2 ) + l -{s, ft 2 - Aoj 2 - 32A}. 
We can combine the first and last terms on the right-hand side of this equation to express H in the form 



H = 16A 



s + 



i{R 2 - 4lu 2 - 32 A) 
32A 



2 ,2 



(# 2 - 4oj 2 ) 
64A 



2' 



(258) 



(259) 



As seen from this relation, the term responsible for the non-Hermiticity of H may be removed by a translation of s, 
namely 

i(ft 2 - 4uj 2 - 32A) 



32A 



(260) 



It is not difficult to see that such a translation is affected by a ^-dependent similarity transformation of the form 
s — > e'Ws e~ 9( ^. Recalling that for any analytic function g : M — > R, 



and comparing this equation with (260), we find 



-m = s _ 



ig'(fi), 



96A 



8A 



+ c. 



(261) 



(262) 
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Here c is an integration constant that we can set to zero without loss of generality. In view of (261), we can map H 
to a Hermitian Hamiltonian h according to 

H^h:= e*Wff e-»W = 16A.s 2 + ^ ~~ ^ - -. (263) 

64A 2 y ' 

This is precisely the equivalent Hermitian Hamiltonian given by (255). Moreover, comparing (263) with the defining 
relation for the equivalent Hermitian Hamiltonian, namely h := pHp^ 1 , and recalling that p := v / ?/ + , where rj + is a 
metric operator associated with the Hamiltonian H, we find rj + = e 29 ^. In light of (262), this coincides with the 
metric operator given by (254). 

Next, we wish to explore the underlying classical system for the pseudo-Hermitian quantum system defined by 
the potential V(z) = uj 2 z 2 — Xz 4 along T 2 . If we view this potential as an analytic continuation of V(x) = lu 2 x 2 — Ax 4 , 
with x denoting the Hermitian position operator, we should identify the Hamiltonian operator for the system with 

Hr 2 := — + n 2 Z 2 - AZ 4 , (264) 
2m 

where £1 € R and A £ M + are dimcnsionful coupling constants, and Z and Pz arc the dimcnsionful coordinate and 
momentum operators along T%. Using an arbitrary length scale £, we can introduce the corresponding dimcnsionlcss 
quantities: 

2m£ 6 A y/2m~£ 2 n iP z Z . „. 

A: =^^' w: =^^' P ' := 1T' Z '^T (265) 

In terms of these the eigenvalue equation Hy 2 ^> = £^f takes the form 



d 2 

u i 2 2 \ 4 

-dz~ 2+UJZ 



V(z)=EV(z), (266) 



where E := 2mZ 2 £/h 2 . 

Now, we are in a position to apply our earlier results. Setting z = £(s) := —H\f\ + is, we can identify (266) 
with the eigenvalue equation for the pseudo-Hermitian Hamiltonian (253). One might argue that because z and p z 
represent dimensionless position and momentum operators, the same should also hold for s and &, respectively. This 
suggests to identify the dimcnsionful position (x) and momentum (p) operators as 

x = a£s, p= — -, (267) 
at 

where a £ IR + is an arbitrary constant. In view of (255), (265), (259), and (267), we obtain the following expressions 
for the dimcnsionful pseudo-Hermitian and equivalent Hermitian Hamiltonians. 



H > ^JL = (p'-SmW) 2 _ hp_ igA 

2m£ 2 64Al 4 W 



i(£- 2 p 2 - 8mQ 2 - 64m£ 2 A 

£x+ - 

64mA 



(268) 



h' := ^=^- 8mW)2 -^ + 16A^, (269) 
2m^ 2 64A£ 4 Ami 

where £ := £/a. Note that unlike h' that only involves the length scale £, H' depends on both £ and £. The same is true 
for the metric operator 77+ . This shows that the pseudo-Hermitian quantum systems defined by the Hamiltonian H' 
and the metric operator rj + with different values of the parameter a are unitary-equivalent to a Hermitian quantum 
system that depends on a single length scale (£). The latter is not, however, fixed by the Hamiltonian (264) and the 
contour T 2 along which it is defined. 61 

Supposing that £ is independent of h, we can take h — > 0, x — > x c , and p — > p c in (269) to obtain the underlying 
classical Hamiltonian. The result is 

H , = (pl-*rn£W? _ wpkxl 
64A£ 4 

We can also introduce the pseudo-Hermitian position and momentum operators (55): 

X := rj-^Wj 2 - * + t( ^ 2 - 8W " 2 7 6W2A) , P = V^Wf = P- (271) 

64mA£ 



31 If f! / 0, we can choose a such that £ = fl/yX. 
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As expected, in terms of X and P, the Hamiltonian H' takes the form 



(P 2 -8m£ 2 n 2 ) 2 _ HP_ 
64A£ 4 Ami 



16A£ 2 X' 



(272) 



Furthermore, the ?y-|--pseudo-Hermitian canonical quantization of the classical Hamiltonian (270) yields (272) except 
for the linear term in P. 

The fact that the pseudo-Hermitian quantum systems defined by the Hamiltonian (264) depend on an arbitrary 
length scale has its root in our identification of s with the relevant dimensionlcss position operator. We will next 
consider an alternative approach that incorporates the spectral equivalence of the Hamiltonians (264) and (257). It is 
based on treating & as the appropriate dimensionless position operator. More specifically, it involves replacing (267) 
with x = fl£& and p = — hs/((3£), where /3 is an arbitrary dimensionless real parameter. If we set (3 := l/(4\/A) = 
(32mA)~ 1 / 2 £~ 3 h, we find the following ^-independent expressions for the equivalent Hermitian Hamiltonian b! and 
the underlying classical Hamiltonian H' c : 



<;2 



ti 
H' 



^+4A 
2m 



x 2 - 



f* 2 
4A 

+ 4A x 2 - — 
2m \ c 4A 



fry — x , 

V m 



(273) 
(274) 



Note, however, that H' still depends on I; 
i 



H' = 



2m 



p — i 



n 2 



, iVXx 2 = - 8VA£ 2 

2 V Vk 



4A 



x 2 - 



4A 



- h 



(275) 



This is also true for rj + . Again the quantum systems determined by H' and r] + with deferent values of £ are 
unitary-equivalent to a system that is independent of £. Similarly to our earlier analysis, we can obtain an £- 
independent expression for H' in terms of the pseudo-Hermitian position and momentum operators: X = x and 

P = p- i A /f (aVAx 2 - - 8\/A£ 2 ) . The result is H' = |^ + 4A ( X 2 - ) - H.f^X. 



6 Complex Classical Mechanics versus Pseudo-Hermitian QM 

6.1 Classical-Quantum Correspondence and Observables 

In Subsection 3.1 we outlined a procedure that assigns an underlying classical system for a given pseudo-Hermitian 
quantum system with reference Hilbert space L 2 (M. d ). According to this procedure, that we employed in Subsec- 
tions 4.2.1, 4.2.2 and 5.3, the classical Hamiltonian H c may be computed using (60). In other words, to obtain 
H Cl we replace the standard position and momentum operators with their classical counterparts in the expression 
for the equivalent Hermitian Hamiltonian h and evaluate its (H — > 0)-limit. We can quantize H c to yield h, if we 
employ the standard canonical quantization scheme. We obtain the pseudo-Hermitian Hamiltonian H , if we use the 
pseudo-Hermitian canonical quantization scheme (58). 63 

By definition, a classical observable O c is a real- valued function of the classical states (x c ,p c ), i.e., points of the 
phase space M. 2d . If we apply the usual (Hermitian) canonical quantization program, the operator associated with a 
classical observable O c (x c ,p c ) is given by o := O c (x,p) where x and p are the usual Hermitian position and momen- 
tum operators. If we apply the pseudo-Hermitian quantization, we find instead O := O c (X,P) where X and P are 
the pseudo-Hermitian position and momentum operators, respectively (55). The common feature of both these quan- 
tization schemes is that they replace the usual classical Poisson bracket, {A c , £? c }pb := EjLi ( y §i JL §jT~ ~ §§^ L §^j > 
of any pair of classical observables A c and B c with (ih) -1 times the commutator of the corresponding operators A 
and B, {A c , i? c }pB — > (ih)^ 1 ^, B]. Let us also recall that given a pseudo-Hermitian quantum system specified 
by the reference Hilbert space H, a quasi- Hermitian Hamiltonian operator H : % — > "H, and an associated metric 
operator rj + : % — > H, the observables of the system are by definition 77 + -pseudo-Hermitian operators O acting in %, 

O f = r] + m 1 . (276) 



62 Eq. (273) was previously obtained in [36] where the term in (273) is attributed to an anomaly in a path-integral quantization of the 

2 

complex classical Hamiltonian + f2 2 z 2 — Az* along the contour F2- For a more careful treatment of this problem, see [121]. 



Clearly this is true up to factor-ordering ambiguities/terms proportional to positive powers of h. 



46 



It is absolutely essential to note that such an operator acquires its physical meaning through the classical-to-quantum 
correspondence: 

O c — > O, (277) 

where O c is the classical observable corresponding to the operator O. For example in conventional quantum mechanics, 
we identify the operator p : L 2 (R) — > L 2 (R) defined by pip = —ihip', with the momentum of a particle moving on R, 
because p corresponds to the classical momentum p c of the underlying classical system. Without this correspondence, 
p is void of a physical meaning. It is merely a constant multiple of the derivative operator acting in a function space. 

The situation is not different in pseudo-Hermitian quantum mechanics. Again the (?7+-pseudo-Hermitian) oper- 
ators that represent observables derive their physical meaning from their classical counterparts through the pseudo- 
Hermitian version of the classical-to-quantum correspondence. This is also of the form (277), but the operator O 
is now selected from among the ?7 + -pscudo-Hermitian operators. As we discussed in Subsection 3.1, if o : H % 
denotes the Hermitian observable associated with a classical observable O c , the corresponding pseudo-Hermitian 
observable is given by Eq. (54), i.e., O := p~ 1 op where p := y/fj,- 

Whenever one deals with a symmetric, P7~-symmetric, diagonalizable Hamiltonian H with a real spectrum, one 
can define the observables as operators O fulfilling the condition [39] 

O t =CVTOCVT, (278) 

where 

T :=TO^T (279) 

stands for the transpose of O, and C, V, T are respectively the charge, parity, and time-reversal operators that are 
assumed to satisfy 

C 2 = V 2 = T 2 = I, [C, VT] = [P,T\ = [C, H] = 0. (280) 
As we explained in Subsection 3.4, we can relate C to an associated metric operator ry + according to 

C = ?/+ 1 'P. (281) 

Because C 2 = V 2 — I, we also have 

C = Vr\ + . (282) 

Inserting (279) in (278) and making use of (280), we obtain TO f T = CVT O CVT = VTCOCVT = TVCOCVT. 
In view of (280), (281), and (282), this relation is equivalent to = VCOCV = V(Vr) + )0 {j)+ x V)V = t] + Ot]- 1 . 
Therefore (278) implies the 77 + -pscudo-Hermiticity of O. The converse is also true for the cases that (278) can be 
applied consistently. This is actually not always the case. For example, the application of (278) for the Hamiltonian 
operator that commutes with CVT gives H T = H. Therefore, unlike the ?7 + -pseudo-Hermiticity conditions (276), 
(278) cannot be employed for non-symmetric Hamiltonians. This shows that (278) has a smaller domain of application 
than (276). 

Another advantage of the requirement of 77 + -pseudo-Hcrmiticity (276) over the condition (278) is that it makes 
the dynamical consistency of the definition of observables more transparent. Recall that the main motivation for 
the introduction of (278) in [39] is that its original variant [37, 38], namely the requirement of CPT-symmetry of 
O, i.e., O = CVTOCVT, conflicts with the Schrodinger time-evolution in the Hcisenberg picture; in general the 
Heisenberg-picture operators On(t) := e - ltH / h Oe ltH / h do not commute with CVT for t ^ 0, even if they do for t = 0, 
[153] . The reason why (278) does not suffer from this problem is that it is a restatement of the ?7 + -pseudo-Hermiticity 
of O. To see why the Heisenberg-picture operators satisfy the latter condition for all t, we first recall that because 
H is ?7 + -pseudo-Hermitian, it is a Hermitian operator acting in the physical Hilbert space 'Hphys defined by the 
inner product (•|-) rj+ - This implies that the time-evolution operator e - ltH / h is a unitary operator acting in H p h ys - 
Therefore, e -*H/h 0e itH/H acts in 

as a Hermitian operator for all t, i.e., it is 77 + -pseudo-Hcrmitian for all t. 
Alternatively, we could argue that because H is ?7 + -pseudo-Hermitian, e - ltH / h is r?+-pseudo-unitary [4, 5, 159] and 
e -itH/hQ e itH/h j g ? 7 + _pseudo-Hermitian. 

6.2 Complex Classical Systems and Compatible Poisson Brackets 

In their pioneering article [32], Bender and Boettcher perform an asymptotic analysis of the spectral properties of 
the complex potentials V v (z) = z 2 (iz) v that makes use of the complex WKB-approximation. This involves the study 
of a certain type of complex classical dynamical system Sbbm that Bender, Boettcher, and Meisenger (BBM) [33] 
identify with the underlying classical system for the quantum system defined by V v . This approach, which has been 
the focus of attention in a number of publications [186, 43, 46, 47], is fundamentally different from the prescription 
we used in Subsections 3.1 and 6.1 to associate a classical system S with a pseudo-Hermitian quantum system S. In 
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this subsection, we examine the structure of the complex classical system Sbbm- For simplicity we consider complex 
potentials V that depend on a single complex variable 3. 64 

According to [33] the dynamics of Sbbm is determined by the Newton's equation 



mi = -V'{l), 



(283) 



where m £ R + , each overdot stands for a time-derivative, and a prime marks the differentiation with respect to 3. 
We can express (283) as a pair of first order differential equations: 



777,3 = P, 



(284) 



This is the Hamiltonian formulation of the dynamics of Sbbm', introducing the complex Hamilton function f) 
+ V(i), we can express (284) as the following pair of Hamilton equations. 



(285) 



The variables 3 and p are the coordinates of the phase space of the system Cp which is as a set identical to C 2 and 
M 4 . This observation suggests that, similarly to the quantum systems defined along a complex contour, the complex 
classical system Sbbm might admit a formulation that involves real variables. This is actually quite straightforward. 
We can define the real variables 



x := K(3), 
H r := = 



»(3), 



2m 



V 



5R(p), g:=9f(p), V r :=X(y), V% := 9(V), 



+ V r (x,y), H i :=$l(t)) = ?Z+V i (x,y), 

m 



(286) 
(287) 



and use the well-known relations 



' JL 

< dx 



dp 



J dqi 



to turn the complex Hamilton equations 

(285) to a system of four real equations. 

It turns out that the resulting system of equations and consequently the complex Hamilton equations (285) are 
not consistent with the standard symplectic structure (Poisson bracket) on the phase space C 2 = K 4 , [68, 166]. To 
see this, let us also introduce W\ := x. W2 '■= p, w 3 :— y, and W4 := q. Then the standard Poisson bracket on R 4 
takes the form 



{A, 73} 



PB 



v (st) OA dB 



j,k=l 



where are the entries of the standard symplectic matrix 



(288) 



j(s t) 



V 



1 \ 

-10 

1 
0-10/ 



(289) 



and A and B are a pair of classical obscrvables (real- valued functions of Wj). 65 Recall that given a Hamilton function 
H on the four-dimensional phase space obtained by endowing IR 4 with the symplectic structure corresponding to 



the standard Poisson bracket (288), we can express the Hamilton equations in the form 



{wj,H} PB . If we 



express (288) in terms of the complex variables 3 and p, we find that {3, f)}ps = {p, t)}pB = 0, [166]. Therefore, it is 
impossible to formulate the dynamics defined by (285) using the standard symplectic structure on C 2 . 

This observation raises the problems of the existence, uniqueness, and classification of the symplectic structures on 
C 2 = K 4 that are compatible with the dynamical equations (285). Ref. [68] gives a family of dynamically compatible 
symplectic structures. Ref. [166] offers a complete classification of such structures. The most general compatible 
symplectic structure is defined by the following non-standard Poisson bracket 



j,k=l 



64 The use of complex phase-space variables in standard classical mechanics is an old idea [137]. 
present study is to consider complex configuration variables. 

65 S1 := fe=1 Jjk dwj A dwh is the standard symplectic form on R 4 , [137]. 



(290) 



See also [222]. The subject of the 
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where Jy are the entries of the symplectic matrix 



J 



/ 1 + c -a -d \ 

-(1 + c) -d -6 

a d -1 + c 

d 6 1-c 



(291) 



and a, 6, c, d are arbitrary real parameters satisfying c 2 + d 2 — ab + 1 1. Regardless of the choice of these parameters, 
we have 3 = -{{3, f)]} and p — §p, F)]}. A particularly, simple choice isa = 6 = c = d = that yields 



J=J :=- 



/ 10 \ 

-10 

-1 

\ 01 / 



(292) 



6.3 Real Description of a Complex Classical System 

Among the basic results of classical mechanics is the uniqueness theorem for symplectic structures on the phase 



p2,/ 



[13]. In order to explain the content of this theorem, first we recall that a symplectic structure on 



p2d 



space 

determined by the corresponding Poisson bracket. Choosing a system of coordinates Wj, we can express the latter 
in the form {A, B}j = Y^j h=x^jk§^: §§~i where J' is a real, antisymmetric, nonsingular 2d x 2d matrix. The 
above-mentioned theorem states that there is always a system of (so-called Darboux) coordinates in which {A, B} j 
takes the form of the standard Poisson bracket. Application of this theorem for the Poisson bracket (290) yields a 
description of the dynamics defined by the complex Hamiltonian f) in terms of a real Hamiltonian K . 

The construction of the Darboux coordinates associated with the most general symplectic matrix (291) is described 
in [166]. These coordinates take the following particularly simple form for a = b = c = d = 0. 



"1 := \/2wi = \2x, pi := \2~W2 = V2p, X2 '■= V2W4 = s/2 



q, P2 



= V2w 3 = s/2 



(293) 



These are precisely the phase space coordinates used in [234] to study the complex trajectories appearing in the 
scmiclassical treatment of the propagator for a quartic anharmonic oscillator. They are subsequently employed in 
the description of 'P7~-symmctric models [124]. 

The use of the coordinates (293) together with the assumption that the potential V is an analytic function, i.e., 

dx 



the Cauchy-Riemann conditions, 7^ — = ^ + ^ = 0, hold, lead to the following remarkable observations. 



The equivalent real Hamiltonian K that describes the dynamics is twice the real part of the complex Hamiltonian 
i), [234, 166], 



K 



pI 



2m 



2V r (^,^)=2H r 



The imaginary part of f), i.e. 



Hi = 



X2P1 
2m, 



(294) 
(295) 



is an integral of motion; {-Hi, K}pb = 0, [166] 



This implies that the dynamical system defined by the Hamilton equations (285) is equivalent to a classical Hamil- 
tonian system with phase space R 4 and the Hamiltonian K = 2id r . Furthermore, in view of Liouville's theorem on 
intcgrablc systems, because the phase space is four-dimensional and Hi is an integral of motion that is functionally 
independent of K, this system is completely integrable [13, 228]. 

The integral of motion Hi generates a certain class of transformations in the phase space that leave the dynamics 
invariant. The infinitesimal form of these symmetry transformations is given by [166]: 



Xx - 




- Sxx, 


6x1 


:= £{a-i,-ffi}pB 


X2 ~ 




~ 5X2, 


5x 2 


:= £,{x2,Hi} PB 


Pi - 




\-Spt, 


5p\ 


= £{Pl,Hi}j> B 


P2 - 


"> P2 - 


I" Sp2, 


5P2 


= £,{P2,H 1 }pb 



%X2 

' 2m' 
d 

« Q Xl V r[ 2 . 2 J' 
2m ' 



(296) 
(297) 
(298) 
(299) 
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where £ is an infinitesimal real parameter. 

The existence of these symmetry transformations is related to the fact that the system involves a first class 
constraint. Choosing a particular value C for Hi, i.e., imposing the constraint $ := Hi — C = 0, and moding 
out the above symmetry transformations by identifying each orbit of these transformations with a single point, we 
can construct a reduced dynamical system © that has a two-dimensional real phase space [166]. This procedure 
has been implemented for a class of monomial potentials in [217] where the above symmetry transformations have 
been examined in the Lagrangian formulation and the difficult issue of the quantization of these systems has been 
addressed in some detail. 66 In general the resulting quantum system depends on whether one imposes the constraint 
before or after the quantization. In the former case the prescription used to obtain the reduced classical system also 
affects the resulting quantum system. For the imaginary cubic potential that allows for a fairly detailed analysis, one 
obtains a variety of quantum systems [217]. But none of these coincides with the pseudo-Hermitian quantum system 
we studied in Subsection 4.2.2. 

In summary the identification of the complex classical system <Sbbm with the classical limit of the pseudo- 
Hermitian quantum system S that is defined by a complex potential V meets two serious difficulties. Firstly, while S 
has a single real degree of freedom (one-dimensional real configuration space and two-dimensional real phase space) , 
<5bbm has two real degrees of freedom (a four-dimensional real phase space). Secondly, the Hamilton equations (285) 
that determine the dynamics of 5bbm are not consistent with the standard symplectic structure on the phase space 
of iSbbm- This in particular means that under the naive correspondence 5bbm ~~ ^ the Hamilton equations (285) 
are not mapped to the Heiscnbcrg equations for S. These problems persist regardless of whether the Hamiltonian 
operator for S is defined on the real line or a complex contour. In fact, the study of the systems defined on a complex 
contour reveals another difficulty with the naive correspondence <Sbbm — ^ namely that while the definition of S 
requires making a proper choice for a contour, Sbbm is independent of such a choice. 67 

A proper treatment of the complex dynamical systems 5bbm requires the investigation of a dynamically compatible 
symplectic structure. Once such a structure is selected one can adopt a corresponding set of Darboux coordinates 
and use them to obtain a standard (real) description of <Sbbm- The use of the real description reveals the curious 
fact that this system admits an integral of motion that is functionally independent of the Hamiltonian. This has two 
important consequences: 

1. The system is completely integrable; 

2. The system has a first class constraint. 

The choice Hi = that is adopted in [43, 46, 47] is just one way of imposing the constraint. It corresponds 
to restricting the dynamics to a three-dimensional subspace of the phase space. On this subspace there act the 
symmetry transformations (296) - (299) that leave the dynamics invariant. Moding out these transformations, one 
finds a reduced dynamical system & with a two-dimensional real phase space. It is the latter that can, in principle, 
be related to the quantum system S. So far the existence and nature of such a relationship could not be ascertained. 
For the simple polynomial potentials that could be studied carefully, the various known ways of constructing and 
quantizing & lead to quantum systems that differ from S. 

Another important issue is that the above procedure of constructing complex dynamical systems <Sbbm an d the 
corresponding reduced systems & may be carried through for any complex analytic potential. But, not every such 
potential defines a unitary pseudo-Hermitian quantum system. A typical example is the exponential potential V(x) = 
e lKX that is defined on K. It is well-known that the spectrum of the Hamiltonian operator H = p 2 / (2m) + e e lKX , with 
m G IR + and e, k G K — {0}, includes spectral singularities [91]. This shows that this operator cannot be mapped 
to a Hcrmitian operator by a similarity transformation, i.e., it is not quasi- Hcrmitian. Hence H is not capable of 
defining a unitary quantum system. Although the potential V(x) — e lKX defines a complex dynamical system <Sbbm, 
it cannot be related to a unitary quantum system. 68 

We conclude this section by underlining a rather interesting parallelism between the quantum and classical 
mechanics of complex (analytic) potentials. Supposing that a complex (analytic) potential V defines a unitary pseudo- 
Hermitian quantum system, we showed in Section 3 that this system admits an equivalent Hermitian representation. 
The above discussion reveals a classical analogue of this equivalence; the complex dynamical system defined by 
the complex Hamiltonian H = p 2 /(2m) + V admits an equivalent description involving a real Hamiltonian. What 
differentiates the pseudo-Hermitian and Hermitian representations of quantum mechanics is the choice of the inner 
product (equivalently a metric operator) on the space of state-vectors. What differentiates the complex and real 

66 The application of this approach for some multi-dimensional systems is studied in [95]. 

67 The phase-space path integral formulation of the wrong-sign quartic potential denned on the contour (247) is consistent with the 
standard Hilbcrt-spacc formulation, if one restricts the complex classical Hamiltonian f) to the contour (247), [36, 121]. Whether 
imposing this restriction would lead to a particular reduced classical system that is identical with the classical system S defined by the 
Hamiltonian (274) is worthy of investigation. 

68 For a study of the classical dynamics generated by this exponential potential, see [68]. 
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representations of the classical mechanics is the choice of the symplectic structure (equivalently Poisson bracket) on 
the phase (state) space. Mathematically the equivalence of the pseudo-Hermitian and Hermitian representations of 
quantum mechanics stems from the uniqueness theorem for separable Hilbert spaces. The classical counterpart of 
this theorem that is responsible for the above-mentioned equivalence of the complex and real representations of the 
classical mechanics is the uniqueness theorem for symplectic manifolds diffeomorphic to R 2d . 

7 Time-Dependent Hamiltonians and Path-Integral Formulation 

7.1 Time-dependent Quasi-Hermitian Hamiltonians 

Time-dependent Hamiltonian operators arise in a variety of applications of conventional quantum mechanics. Their 
time-dependence does not cause any difficulties except that for the cases that the eigenvectors of the Hamiltonian arc 
time-dependent, the time-evolution operator takes the form of a time-ordered exponential involving the Hamiltonian 
[54]. 69 The situation is quite different when one deals with time-dependent quasi- Hermitian Hamiltonians. 70 As the 
following no-go theorem shows the observability of the Hamiltonian and the unitarity of the time-evolution put a 
severe restriction on the way a quasi- Hermitian Hamiltonian can depend on time [167]. 

Theorem 2: Let T G K + and for all t G [0, T], H(t) be a time- dependent quasi- Hermitian operator 
acting in a reference Hilbert space H. Suppose that H{t) serves as the Hamiltonian operator for 
a pseudo-Hermitian quantum system with physical Hilbert space H p h ys . If the time- evolution of the 
system, that is determined by the Schrodinger equation: ihip(t) = H(t)tp(t), is unitary and H(t) is 
an observable for all t <= [0,T], then the metric operator defining Hphys does not depend on time, i.e., 
there must exist a time-independent metric operator rj+ such that H(t) is rj + -pseudo-Hermitian for all 
t£[0,T\. 

Following [167] we call a time-dependent quasi- (pseudo-) Hermitian Hamiltonian admitting a time- independent 
metric (pseudo-metric) operator quasi- stationary. Theorem 2 states that in pseudo-Hermitian quantum mechanics 
we are bound to use quasi-stationary Hamiltonians. To demonstrate the severity of this restriction, consider two-level 
quantum systems where the Hamiltonian H{t) may be represented by a 2 x 2 complex matrix H(t) with possibly 
time-dependent entries. The requirement that H(t) is quasi-Hermitian implies that H(t) involves 6 independent 
real- valued functions (because its eigenvalues are real). The additional requirement that H(t) is quasi-stationary 
reduces this number to 4. 71 This is also the same as the maximum number of independent real-valued functions that 
a general time-dependent Hermitian Hamiltonian can include. 

A simple implication of Theorem 2 is that the inner product of the physical Hilbert space cannot depend on 
time, unless one defines the dynamics of the quantum system by an operator that is not observable or allows for 
nonunitary time-evolutions. In other words, insisting on observability of the Hamiltonian operator and requiring 
unitarity prohibit scenarios involving switching Hilbert spaces as proposed in [41]. 

7.2 Path-Integral Formulation of Pseudo-Hermitian QM 

Among the original motivations to consider "PT-symmetric quantum mechanical models is the potential applications 
of their relativistic and field theoretical generalizations [37, 38] in elementary particle physics. A necessary first 
step in trying to explore the relativistic and field theoretical generalizations of PT-symmetric or more generally 
pseudo-Hermitian QM is a careful examination of its path-integral formulation. In this section we use the approach 
of [169] to elucidate the role of the metric operator in the path-integral formulation of pseudo-Hermitian QM and 
demonstrate the equivalence of the latter with the path-integral formulation of the conventional QM. 

We shall first review the emergence of path integrals in dealing with a simple conventional (Hermitian) quantum 
system. This requires a brief discussion of the trace of a linear operator. 

Let L : H — >• H be a linear operator acting in a separable Hilbert space H with inner product (-|-). Then the 
trace of L is defined by 

TV 

tr(L):=J2(^m n ), (300) 
n=l 

69 A rather common misconception in dealing with time-dependent Hamiltonians is to think that the time-reversal operator T changes 
the sign of the time variable t in the Hamiltonian, i.e., H(t) — > TH(t)T = H(—t) which in view of the definition of T, namely for 
all ip 6 L 2 (R), (J~ij))(x) := if>(x)* , is generally false. Sec for example [237], where the author considers a trivial non-Hcrmitian time- 
dependent Hamiltonian that is obtained from a constant 'PT-symmctric Hamiltonian through a timc-dcpcndcnt point transformation and 
a time-reparamctrization. 

7c) Timc-dcpcndent quasi-Hermitian Hamiltonians arise naturally in the application of pseudo-Hermitian quantum mechanics in quantum 
cosmology [150, 151]. See also [156, 88]. 

71 This can be easily inferred from the results of [178]. 
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where is an arbitrary orthonormal basis of W, [195]. Obviously, for N = oo, the right-hand side of (300) may 
not converge, and tr(L) may not exist. 

Suppose that K and L are linear operators for which tx{KL) < oo. Then invoking the completeness relation for 
£ n and using Dirac's bra-ket notation, we can show that 



N N N 



tr(LK) = (UL\U)(U\K\£ n ) = J2 (U\K\£ n )(UL\U) = ^ MU) =tr(KL). (301) 

m,n— 1 m,n— 1 ra—1 

A simple implication of this identity is that the right-hand side of (300) is independent of the choice of the orthonormal 
basis {6i}- 72 

For a linear operator L acting in L 2 (R), we can use the position basis {\x)} to compute tr(L). To demonstrate 
how this is done, let {£„} be an orthonormal basis of L 2 (M.). Using (300), the completeness relation for \x) and £„, 
and Dirac's bra-ket notation, we have 

°° poo poo poo poo 00 

tr(L) = E/ / dxdx'(Z n \x)(x\L\x')(x'\Z n )= / dx dx' (x\L\x')Y,( x '\U(U^) 

n=Q J-oo J -00 J -00 J -00 n=Q 

^00 poo 

dxdx',(x\L\x'} (x'\x) = dx(x\L\x). (302) 



— 00 •/ — 00 



Now, consider a quantum system defined by the Hilbert space T~L = L 2 (R) and a Hcrmitian Hamiltonian H that 
is an analytic (or piccewisc analytic) function of the usual (Hcrmitian) position operator (x), momentum operator 
(p), and possibly time (t). The generating functional (also called partition function) for the n-point (correlation) 
functions of the system is given by 



( 1 1 

h 



Z[J] = tr I Texp i ~ I dt[H-Jx]\), (303) 



where "Texp" denotes the time-ordered exponential, t\ and are respectively the initial and final times for the 
evolution of the system that are taken to be —00 and 00 in the scattering setups used particularly in quantum field 
theory, and J stands for the (possibly time-dependent) coupling constant for the source terms Jx, [70, 231]. The 
latter is by definition an observable [75]. In view of the fact that x is also an observable, this implies that J must be 
real- valued. 

One can easily justify the condition of the observability of the source term by noting that Z[J 
compute the rt-point functions according to [70]: (%[x{ti)x(t2) ■ ■ ■ a;(r n )]) = ^~z{i\ 5j(Ti)Sj(rly--sj(T ) 



is used to 
, where 

j=o 

Ti, T2, • • • , t„ £ [ti, £2]) x(t) denotes the position operator in the Heisenberg picture, i.e., x(t) := U 7 (t, t\) x U J {t\ 1 r), 
and Uj is the time-evolution operator associated with the interacting system; for all ti , t2 € K, 



C/ 7 (ti,t 2 ) :=Xcxp j-^ 2 cft [H- Jx]|. (304) 

The n-point functions are essentially the expectation values of the observables 1[x(ti) ■ ■ ■ x(T n )] (in the ground state 
of the system if ti = —t\ — > 00.) Therefore, the observability of the source term Jx is linked to the observability of 
the Hcisenberg-picture position operators x(n). The latter is ensured by the unitarity of the time-evolution operator 
Uj{h, ti) and the observability of the Schrodinger-picturc position operator x. 
In view of (302), we can express the partition function (303) in the form 

Z[J] = j dx (x\1 exp| — — J dt [H — Jx]j \x) — J dx(x, t\\x, ti), (305) 

where 

\x,t) := C/ J (0,t) t |x), (306) 
are the (generalized) eigenfunctions of the Heisenberg-picture position operator x(t). In light of (304) and (306), we 
also have, for all x\, x% € R, (xi,ii|x2,t2) = (xi\U,(ti , ^2) |^2) = (xi|Texp | — j- dt [H — Jx]| \x 2 ). Computing 

72 To see this, let {Cn} be another orthonormal basis of H. Then as we described in Section 2.3, £ n are related to £ n by a unitary operator 
U : H -+ Ti, Cn = Ui n . This in turn implies £^=i<Cn|£Cn> = <t/£„ \LU^ n ) = T,n^n\U^(LU)U) = T,n=i(^\(LU)U^ n ) = 

5Z^=i(?n|^?n)i where we have used (301) and t/t/t = /. 
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this quantity as a phase-space path integral and substituting the result in (305), we find the following phase-space 
path integral expression for the generating functional. 



Z[J] = V(x)V(p)e-*ti? dt[H(x > p > t) - J( - t)x] . (307) 



If H is a quadratic polynomial in p, we can perform the momentum path integral in (307) and convert it into a 
configuration-space (Lagrangian) path integral. This yields 

Z\J\ = [ V(x) 5tl (308) 



where Lj(x,x;t) := xp — H(x,p;t) + J(t)x, and p is to be identified with its expression obtained by solving x = 
dH(x,p] t)/dp for p as a function of x, x, and t. 

Next, we consider the extension of the above constructions to a system defined by the reference Hilbert space 
% = L 2 (R), a metric operator rj + : i 2 (K) — > L 2 (R), and an r] + -pseudo- Hermitian Hamiltonian H : L 2 (R) — > i 2 (R) 
that is again a (piecewise) analytic function of x, p, and possibly t. According to the Theorem 2, in order for H 
to be an observable that generates a unitary time-evolution, it must be quasi-stationary. As discussed in [167], this 
puts a severe restriction on the form of allowed time-dependent Hamiltonians. 73 

We can certainly work in the Hermitian representation (H, h) of the system where h := pH p _1 (with p := y / J77/) 
is the equivalent Hermitian Hamiltonian (59). In this representation the generating functional has the form 



Z[J] =tr ^Xexpj-~jf dt[h-Jx]^j. (309) 

As we showed above this quantity admits a phase-space path integral expression. But even for the case that H is 
a quadratic polynomial in p, the equivalent Hermitian Hamiltonian h docs not share this property, and we cannot 
convert the right-hand side of (309) into a Lagrangian path integral in general. This provides a concrete motivation for 
the derivation of a the path-integral expression for the generating functional in the pscudo-Hcrmitian representation 
of the system, i.e., CH p h ys ,H). 

In [44] the authors use the expression (308) (with t% = —t\ — > oo) to perform a perturbative calculation of the 
generating functional and the one-point function for the "PT-symmetric cubic anharmonic oscillator (161). 74 As a 
result they find an imaginary value for the one-point function. This is simply because the one-point function they 
calculate corresponds to the ground state expectation value of the usual position operator that is indeed not an 
observable of the system. The physically meaningful generating functional is [169] 

Z[J] =ti v+ (Xexpj— ^ dt[H- JX]^j , (310) 

where for every linear operator K : % — > 

N N 

tr r;+ (K) := J2^n\Kij n ) v+ = J2(^n\v + K^„), (311) 

n—1 n—1 

{ipn} is an arbitrary orthonormal basis of "Hphys, and X is the 77,-pseudo-Hermitian position operator (55). The 
n-point functions generated by (310) correspond to the expectation value of time-ordered products of the physical 
position operators and the resulting numerical values are necessarily real. 

It is not difficult to show that tr I)+ = tr. To do this, first we recall that because p :~ y/rj : "Hphys — ► H is a 
unitary operator, it maps orthonormal bases of "Hphys onto orthonormal bases of %. In particular, := pV'n form 
an orthonormal basis of W. This together with p 2 = r) + , p^ = p, (311), (300), and (301) imply 

N N N 

tr v+ (K) = J2(i'n\p 2 K^n) = J2(P^n\pK^ n ) = £ (£„ \pKp~^ n ) = tripKp- 1 ) = tr(K). (312) 

n—1 n—1 n—1 

This relation allows us to express (310) in the form 



ti 



Z[J] =tr ( £exp{ -- / dt [H - JX] } ) . (313) 



73 In rclativistic field theories, H is obtained by integrating the Hamiltonian density over a space-like hypersurface. This makes the 
time-dependence of H quite arbitrary and renders the imposition of the condition of quasi-stationarity of H an extremely difficult task. 
74 See also [34]. 
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Next, we employ the definitions of h and A, namely h := pH p 1 and X := p 1 xp, and the fact that r] + and 

consequently p do not dependent on time, to establish Texp | — jr [h — Jx]| = p Texp j — j- [H — JX] \ p^ 1 . 

In view of this relation and (301) the right-hand sides of (309) and (313) coincide. This is another manifestation of 
the physical equivalence of Hermitian and pseudo-Hcrmitian representations of the system. 

As we emphasized in the preceding sections the metric operator plays a fundamental role in the operator formu- 
lation of pseudo-Hermitian quantum mechanics. The same is true about the path-integral formulation of this theory. 
To elucidate this point we examine the nature of the dependence of the generating functional on the choice of a 
metric operator r/ + . 

A simple consequence of (302) and (313) is 

Z[J]=J fx (xlTexpj-^ dt [H - JA]| \x). (314) 

Clearly, Z[0] does not depend on rj + , [169]. This explains the results of [113] pertaining the metric-independence of 
thermodynamical quantities associated with non-interacting pseudo-Hermitian statistical mechanical models. How- 
ever, in contrast to the view expressed in [122], the metric-independence of Z[0] does not extend to Z[J] with Z^O. 
This is actually to be expected because the knowledge of Z[J] allows for the calculation of the n-point functions 
that are expectation values of the time-ordered products of the Heisenberg-picture rj + -pseudo-Hermitian position 
operators Xfa). 

The dependence of Z [J] on the choice of rj + is rather implicit. In the Hermitian representation, i] + enters the 
expression for Z[J] through the equivalent Hermitian Hamiltonian h. In the pseudo-Hermitian representation, this 
is done through the source term JX. The presence of X in (314) prevents one from obtaining a Lagrangian path 
integral for Z[J] even for the cases that H is a quadratic polynomial in p. Therefore, contrary to claims made in [44], 
in general, the pseudo-Hermitian representation is not practically superior to the Hermitian representation. There 
are certain calculations that are performed more easily in the pseudo-Hcrmitian representation, and there are others 
that are more straightforward to carry out in the Hermitian representation [154]. 



8 Geometry of the State- Space and the Quantum Brachistochrone 

8.1 State-Space and Its geometry in Conventional QM 

In conventional quantum mechanics the states are not elements of the Hilbcrt space "H, but the rays (one-dimensional 
sub-spaces) of the Hilbert space. 75 The space of all rays that is usually called the projective Hilbert space and denoted 
by 'P('H) has the structure of a manifold. For an A-dimcnsional Hilbert space T~L, 'P('H) is the complex projective 
space CP^" 1 . This is a compact manifold for finite N and a well-known infinite-dimensional manifold with very 
special and useful mathematical properties for infinite N, [54]. 

The projective Hilbert space V{H) is usually endowed with a natural geometric structure that is of direct relevance 
to physical phenomena such as geometric phases [189] and optimal-speed unitary evolutions in quantum mechanics 
[9]. To describe this structure, we need an appropriate representation of the elements of 'P('H). This is provided by 
the projection operators associated with the states. 

Consider a state A^, represented by a state- vector ip G H — {0}, i.e., = {c0|c e C}, and the projection operator 

^ := 7v#T (315) 

Clearly the relation between states A^, and state-vectors cip is one to (infinitely) many. But the relation between the 
states A,/, and the projection operators A^, is one-to-one. This suggests using the latter to identify the elements of the 
projective Hilbert space 'P('H). This parametrization of V(J€) has the advantage of allowing us to use the algebraic 
properties of the projection operators (315) in the study of states. 

An important property of (315) is that it is a positive operator having a unit trace. The positivity of A^, is a 
simple consequence of the identities 

A;=A^, A2=A^. (316) 

We recall from Subsection 7.2 that the trace of a linear operator J : H — > H is defined by 

N 



tr(J) :=J2(UJU, (317) 



3 Throughout this article the word "state" is used to mean "pure state" 
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where is an arbitrary orthonormal basis of H, [195]. If L : H — >• H is a linear operator such that tr(L^i) < oo, 
L is said to be a Hilbert- Schmidt operator. In view of (315), (316), (317) and (19), we have 

tr(A^) - tr(A^) - tr(A^) - ^ ^ ^ ^ - 1. (318) 

Therefore, A^ is a Hilbert-Schmidt operator with unit trace. 

The set 032(H) of Hilbert-Schmidt operators forms a subspacc of the vector space of bounded linear operators 
acting in W. We can use "tr" to define the following inner product on 032(H)- 

(L\J) := tr(£t J) for all L,Je 03 2 (H). (319) 

This is called the Frobenius or Hilbert-Schmidt inner product [106, 195]. It has the appealing property that given 
an orthonormal set {xn} of state- vectors the corresponding set of projection operators {A Xii } forms an orthonormal 
subset of 03 2 (H); (xm\Xn) = S mn implies (A x JA Xm ) = S mn . 

The set i5 2 (H) of Hermitian Hilbert-Schmidt operators to which the projection operators A^ belong is a subset 
of 032 (H) that forms a real vector space with the usual addition of linear operators and their scalar multiplication. 
It is not difficult to see, with the help of (301), that (319) reduces to a real inner product on Sj 2 CH), namely 

(L\J) := tr(ZJ) for all L, J e £ 2 (H). (320) 

Therefore endowing ij 2 (H) with this inner product produces a real inner product space. 

By identifying states with the projection operators A^, we can view the projective Hilbert space V{T-L) as a 
subset of -f)2("H) and use the inner product (320) to endow 'P('H) with a natural metric tensor. The corresponding 
line element ds at A^ is given by 

- \ c«*m*) - 

where we have inserted a factor of \ to respect a mathematical convention and used (315), (317), (320), and the fact 
that is an orthonormal basis of %. 

For N < oo, we can identify ip with a nonzero complex column vector J* with components 31,32, ■ ■ ■ ,}n- In terms 
of these we can express (321) in the form ds 2 = J^o 6=1 9ab*d% a d%l where 

1 3*1 $ab ~ iabb /-oil 1 ! 
gab* ■■= j~4 • (322) 

This is precisely the Fubini-Study metric on the complex projective space CP^ -1 , [82]. 

As a concrete example, consider the case that H is two-dimensional, i.e., N = 2. Then using the standard basis 
representation of operators acting in C 2 , and noting that in this case all operators arc Hilbert-Schmidt, we can infer 
that i}2(H) is equivalent to the set of all Hermitian matrices. This is a four-dimensional real vector space which we 
can identify with R 4 . Specifically, we can represent each J £ .f)2(H) using its standard matrix representation: 

J = ( Z x-'y) (323) 
— V x + iy w J ' 

and observe that these matrices are in one-to-one correspondence with (x,y,z,w) G R 4 . 

The projective Hilbert space V(H) is a two-dimensional subset of the four-dimensional real vector space i5 2 (H)- 
To see this let us choose an arbitrary state-vector -0 6 C 2 — {0}. Then ijj = ^ ) for some 31,32 G C such that 
1 3 1 1 2 + 1 3 2 1 2 7^ 0, and in view of (315) we can represent A^, by 



^ - , l2 1 I2 ( Si' £51 ) ■ 



I3i| 2 + l32| 2 V 3l32 |32 
Using the parametrization (323), we find that for J = A^, 

3i32 +3i32 ^(3i32 ~3i32) l3i| 2 I32I 2 . , 

* 2(|3i| 2 + |3 2 | 2 )' y 2(|3i| 2 + |3 2 | 2 )' Z l3i| 2 + l32| 2 ' W l3i| 2 + l32| 2 " 1 ' 

Therefore, as expected w = 1 — z, so that 
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Figure 2: Stereo-graphic projection of the sphere S 2 defined by x 2 + y 2 + (z — i) 2 = j. N and S are 
respectively the north and south poles of S 2 . IIn is the tangent plane at N. p is a point on S 2 — {S}, 
and p is its stereographic projection on IIn- 



and the condition A 2 = takes the form 

x 2 +y 2 + (z- 1 -) 2 = \. (327) 

This defines a two-dimensional sphere S 2 that we can use to represent V(Ti). 

If we endow R 3 , that is parameterized by the Cartesian coordinates (x, y, z), with the Euclidean metric, we can 
identify S 2 with a round sphere of unit diameter. We will next obtain an expression for the metric induced on S 2 by 
the embedding Euclidean space R 3 . 

Let N and S respectively denote the north and south poles of S 2 , i.e., N := (x = 0, y = 0, z = 1) and S := (x = 
0,y = 0,z = 0), and consider the stereographic projection of S 2 onto the tangent plane IIn at N, as depicted in 
Figure 2. The line connecting S to a an arbitrary point p = (x, y, z) on S 2 — {S} intersects IIn at a point p. If we 
set up a Cartesian coordinate system in IIn with N as its origin and axes parallel to the x- and j/-axes, and denote 
by (x, y) the coordinates of p in this coordinates system, we can uniquely identify the points p G S 2 — {S} with 
(x, y) el 2 . Using simple methods of analytic geometry, we can easily verify that 



x y x y 1 

z' y ~z' x "l + x 2 + y 2 ' y ~l + x 2 + y 2 ' z ~l + x 2 + y 2 



We can employ the last three of these relations to compute the line element over the sphere in the (x, y)-coordinates. 
A rather lengthy but straightforward calculation yields 



dx 2 + dy 2 
(l + x 2 + y 2 ) 2 



ds 2 = dx 2 + dy 2 + dz 2 = Fl ■ ^ I . (329) 



This relation together with ds 2 = Y^,ij=i 9lj dlCdx? , where X 1 := X and x 2 := y, gives the following local coordinate 
expression for the Fubini-Study metric tensor [82], 

(FS) = Sjj 

13 [1 + (xi) 2 + (x 2 ) 2 ] 2 
Expressing X and y in terms of the spherical coordinates, 

^:=tan _1 (|), 9 := cos" 1 (2z - 1), (331) 
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of S 2 , we find 

X = tan (|) coscp, J = tan (§) sinip. (332) 
Substituting these in (329) leads to the following familiar relation for the line element of a sphere of unit diameter. 

ds 2 = - (d9 2 + sin 2 9 dip 2 ) . (333) 

In order to make the relationship between (321) and (329) more transparent, we recall that the south pole S of 
S 2 corresponds to the projection operator A e2 represented by 



and the state A e2 := {ce2 | c € C}, where e-i :=( J J . The coordinates (x, y) parameterize all the states except A 
These are represented by the state-vectors 



lf-2 ■ 



ip = ( ^ ) with 3i ^ 0. (335) 



Introducing Z := we can simplify the expression (324) for the corresponding projection operator. In terms of Z, 

3?(Z) _ 3(Z) _ 1 
r+|zp' » ~~ 1+|Z| 2 ' Z ~ TT|Zp 



the coordinates x, y, z appearing in (326) take the form x = m§p i 2/ = i+izl 3 ' z = T+Tzp' Comparing these with 



(328) reveals 

x = 5R(z), y = 3?(z)- (336) 

Now, we are in a position to compute the line element (321). Inserting (335) in (321), setting j 2 = 3i z j and using 
(336) and (332), we find 

, 9 dz*dz dx 2 + dy 2 1 / ,„9 9 „ , o\ , 

ds = 77 7^77 = 77 7-^-977 = 7 ( d0 + sin dl f ) ■ 337 

(l + |z| 2 ) 2 (l + x 2 + y 2 ) 2 4 V ; v ; 

Therefore, as a Ricmannian manifold the state-space V(H) is identical to a two-dimensional (round) sphere of unit 
diameter. 76 



8.2 State-Space and Its geometry in Pseudo-Hermitian QM 

The construction of the space of states in pseudo-Hermitian quantum mechanics is similar to that in conventional 
quantum mechanics. Again, the states are rays in the (reference) Hilbert space H which are identical with the rays 
in the physical Hilbert space "Hphys- The only difference is in the way one associates projection operators to states 
and defines an appropriate notion of distance (metric tensor) on the state-space. 

In the following we shall use 'P('Hphys) to denote the state-space of a pseudo-Hermitian quantum system with 
physical Hilbert space 'Hphys- The latter is obtained by endowing the underlying vector space of % with the inner 
product (-|-)r ;+ , where rj + : % — > % is a given metric operator rendering the Hamiltonian of the system ^-pseudo- 
Hermitian. We shall also introduce the following notation that will simplify some of the calculations: For all 

\(y.= |C> = C. <C\ ■■= (C\v + , -<t\C>--= (t\v + \0 = (t\v + = (£IC)r, + . 

First we define, for each pair of linear operators L, J : H — > H, 

(L\J) n+ := tr„ + {lKj) = tr(i« J), (338) 
where iJ stands for the r] + -pseudo-adjoint of L: 

LS^^-^V, (339) 

tr^ is defined by (311), and we have used (312). The linear operators A : Hphya — > "Hphys for which (A|A)^ + is 
finite together with (338) form the inner product space ^("Hphys) of Hilbert-Schmidt operators acting in "H p hys- 
Substituting (339) in (338) and using p 2 = 77 and = p, we also have 

(L\J) V+ = tT{ n+ - l L^ + J) = tr(p- 1 tfp 2 Jp- 1 ) = trdpLp-^pJp- 1 ) = (pLp-^pJp- 1 ). (340) 

76 The above calculation of the metric tensor on S 2 makes use of the stereographic projection of S 2 — {S} onto the plane I1n which 
is a copy of M 2 = C. We could also consider the stereographic projection of S 2 — {N} onto the tangent plane lis at S. Using both the 
projections we are able to describe all the points on S 2 . This is a manifestation of the manifold structure of S 2 . 
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These calculations show that p : H p h y s — > "H induces a unitary operator U p : ^("Hphys) — > 032(H) defined by 

U P {L) := P Lp-\ for all L £ B a (?W). (341) 

Now, consider a state A^ := {cip\c £ C} for some ip £ "Hphys — {0}. Because -< -|- >- is the inner product of H p hy S , 
the orthogonal projection operator onto A^ is given by 

A (r,+) = h/^^| = \WMv+ = A V- ^ ,o 42 n 
* • (i&M) (V>h+^> ' 1 ' 

(v ) 2 (»; ) 

A quick calculation shows that = . Furthermore, using the arguments leading to (318), we have 

tr fA^A^ -tr (\ {rl + } ) - ^1 ^V^t^H^ _ £n=l I | 2 _ 



This shows that A^ + ^ £ ^2 + \ and also, because tr^ + = tr, 



tr(Aj +) ) = l. (343) 

Another direct implication of (342) is A^' + ^ = r] + A,^' + ' i] + ~ 1 . Hence A? + belongs to the subset ^(Hphys) 
of r), -pseudo- Hermitian elements of S2 ("Hphys)- We can view this as a real vector space. (338) defines a real inner 
product on this space, and the restriction of (341) onto -^(Hphys) that we also denote by U p yields a unitary operator 
that maps -f)2(H p hy S ) onto PbCH)- In fact i^CHphys) and f)2(H) are real separable Hilbert spaces and the existence 
of the unitary operator U p : ifoCHphys) — > ^(H) is a manifestation of the fact that real separable Hilbert spaces of 
the same dimension are unitary-equivalent. 

The set of the projection operators (342) that is in one-to-one correspondence with the projective Hilbert space 
"P("H p h ys ) is a proper subset of -^(Hphys)- Similarly to the case of conventional quantum mechanics, we can define a 
natural metric on this space whose line element has the form 

del (A,) := VJ+X^) - -I I 2 



(344) 



It is important to note that as smooth manifolds V(H) and "P(H p hys) are the same, but as Riemannian manifolds 
they are different. While V(H) is endowed with the Fubini-Study metric, "P("H p hys) is endowed with the metric 
corresponding to (344). For N < 00 we can obtain a global expression for the latter in terms of the coordinates 
3i : 32 , ■ ■ • ,dN of the state-vectors if) = 3 £ C N . This yields the following generalization of (322) that satisfies 

ds 2 _y-N >+) d * 

(V+) Ec.fcltf+cii V+ba — V+ca V+bd) Icld m te\ 

9 ab , := — —2 • (345) 

W-ir,s=l V+rs iris I 

Here -q +a b are the entries of the standard representation -q + of rj + , [168]. 

For two-level systems where N = 2, we can easily obtain an explicit expression for the line element (344). In 
general the metric operator ij + is represented by 

/ a bi - ib 2 \ ,„._<. 
V±={b 1+ ib 2 c J' ^ 

where a, b\, 62, c € R are such that 

a + c = tr(i7 + ) > 0, d := ac - (6 2 + b\) = det(?y + )> 0. (347) 

In view of (342) we can again parameterize A^ + using the cartesian coordinates x, y, z of the sphere S 2 defined by 

(327). For the states differing from A e2 , we can alternatively choose the coordinates X and y of (336) and show using 
(335), (342), and (346) that 

dsl = d(rfx2 + dy2) (348) 

"+ [a + 2(& 1 X + & 2 y) + c(x 2 +y 2 )] 2 
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It is easy to see that for ij + = I where a = c = 1 and b 1 = b 2 = 0, (348) reproduces (337). 

To gain a more intuitive understanding of (348), we next express its right-hand side in terms of the spherical 
coordinates (331). Inserting (332) in (348) and carrying out the necessary calculations, we find 

2 fci(d6> 2 + sin 2 9 dp 2 ) _ fa (d8 2 + sin 2 9 dp 2 ) 

+ [1 + &2 cos# + g(p) sinf] [1 + k 2 cos 9 + k$ cos p> sin 9} 

where we have introduced 



d det(7jj a - c 2^bJTq 

ki:= Ja-Tc^ = ^J' k2: =—> fc3:= ^T^' (350) 

2frooB y + ferin y ) = ^ :=tan - 1( ^). (351) 

a + c 61 

Note that the change of coordinate tp —> (p corresponds to a constant rotation about the z-axis, and because of (347) 
we have k\ > 0, — 1 < k 2 < 1 and < k% < 1. 

The projective Hilbert space "P("H p hys) is the Riemannian manifold obtained by endowing the sphere S 2 with the 
metric g^+^ corresponding to the line element (349). 

Next, consider the general case where N need not be two. In order to compare the geometric structure of 
^(^phys) and V(H), we recall that 7 , ('Hphy S ) C S 2 ('Hphys) ! V(H) C < B 2 ('H), and that the linear operator U p of 
Eq. (341) maps < B 2 {'H P hys) onto ( B 2 (H). It is easy to check that the restriction of U p on 'P('Hphys), i.e., the function 
Up : V(Hp hys ) -> defined by 

u,(Aj +) ) := P I$ +) p-\ for all A^ +) G P(H phy8 ), (352) 

is a diffeomorphism. Furthermore, in view of (315), (342), and = p = -JffT, we have 



where * := pip G "H. A straightforward consequence of (321), (320), (340), (344), (352) and (353) is 

ds 2 {u p {K { ^ ] )) = ds 2 (A*) = l - (d^ |dA*) = I (pdA^p-'lpdA^p- 1 ) 

= i(dA ( ;+VA? +) ), + =d S) 2 )+ (A ( ; +) ). (354) 

This shows that u p : T'('Hphys) — > 'P('H) is an isometry, i.e., it leaves the distances invariant. Therefore, ^('Hphys) 
and V(7i) have the same geometric structure. In particular, for N = 2, we can identify 'P('Hphys) with a sphere of 
unit diameter embedded in R 3 with its standard geometry. 

This result is another manifestation of the fact that pseudo-Hermitian quantum mechanics is merely an alternative 
representation of the conventional quantum mechanics. Because the geometry of the state-space may be related to 
physical quantities such as geometric phases, we should not have expected to obtain different geometric structures 
for T'(Hphys) and T(U). 

8.3 Optimal-Speed Evolutions 

Let H : % — > H be a possibly time-dependent Hermitian Hamiltonian with a discrete spectrum. Suppose that we 
wish to use H to evolve an initial state X^j G V(H) into a final state X^ F G V(Ji) in some time r. Then the evolving 
state- vector ip(t) G T~L satisfies 

Hi>{t) = ihi>{t), *l>(Q)=1>i, 1>(t)=1> p , (355) 

and the corresponding state X^m traverses a curve in the projective Hilbert space V{'H). 
The instantaneous speed for the motion of A^,( t ) in V(J-L) is 

v., f , := (356) 



where ds is the line element given by (321). In view of this equation 

2 {m\m){m\m) -\mm 
u * (m\m) 2 v 



2 m)\m)(m\m) -ummw ^m? n , 7 , 

H = 7Z77u377ax2 = . ( 357 ) 
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where 

AE ^ t} : ~ (m\m) m)\m) 2 • (358) 

is the square of the energy uncertainty, and we have used (355) and the Hcrmiticity of H. We can employ (356) and 
(357) to express the length of the curve traced by X^u) m Viji) in the form [9]: 



AE^(t)dt. (359) 



Because AE^ is non-negative, s is a monotonically increasing function of r. This makes r a monotonically increasing 
function of s. Therefore, the shortest travel time is achieved for the paths of the shortest length, i.e., the geodesies 
ofV(H). 

For a time-independent Hamiltonian H we have ij)(t) = e~ ltH ^ h ipi and as seen from (358), AE^ is also time- 
independent. In this case, (359) yields 

T = xl/ (360) 

and the minimum possible travel time is achieved when s is identified with the geodesic distance between tpj and ipF- 
Because of the particular structure of V(H) the geodesic (s) connecting any two states A,/,, and A^ F lie entirely 
on a two-dimensional submanifold of V(H) that is actually the projective Hilbert space V(Hi,f) for the subspace 
Wi,f of "H spanned by ipi and ipF- If the time evolution generated by the Hamiltonian minimizes the travel time 
the evolving state should stay in V(Hi,f) during the evolution. This means that the problem of determining 
minimum-travel-time evolutions [90, 227, 135, 55, 64] reduces to the case that H is two-dimensional [56]. As we saw 
in Subsection 8.1, in this case 'P('H) is a round sphere of unit diameter and the geodesies are the large circles on this 
sphere. 

Next, we study, without loss of generality, the case N = 2. It is easy to show the existence of a time-independent 
Hamiltonian that evolves to X^, F along a geodesic. We will next construct such a Hamiltonian. 

Consider a time-independent Hamiltonian H acting in a two-dimensional Hilbert space. We can always assume 
that H has a vanishing trace so that its eigenvalues have opposite sign, E2 = —E\ =: E. 77 Because AE^, is time- 
independent we can compute it at t = 0. Expanding tjj(0) = tpi in an orthonormal basis {ipi, ^2} consisting of a pair 
of eigenvectors of H, i.e., writing it in the form 

tpi = Cllpl + C 2 -02, Ci, c 2 e C, (361) 

we find [1711 



12 _ U„|2 
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AE *' E r-{wT&?) iK <362) 

Therefore, the travel time r satisfies 

hs 

t > T min := — , (363) 

where s is the geodesic distance between A^, z and A^ in 'P('H). (363) identifies r m ; n with a lower bound on the travel 
time. Next, we shall construct a Hamiltonian H+ with eigenvalues ±E for which r = r m ; n . This will, in particular, 
identify T m ; n with the minimum travel time. 

In order to determine H+ we only need to construct a pair of its linearly independent eigenvectors ip\ and tp2 and 
use its spectral resolution: 

H it =E(-\^)(^\ + \rh)(H)- ( 364 ) 

As seen from (362), to saturate the lower bound on r, we must have |ci| = \c2\. In view of time-independence of 
AE^, we could also use ipF to compute this quantity. If we expand the latter as 

V'f = diipi + d 2 *p2, di,d 2 £ C, (365) 

we find that AE^p satisfies (362) with (ci, C2) replaced by (di, d 2 )- Therefore, by the same argument we find \d\\ = \d 2 \- 
Therefore, there must exist (3i,(3f £ R such that 

c 2 = e^ci, d 2 = e^di. (366) 



77 This is true for general possibly time-dependent Hamiltonians H(i). Under the gauge transformation ip[t) — > ip'{t) — > e ta ^/ n, 4){t) 
with a(t) = N^ 1 f tr[H(s)]ds, the Hamiltonian H(t) transforms into the traceless Hamiltonian H'(t) := H{t) — N ~ 1 tr[H (t)]I . 
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Inserting these relations in (361) and (365) gives 

V'i + e if> 'ih = c{ Vi, 01 + = CVf- (367) 

Solving these for ?/>i and 02 j we obtain 



1 - e* 



V'i = ; 3* . ^2 = ; , (368) 



where we have introduced 



0:=j8j-ftr, ^/ --ft-, ^-^jff- ( 369 ) 

V2ci \J ld\ 

Clearly i/> n determines the same state as n for n £ {1, 2}. Also in view of (364), the presence of /3p in (368) does 
not affect the expression for H+. The only parameter that has physical significance is the angle The orthonormality 
of ipi and ip2 implies that ipi and ipp have unit norm and more importantly that 1? fulfils 

$i\$f) = cos(|). (370) 

In terms of ipi and ipp this equation takes the form 

(^/|^) = ^di(l + e- i *). (371) 

Note that because tpi and tpF have unit norm, according to (369), {ipi\ipi) = 2|ci| 2 and (i/>f|0f) = 2|di| 2 . These 
relations together with (369) and (370) imply 



cos 2 (f ) = x . (372) 



Therefore whenever ipj and ipF are orthogonal, $ = tt. Furthermore, as discussed in [9], (372) shows that 1? is related 
to the geodesic distance s between A^ and X^, F according to 

i9 = 2s. (373) 

In the case that ipi and ipF are orthogonal, and A.^ F are antipodal points on V(7i). Therefore their geodesic 
distance s is half of the perimeter of a large circle. Because V(7i) is a round sphere of unit diameter, we have s = ir/2, 
which is consistent with (373). 

Having calculated the eigenvectors f/'i and '02, we can use (364) to obtain an explicit expression for the Hamiltonian 
H+ that evolves A^, 7 into A^, F in time T min . Substituting (368) in (364) and using (369), (370) and (373), we find 
[64, 56, 171] 

= iE{\^ F ){M-\i>i){j>F\) iEcot(s) ( \i> F )(^i\ _ ( n 74) 

4sin(|) 4 ^ ^j\i> F ) (Vf|0i) )' 1 ' 

The last equation shows that depends only on the states A,^ and A^, F and not on the particular state- vectors one 
uses to represent these states. Note also that this equation is valid generally; it applies for quantum systems with an 
arbitrary finite- or infinite-dimensional H. 

This completes our discussion of the quantum Brachistochronc problem in conventional quantum mechanics. We 
can use the same approach to address this problem within the framework of pseudo-Hcrmitian quantum mechanics 
[168]. This amounts to making the following substitutions in the above analysis: \ip n ) \ipn . >~, (^Pnl — > -<^n\, and 
s — > s, ;+ . In particular, the minimum travel time is given by 

= (375) 

and the 7y + -pseudo-Hermitian Hamiltonian that generates minimal-travel-time evolution between X^, I and A^, F has 
the form 

( V+ ) ^Ecot( Sv+ ) / typ^frl _ > — \ 

where 

>+' ^j^j^^f^F^ 



78 Note that the metric on V^H) that is used in [9] differs from our metric by a factor of v2. 



61 



Eq. (376) gives the expression for the most general time-independent optimal-speed quasi-Hcrmitian Hamiltonian 
operator that evolves 1^1 into tpp. Similarly to its Hcrmitian counterpart, it applies irrespective of the dimensionality 
of the Hilbert space. 

In [41], the authors show using a class of quasi- Hermitian Hamiltonians that one can evolve an initial state 
into a final state X^ F in a time t that violates the condition r > r m ; n . They actually show that by appropriately 
choosing the form of the quasi-Hcrmitian Hamiltonian one can make r arbitrarily small. This phenomenon can be 
easily explained using the above treatment of the problem. The minimum travel time for an rj -pseudo-Hcrmitian 
Hamiltonian is given by (375). In view of (349), depending on the value of k\ = det(r] + )/tr(r] + ) one can make s v 

and consequently r mi ^ as small as one wishes. This observation does not, however, seem to have any physically 
significant implications, because a physical process that involves evolving X^ )I into A^ F using an 77, -pseudo-Hcrmitian 
Hamiltonian H can be described equally well by considering the evolution of X p ^ I into A P ^ F using the Hcrmitian 
Hamiltonian h :— pH p _1 . In light of the existence of the isometry u p : 'P('Hphys) — > 'P('H), the distance between X^, 
and \,,p F in 'P('Hphys) is equal to the distance between A p ^, 7 and X P ^ F in V(H). It is also easy to show that the travel 
time for both the evolutions are identical. Therefore as far as the evolution speed is concerned there is no advantage 
of using the ?7 + -pseudo-Hermitian Hamiltonian H over the equivalent Hermitian Hamiltonian h, [168]. 

We wish to emphasize that the existence of a lower bound on travel time is significant, because it limits the 
speed with which one can perform unitary transformations dynamically. Such transformations play a central role in 
quantum computation. For example the construction of efficient NOT-gates involves unitary transformations that 
map a state into its antipodal state. The distance between antipodal states in 'P('Hphys) is the same as in V(H). 

Hence, for such states tJ 7 ^ = r m j n . 

The situation is quite different if we consider the evolution generated by the rj + -pseudo-Hcrmitian Hamiltonian 

H± + in the standard projective Hilbert space V(7i). In this case, we can indeed obtain arbitrarily fast evolutions, 
but they will not be unitary [16]. The possibility of infinitely fast non- unitary evolutions is actually not surprising. 
What is rather surprising is that one can achieve such evolutions using quasi-Hermitian Hamiltonians; there are 
arbitrarily fast quasi-unitary evolutions [171]. 

A scenario that is also considered in [41] is to use both Hermitian and quasi-Hcrmitian Hamiltonians to produce 
an arbitrarily fast evolution of into X^, F . This is done in three stages. First, one evolves the initial state X^, I into 
an auxiliary state X^ using a Hcrmitian Hamiltonian hi in time n, then one evolves X^ into another auxiliary state 
X^> F using a quasi-Hermitian Hamiltonian H in time r', and finally one evolves X^i into the desired final state A^ F 
using another Hermitian Hamiltonian /12 in time T2 . By choosing the intermediate states X^> and X^> appropriately 
one can make n and T2 as small as one wishes. By choosing H to be an rj + -pseudo-Hcrmitian operator of the form 
(376) with the parameter ki of 77, sufficiently small one can make the total travel time r := n + r' + T2 smaller than 
r m i n . In this scenario both the initial and final states belong to V(H), but to maintain unitarity of the evolution 
one is bound to switch (the defining metric of) the physical Hilbert space at t = t\ and t = t\ + r'. Therefore, this 
scheme involves a physical Hilbert space with a time-dependent inner product. As discussed in Subsection 7.1, the 
latter violates the condition that the Hamiltonian is an observable. Therefore, there seems to be no legitimate way 
of lowering the bound on travel time between two states of a given distance except allowing for non-unitary (possibly 
quasi-unitary) evolutions. 

9 Physical Applications 

Since its inception in the form of 'PT-symmctric models in the late 1990's and later as a consistent quantum mechanical 
scheme [175] pseudo-Hcrmitian QM has been the subject of extensive research. The vast majority of the publications 
on the subject deal with issues related to formalism or various (quantum mechanical as well as field theoretical) toy 
models with mostly obscure physical meaning. There are however a number of exceptions to this general situation 
where concrete problems are solved using the methods developed within the framework of pseudo-Hcrmitian QM. In 
this section we outline the basic ideas upon which these recent developments rest. Before engaging into a discussion 
of these, however, we wish to list some of the applications that predate the recent activities in the field. 

9.1 Earlier Applications 
9.1.1 Dyson Boson Mapping 

Among the earliest manifestations of pseudo-Hermitian operators is the one appearing in the context of the Dyson 
mapping of Hermitian Fermionic Hamiltonians to equivalent quasi-Hermitian bosonic Hamiltonians [81]. Dyson 
mapping has subsequently found applications in nuclear physics [115] and provided the basic idea for the formulation 
of quasi-Hermitian QM [207]. For a brief review of the Dyson mapping method see [92]. 
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9.1.2 Complex Scaling and Resonances 

Consider the one-parameter family of operators: ru a = exp (|?| {x,p}^j with a £ C, that act in the Hilbcrt space 
L 2 (R). We can easily use the Backer-Campbell-Hausdorff identity (150) together with the canonical computation 
relation [x,p] = ih to show that U Q induces a scaling of the position and momentum operators: x —> U a x U" 1 = e a x 
and p — > ll Q p U" 1 = e~ a p. For a g R, U Q is a unitary transformation, and one can use the latter property to show 
that (U a ip)(x) = e-^ip{e a x) for all ip g L 2 (R). 

For a g C — R, the transformation -0 — >• U Q '0 is called a complex scaling transformation. In this case, U„ is 
no longer a unitary operator. In fact, neither \l a nor its inverse is bounded. This implies that its action on a 
Hcrmitian Hamiltonian H, namely H —> H' = Ua-HU" 1 , that (neglecting the unboundedness of U Q and U^ 1 ) maps 
H into a quasi- Hcrmitian Hamiltonian H' , can have dramatic effects on the nature of its continuous spectrum. This 
observation has applications in the treatment of resonances (where one replaces x with the radial spherical coordinate 
in R 3 ). The main idea is to perform an appropriate complex scaling transformation so that the non-square-intcgrablc 
wave functions representing resonant states of H are mapped to square-integrable eigenfunctions of H' . For details, 
see [214, 133, 12] and references therein. 

9.1.3 Vortex Pinning in Superconductors 

Consider the Hamiltonian: H g = fe±M^ — \- v(x), where g is a real constant and v is a real- valued potential. This 
Hamiltonian can be mapped to the Hcrmitian Hamiltonian Hq by the similarity transformation: 

H g -> e~ 9X f n H g e ax ' n = H . (378) 

This in turn implies that H g is rj + -pseudo- Hermitian for the metric operator rj + := e~ 2 9 x / ft . This is to be expected, 
because (378) is an example of the quasi-Hcrmitian Hamiltonians of the form (200) that admit x-dependent metric 
operators. In [101, 102] the Hamiltonian H g (with random potential v) is used in modeling a derealization phe- 
nomenon relevant for the vortex pinning in superconductors. A review of the ensuing developments is provided in 
[100]. 

9.2 Relativistic QM, Quantum Cosmology, and QFT 

The issue of constructing an appropriate inner product for the defining Hilbert space of a quantum mechanical system, 
which we shall refer to as the Hilbert-space problem, is almost as old as quantum mechanics itself. Probably the first 
serious encounter with this problem is Dirac's attempts to obtain a probabilistic interpretation of the (first-quantized) 
Klein-Gordon fields in the late 1920's. The same problem arises in the study of other bosonic fields and particularly 
in the application of Dirac's method of constrained quantization for systems with first class constraints [78]. This 
method defines the "physical space" V of the state-vectors as the common null space (kernel) of the constraints, but 
it does not specify the inner product necessary to make V into a Hilbert space. Often the inner product induced from 
the auxiliary Hilbert space of the unconstrained system is not physically admissible, and one must find an alternative 
method of constructing an appropriate inner. 

In trying to deal with the Hilbert-space problem for Klein-Gordon fields, Dirac was led to the discovery of the wave 
equation for massive spin-half particles and the antimatter that earned him the 1933 Nobel prize in physics. Another 
major historical development that has its root in attempts to address the Hilbert-space problem is the discovery of 
the method of second quantization and eventually relativistic quantum field theories. These developments did not 
bring a definitive resolution for the original problem, but diminished the interest in its solution considerably. 

In the 1960's the discovery of the Hamiltonian formulation of the General Theory of Relativity [14] provided 
the necessary means to apply Dirac's method of constrained quantization to gravity. This led to the formulation of 
canonical quantum gravity and quantum cosmology [74, 232] and brought the Hilbert-space problem to forefront of 
research in fundamental theoretical physics for the second time. In this context it emerges as the problem of finding 
an appropriate inner product on the space of solutions of the Wheeler-DeWitt equation. Without such an inner 
product these solutions, that are often called the "wave functions of the universe," are void of a physical meaning. 
The lack of a satisfactory solution to this problem has been one of the major obstacles in transforming canonical 
quantum gravity and quantum cosmology into genuine physical theories [130, 107]. 

A widely used approach in dealing with the Hilbert-space problem for Klein-Gordon and Proca fields is to use 
the ideas of indefinite-metric quantum theories. These fields admit a conserved current density whose integral over 
space-like hypcrsurfaces yields a conserved scalar charge. This is however not positive-definite and as a result cannot 
be used to define a positive-definite inner product and make the space of all fields V into a genuine Hilbcrt space 
directly. This makes one pursue the following well-known scheme [230]. First, one uses the conserved charge to 
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define an indefinite inner product on V, and then restricts this indefinite inner product to the (so-called positive- 
energy) subspace of V where the inner product is positive-definite. The common practice is to label this subspace as 
"physical" and define the Hilbert space using the fields belonging to this "physical space." 

This approach is not quite satisfactory, because even for physical fields the above-mentioned conserved current 
density can take negative values [23]. Therefore it cannot be identified with a probability density. There are also 
other problems related with the observables that mix "physical fields" with "unphysical fields" or "ghosts." 

The application of pseudo-Hcrmitian QM in dealing with the Hilbert space problem in relativistic QM and 
quantum cosmology [150, 151, 161, 179, 114, 238], and the removal of ghosts in certain quantum field theories 
[35, 119] relies on the construction of an appropriate (positive-definite) inner product on the space of solutions of the 
relevant field equation. 79 

The basic idea behind the application of pseudo-Hermitian QM in dealing with the Hilbert-space problem in 
relativistic QM and quantum cosmology is that the relevant field equations whose solutions constitute the state- 
vectors of the desired quantum theory are second order differential equations in a "time" variable. 80 . These equations 
have the following general form. 

^(t)+£ty(t) = 0, (379) 

where t denotes a dimensionless time variable, ip : K — > C is a function taking values in some separable Hilbert space 
C, and D : C — > C is a positive-definite operator that may depend on t. 
We can express (379) as a two-component Schrodinger equation [87], 

ij^(t) = Hffi(t), (380) 
where * : R -> C 2 and H : C 2 -> C 2 are defined by [140, 141] 

U(*)-#) )' ~2{-D + l -D-l)> (381) 

£ 2 stands for the Hilbert space C © £, and a dot denotes a ^-derivative. The Hamiltonian (381) can be easily shown 
to be quasi- Hcrmitian [150]. 

In Subsection 3.5, we examined in detail the quantum system defined by the Hamiltonian (381) for the case that 
C is C with the usual Euclidean inner product and D is multiplication by a positive number. In this case, the field 
equation (379) is the classical equation of motion for a (complex) harmonic oscillator with frequency s/D. It turns 
out that most of the practical and conceptual difficulties of addressing the Hilbert-space problem for Klein-Gordon, 
Proca, and Wheeler-DeWitt fields can be reduced to and dealt with in the context of this simple oscillator. 81 In 
particular, the cases in which D is ^-dependent (that arises in quantum cosmological models) require a more careful 
examination. We will not deal with these cases here. Instead, we refer the interested reader to [151] where a 
comprehensive discussion of these issues and their ramifications is provided. 

Following the approach taken in Subsection 3.5, one can construct a metric operator rj + : C 2 — > C 2 and a new 
inner product on C 2 that renders H Hermitian. This defines a physical Hilbert space JC of two-component 

fields ^(t). Because H : JC — > JC is Hcrmitian, it generate a unitary time-evolution in JC. In particular, for every 
initial time to & every pair and ^ of solutions of the Schrodinger equation (380), and all t € M., 

(*i(t)|* 2 (t)>„ + = (*i(*o)|* 2 (*b)>, + . (382) 

As a vector space JC (and C 2 ) are isomorphic to the space of solutions of the single-component field equation (379), 
i.e., V := {ip : K — > C | ijj(t) + Dip(t) = for all (el}. We can obtain an explicit realization of this isomorphism 
as follows. Let to be an initial time and U to : V — > JC be defined by 

U t M) :=*(*o)- (383) 

According to (381) and (383), the effect of U to on solutions ip of the field equation (379) is to map them to the 
corresponding initial conditions tj}(to) and ip(to)- Because the field equation is linear and second order, this mapping 
is a linear bijection. Therefore, IAt is a vector space isomorphism. This is an important observation, because it 



79 This should be distinguished with the treatment of the Pais-Uhlenbeck oscillator proposed in [49], because the latter involves changing 
the boundary conditions on the field equation which in turn changes the vector space of fields. 

80 This is the physical time variable in an inertial frame in relativistic QM or a fictitious evolution parameter in quantum cosmology 
which may not be physically admissible [151]. 

81 For a discussion of this particular quantization of the classical harmonic oscillator, see [159]. 
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allows us to use U ta to induce a positive-definite inner product (• 
induced inner product is denned by 



(''Pulp: 



■) n+ on V form the inner product 



for all ipi,ip2 £ V. 



on K. The 



(384) 



Note that in view of (383) and (382), the right-hand side of (384) is independent of the value of to. This makes 
(•, -) r;+ into a well-defined inner product on V and gives it the structure of a Hilbert space. 82 

In principle different choices for r] + give rise to different inner products (•,•)?; > but the resulting Hilbert spaces 
H r)+ := (V, (v)»7 ) are unitary-equivalent. The arbitrariness in the choice of r/ + can be restricted by imposing 
additional physical conditions. For example in the case of Klein-Gordon and Proca fields, the requirement that 
(•, - ) n be Lorentz- invariant reduces the enormous freedom in the choice of rj + to a finite number of free numerical 
parameters [150, 238]. 

The most general Lorcntz-invariant and positive-definite inner product on the space of (real or complex) Klein- 
Gordon fields has the following form [179] 



(i>i,i>2)r,. 



ihK 
2mc 



da 1 " 



ipi(x)* d^C^x) +aipi(x)* <9 M ip2(x) 



(385) 



where £ is a spacelike Cauchy hypersurface in the Minkowski spacetime, x := (x , x , x 2 , x 3 ) are the spacetime 
coordinates in an inertial frame, ipi and ip2 are a pair of solutions of the Klein-Gordon equation: ti 2 [—Oq + V 2 ] ip(x) = 
m 2 c 2 ip(x), such that for all x° G R, ip(x°,x) and doip{x°, x) define square-integrable functions of x := (x 1 , x 2 , x 3 ), 



c?3, for any pair of differentiable functions / and g, f d^g := fd^g — gd^f, k G 



8,, := d/dx", V 2 := d 2 

and a G (—1,1) are arbitrary dimensionless free parameters demonstrating the arbitrariness in the choice of rj + , and 
C is the grading operator defined by 



(W)(x) 



h 2 



-1/2 



ip{x) 



dk 3 dy 3 



(386) 



Note that —V 2 + H 2 /(m 2 c 2 ) is a positive operator acting in C = L 2 (R 3 ) and that C is Lorentz-invariant [179]. 

According to (384), as a linear operator mapping H n+ to /C, Ut is a unitary operator. Similarly p := ^Jff^ is a 
unitary operator mapping K. to C 2 . Therefore pU v+ : TL V+ — > C 2 is also unitary. Usually C is an L 2 -space with well- 
known self-adjoint operators. This allows for a simple characterization of the self-adjoint operators o acting in C 2 . 
We can use these operators and the unitary operator pU v+ to construct the self-adjoint operators O : H v+ —> W. n+ 
that serve as the obscrvables of the desired quantum theory. This is done using 



o = ( P u„ Y 1 o P u n 



(387) 



The application of this construction for Klein-Gordon [161, 179] and Proca [238] fields yields explicit expressions for 
the corresponding relativistic position operators and localized states, a problem that has been a subject of ongoing 
research since the 1940's [193, 188]. A natural consequence of these developments is the construction of a set of 
genuine relativistic coherent states for Klein-Gordon fields interacting with a constant magnetic field [180]. 



9.3 Electromagnetic Wave Propagation 

An interesting application of pseudo-Hermitian QM is its role in dealing with the centuries-old problem of the 
propagation of electromagnetic waves in linear dielectric media [176]. Unlike the applications we discussed in the 
preceding section, here it is the spectral properties of quasi-Hcrmitian operators and their similarity to Hcrmitian 
operators that plays a key role. 

Consider the propagation of the electromagnetic waves inside a source-free dispcrsionless (linear) dielectric medium 
with dielectric and permeability tensors £ = £(x) and \i = \X[x) that may depend on space igF but not on time 
ieK. Maxwell's equations in such a medium read [111] 

V • D = 0, V ■ B = 0, (388) 

B + VxE = 0, D-VxH = 0, (389) 

82 Strictly speaking one must also perform a Cauchy completion of the inner product space obtained by endowing V with (■, -) v . 
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where E and B arc the electric and magnetic fields, a dot means a time-derivative, and 

D:=EE, H:=*[i B. (390) 

Eqs. (388) and (389) arc respectively called the constraint and dynamical equations. The former may be viewed 
as conditions on the initial values of the electromagnetic field, because once they are satisfied for some initial time 
the dynamical equations ensure their validity for all time. 

Similarly to Klein-Gordon equation, we can express the dynamical Maxwell equations (389) as first order ordinary 
differential equations for state-vectors belonging to a separable Hilbert space. To achieve this we introduce the 
complex vector space V of vector fields F : R 3 —> C 3 and endow it with the inner product (-F1I-F2) := L3 d xFx(x)* ■ 
F 2 (x), for all Fi,F 2 G V, to define the Hilbert space of square-integrable vector fields: H := {F : M 3 -> C 3 | (F\F) < 
00 }. The operation of computing the curl of the (differcntiablc) elements of this Hilbert space turns out to define 
a linear Hermitian operator ID : T~L — > H according to (T)F)(x) := V x F(x). We can use D to write (389) in the 

form: B(t) + DE(t) = and D(t) — DH(t) = 0. Evaluating the time-derivative of both sides of the second of these 
equations and using the first of these equations and (390), we find 

E{t) + n 2 E(t) = 0, (391) 

where tt 2 : H -> H is defined by tt 2 := s £ /! S. 

In view of the fact that £, fi, and consequently Q 2 are time- independent, we can integrate (391) to obtain the 
following formal solution 

E(t) = cos(D,t)E + ft -1 sm(nt)E , (392) 
where E := E(0), E := E(0) = ^D^L ^(0), and 

00 (_-\\n 00 (_\\n 

cos(m) := ]T \z±r (t 2 n 2 r, O- 1 sin(^) := t ^ j±-L- {t 2 n 2 T . (393) 



n=0 v ' n=0 



Given initial values of the electric and magnetic fields E(0) and B(0), we can use (392) and (393) to obtain a 
series expansion for the evolving electric field. One can select specific initial fields so that this expansion involves 
a finite number of nonzero terms, but these are of little physical significance. In general the resulting solution is 
an infinite derivative series expansion that is extremely difficult to sum or provide reliable estimates for. A crucial 
observation which nonetheless makes this expansion useful is that, for the cases that £ and jJi are Hermitian, the 

operator Q : T-L — > H. is £-pseudo-Hcrmitian: Q. = £ fl £ .In particular, for lossless material where £ is a 

positive x-dependent matrix, Q 2 is a quasi-Hcrmitian operator 83 that can be mapped to a Hermitian operator h by 

1 

a similarity transformation, namely h = p Q, 2 p^ 1 where p := £ . 

In terms of h the solution (392) takes the form: E(t) = p~ 1 [cos{h^t)p Eq + h~i sm(h^t)p Eq\. Therefore, 

E(x,t) = (x\E(t)) = p- l {x) I d 3 y \c(x,y;t)p(y)E Q (y)+ < S(x,y;t)p(y)E (y)} , (394) 

where ^ ^ 

C{x,y;t) := (x\cos(hh)\y), S{x, y; t) := (x\h~i sin(hh)\y). (395) 

<-> 

The fact that h is a Hermitian operator acting in % makes it possible to compute the kernels C {x, y; t) and 
S(x,y;t) of the operators cos(h^t) and h^^ sin(h^t) using the spectral representation of h: 

N 

n—l a 

where the sum over the spectral label n should be identified with an integral or a sum together with an integral 
whenever the spectrum of h has a continuous part, E n and \ip n .a) denote the eigenvalues and eigenvectors of h 
respectively, and a is a degeneracy label. In view of (396), for every analytic function F of h, such as cos(h^t) 



83 As an operator acting in H the dielectric tensor £ plays the role of a metric operator. This is one of the rare occasions where a 
metric operator has a concrete physical meaning. 
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and h a sin(/i2t), we have (x|F(/i)|y) = X)^Li S a F{E n ) ipn^ffiipn^iv)* ■ In the scattering setups where E and 
// tend to constant values, h has a continuous spectrum and one finds integral representations for the kernels (395) 
that reduce the solution (394) of Maxwell's equations into performing certain integrals (after solving the eigenvalue 
problem for h.) 

Rcf. [176] outlines the application of this method for the cases that the medium is isotropic, the initial fields as 
well as the dielectric and permeability constants change only along the ^-direction, and the WKB approximation is 
applicable in dealing with the eigenvalue problem for h. Under these conditions, one can compute the kernels (395) 
analytically. This allows for the derivation of the following closed form expression for the propagating electromagnetic 

field in terms of the initial fields E , E , and the z-dependent dielectric and permeability constants e(z) and fJ,(z). 



E(z,t) 



1 


'n(z)~ 


'{ 


e(w-(z,t)) 


2 


[<*). 




_fi(w-(z,t))_ 



E (w-(z,t)) + 



■(w+(z,t)) 



fx(w + (z,t)) 



E (w+{z,t)) 



w+ {z,t) 



W- (z : t) 



dw n(w) i e(w) 1 Eo(w) 



where w±(z,t) := u^ 1 (u(z) ± t) and u(z) :— J Q d% y/ 'e{l) , and u^ 1 stands for the inverse function for u, [176]. 

The possibility of the inclusion of dispersion effects in the above approach of solving Maxwell's equations is 
considered in [174]. 



9.4 Other Applications and Physical Manifestations 

The following arc some other areas where pseudo-Hcrmitian operators arise and/or the methods of pseudo-Hcrmitian 
QM arc used in dealing with specific physics problems. 



9.4.1 Atomic Physics and Quantum Optics 

Effective quasi-Hcrmitian scattering Hamiltonians arise in the study of the bound-state scattering from spherically 
symmetric short range potentials. As shown by Matzkin in [136], the use of the machinery of pseudo-Hcrmitian QM 
in the study of these Hamiltonians leads to a more reliable quantitative description of the scattering problem. It also 
provides a better understanding of the approximation schemes used in this context in the past and allows for their 
improvement. 

The relevance of pseudo-Hcrmitian operators to two-level atomic and optical systems has been noted in [27, 26, 
204], and their application in describing squeezed states is elucidated in [73, 25]. The optical systems provide an 
important arena for manufacturing non-Hcrmitian and in particular pseudo- and quasi-Hermitian effective Hamil- 
tonians. Recent experimental studies of 'PT-symmetric periodic potentials that make use of "PT-symmetric optical 
lattices is based on this observation [134]. See also [52]. 



9.4.2 Open Quantum Systems 

The emergence of non-Hcrmitian effective Hamiltonians in the description of the resonant states, that is based on 
Feshbach's projection scheme [86] is a very well-known phenomenon [181]. The application of a similar idea, that 
replaces the projection scheme with an averaging scheme, for open quantum systems also leads to a class of non- 
Hermitian effective Hamiltonians (usually called Liouvillian or Liouville's super operator) [57]. These Hamiltonians 
that determine the dynamics of the reduced density operators can, under certain conditions, be pseudo-Hcrmitian or 
even quasi-Hcrmitian. In [220, 112, 221], Stenholm and Jakob explore the application of the properties of pseudo- and 
quasi-Hermitian operators in the study of open quantum systems. The key development reported in these articles 
is the construction of a metric operator, that uses the spectral method we discussed in Subsection 4.1, and the 
identification of the corresponding norm with a viable candidate for a generalized notion of entropy. 



9.4.3 Magnetohydrodynamics 

Pseudo- Hermitian effective Hamiltonians arise in the study of the dynamo effect in magnetohydrodynamics [96, 97]. 
These Hamiltonians are typically non-quasi-Hcrmitian and involve exceptional points. Therefore, they can only be 
treated in the framework of indefinite-metric theories and using the properties of Krein spaces [18]. 
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9.4.4 Quantum Chaos and Statistical Mechanics 

In [71], Date el al study the spectrum of the Hamiltonian operator: H = | (p x + pr) 2 + \ (p y — ) 2 , where a 
is a real coupling constant and r := \/ x 1 + y 2 . This Hamiltonian that is Hcrmitian and 'PT-symmetric describes 
a rectangular Aharonov-Bohm billiard. Here Note that the spectrum is obtained by imposing Dirichlet boundary 
condition on the boundary of the rectangular configuration space, that is defined by |a;| < a and \y\ < b for some 
a,b E R + , and also at the location of the flux line, namely x = y = 0. The main result of [71] is that the nearest 
neighbor spacing distribution for this system has a transition that interpolates between the Poisson (level clustering) 
and Wigner (level repulsion) distributions. In an attempt to obtain a random matrix model with this kind of behavior, 
Ahmed and Jain constructed and studied certain pseudo-Hermitian random matrix models in [4, 5]. 



9.4.5 Biophysics 

In [84], Eslami-Moossallam and Ejtehadi have introduced the following effective Hamiltonian for the description of 
the dynamics of an anisotropic DNA molecule. 

H: =A + ^k + i +iuJoJ3 - f ' cos ^ (397) 

where Ax, A2, C, u)$, and / are real coupling constants, a, /?, 7 are Euler angles, and J\,J2,J3, that satisfy the com- 
mutation relations for angular momentum operators, arc defined by J± := —i ^— ^ + sin7 + cot /3cosj J^J , 

J2 := —i ^fj^a + cos 7 -§p ~ c °t /3sin7 J^j , and J3 := —ijpj- Clearly the Hamiltonian (397) is non- Hcrmitian, 
but it is at the same time real, i.e., it commutes with the time- reversal operator T . In light of the fact that T is an 
antilinear operator, this implies that H is a pseudo-Hermitian operator [145]. It would be interesting to see if this 
observation has any physically interesting implications, besides the restriction it puts on the spectrum of H. 



10 Summary and Conclusions 

In this article, we addressed various basic problems related to what we call pseudo-Hermitian quantum mechanics. 
Our starting point was the observation that a given quantum system admits an infinity of unitary-equivalent repre- 
sentations in terms of Hilbert space-Hamiltonian operator pairs. This freedom in the choice of representation can be 
as useful as gauge symmetries of elementary particle physics. 

Wc have surveyed a variety of mathematical concepts and tools to establish the foundations of pseudo-Hermitian 
quantum on a solid ground and to clarify the shortcomings of the treatment of the subject that is based on the so- 
called charge operator C. We showed that it is the metric operator 77, that plays the central role in pseudo-Hermitian 
quantum mechanics. Although, one can in general introduce a C operator and express r/ + in terms of C, the very 
construction of observables of the theory and the calculations of the physical quantities requires the knowledge of r? + . 
This motivates addressing the problem of the computation of a metric operator for a given quasi-Hermitian operator. 
Wc have described different approaches to this problem. 

We have discussed a number of basic issues related to the classical-to-quantum correspondence to elucidate the 
status of the classical limit of pseudo-Hermitian quantum mechanics. We have also elaborated on the surprising 
limitation on the choice of time-dependent quasi-Hermitian Hamiltonians, the role of the metric operator in path- 
integral formulation of the theory, a treatment of the systems defined on complex contours, and a careful study 
of the geometry of the space of states that seems to be indispensable for clarifying the potential application of 
quasi-Hermitian Hamiltonians in generating fast quantum evolutions. 

Finally, we provided a discussion of various known applications and manifestations of pseudo-Hermitian quantum 
mechanics. 

Among the subjects that we did not cover and suffice to provide a few references for are pseudo-supersymmetry 
and its extensions [146, 215, 203, 216], weak pseudo-Hcrmiticity [219, 20, 242, 164], and the generalizations of VT- 
symmetry [31, 170]. This omission was particularly because of our intention not to treat the results or methods with 
no direct or concrete implications for the development of pseudo-Hermitian quantum mechanics. We particularly 
avoided discussing purely formal results and speculative ideas. 
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Appendix: Reality of Expectation values Implies Hermiticity of the Ob- 
servables 

Theorem 3: Let TL be a Hilbert space with inner product (-|-) and A : % — > % be a (densely-defined, closed) linear 
operator satisfying T>(A) = T>(A'), i.e., A and its adjoint A^ have the same domain D(A). Then A is a Hermitian 
operator if and only if (ip\Aip) is real for all tp € T> (A). 

Proof: If A is Hermitian, we have (</)\Aip) = (A<p\ip) for all xj),<j> G D(A). Then according to property (ii) of 
Subsection 2.1, (tp\Aip) G R for all V G V(A). Next, suppose that for all ip € V(A), (ip\Aip) G R. We will show that 
this condition implies the Hermiticity of A in two steps. 

Step 1 : Let A + := + A^) and A_ := ±(A- A*). Then A = A + + iA_, V{A±) = V(A), and A± are Hermitian 
operators. In view of the first part of the theorem, this implies that 

(ip\A±ip) G R, for all V G V(A). (398) 

Furthermore, according to A = A + +iA^ and the hypothesis of the second part of the theorem, (tp\A + ip)+i(ip\A-ip) = 
(i>\Aip) G R. This relation and (398) show that 

(ip\A_tp) = for all tp € D(A). (399) 

5tep 2: Let </>, ^ be arbitrary elements of T>(A), £± := <p ± ^, and £± = ± i^. Then a direct calculation, using the 
property (Hi) of Subsection 2.1, shows that = | ((^ + |A_^ + ) - - t(C+|A_C+) +«(C-I^-C-)) = 0, 

where the last equality follows from (399) and the fact that £±,C± G T>(A). This establishes <f>\A-tp) = for all 
</),ip G 25(A). In particular setting = A-ip, we find = which in view of the property (i) of 

Subsection 2.1 implies A_ip = for all ip € D(A). Hence A_ = 0, and according to A = A + + iA-, we finally have 
A = A + . But A + is Hermitian. □ 
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